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1  INTRODUCTION 


An  intensive  program,  comprised  of  modeling  efforts  and  intensive  exper¬ 
imental  measurements,  was  conducted  over  the  past  three  years  aimed  at 
improving  our  understanding  of  the  radar  response  of  terrain  at  millimeter- 
wave  frequencies.  The  experimental  measurements,  which  were  conducted 
primarily  at  35,  94,  and  140  GHz,  included  laboratory  investigations  of  spe¬ 
cific  targets-such  as  small  trees-as  well  as  field  observations  of  trees  and 
snow-covered  terrain.  Of  particular  interest  was  to  model  and  character¬ 
ize  the  polarimetric  radar  backscatter  of  statistically  homogeneous  targets. 
To  this  end,  polarimetric  radar  systems  were  constructed  and  appropriate 
calibration  techniques  were  developed. 

This  report  provides  a  summary  of  the  results  realized  under  this  pro¬ 
gram,  with  the  details  provided  in  Appendix  A  in  the  form  of  reprints  of 
articles  published  in  scientific  journals  and  symposia  proceedings.  Not  in¬ 
cluded  in  the  Appendix  are  two  long  reports  [24,  25]  documenting  MM- Wave 
radar  observations  acquired  under  this  program,  as  well  as  data  reported 
by  other  institutions.  These  reports  have  been  requested  by  and  provided 
to  numerous  U.S.  Army  and  other  DOD  groups  and  laboratories.  The  ref¬ 
erences  cited  in  the  body  of  this  report  refer  to  the  publication  list  given  in 
Section  3. 

2  SUMMARY  OF  RESULTS 

2.1  MM- Wave  Polarimetric  Radar  Systems 

Under  a  previous  ARO  contract,  we  had  developed  radar  scatterometer 
systems  at  35,  94,  and  140  GHz.  In  order  to  investigate  the  polarimetric  re¬ 
sponse  of  terrain,  we  needed  to  modify  the  system  and  develop  appropriate 
calibration  techniques  to  allow  us  to  make  accurate  measurements  of  the 
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backscattered  signed  in  a  polarimetric  mode.  For  indoor  measurements  in 
an  anechoic  chamber,  the  fully  coherent  polarimetric  technique  can  be  used 
to  measure  the  scattering  matrix  of  the  target,  from  which  both  the  mag¬ 
nitude  and  phase  statistics  can  be  obtained  directly.  Under  outdoor  condi¬ 
tions,  however,  wind-caused  fluctuating  movements  of  the  radar  antenna  or 
movements  of  the  illuminated  targets  (such  as  leaves  and  branches  on  trees) 
can  lead  to  large  errors  in  the  polarization  phase  measurements  because  the 
spatial  variations  are  comparable  to  the  wavelength  at  millimeter-wave  fre¬ 
quencies.  Hence,  it  was  necessary  to  develop  an  alternate  measurement 
technique  that  is  not  vulnerable  to  spatial  movements  of  the  radar  antenna 
relative  to  the  illuminated  target.  The  coherent-on-receive  polarimetric 
measurement  technique  satisfies  this  condition.  To  evaluate  the  accuracy  of 
this  technique,  the  35-GHz  radar  system  was  configured  to  operate  in  both 
the  fully  coherent  mode  and  the  coherent-on-receive  mode,  and  then  it  was 
used  to  measure  a  target  placed  on  a  rotating  platform  in  both  modes. 
Having  verified  that  the  two  techniques  indeed  provide  the  same  Mueller 
matrix,  modifications  were  then  made  to  convert  the  94  and  140  GHz  sys¬ 
tems  into  coherent-on-  receive  polarimetric  configurations.  These  results 
are  discussed  in  [19,  21,  and  22]. 

2.2  Modeling  Polarimetric  Backscatter  From  Terrain 
2.2.1  Tree  Canopies 

A  radiative  transfer  model  was  developed  for  characterizing  the  polarimetric 
radar  response  of  tree  canopies  at  millimeter  wavelengths  [1].  The  model 
uses  a  phase  matrix  whose  functional  form  is  based  on  experimental  bistatic 
measurements  obtained  under  laboratory  conditions.  The  model  was  found 
to  be  in  very  good  agreement  with  experimental  data  acquired  under  field 
conditions.  One  of  the  major  conclusions  of  this  study  is  that  the  first-order 
solution  of  the  radiative  transfer  equation  provides  reasonable  agreement 
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with  data  for  the  like-polarized  scattering  coefficient,  but  not  for  the  cross- 
polarized  component.  However,  the  second  order  solution  provides  very 
good  agreement  for  both.  Related  contributions  on  scattering  from  foliage 
are  given  in  references  [6,  14,  21]. 

2.2.2  Snow-Covered  Terrain 

A  radiative  transfer  model  with  tthe  quasi-crystalline  approximation  was 
developed  for  snow-covered  terrain.  The  model  provided  very  good  agree¬ 
ment  with  experimental  observations  of  dry  snow  at  35,  94,  and  140  GHz. 
To  incorporate  the  dependence  on  the  wetness  profile  of  snow  under  wet 
snow  conditions,  a  hybrid  first-order/numerical  model  was  developed  and 
used  to  predict  the  temporal  variation  of  the  backscatter  from  snow.  The 
results  are  available  in  references  [3,  4,  15,  23]. 

2.2.3  Polarization  Phase  Statistics 

The  co-polarized  and  cross-polarized  phase  distributions  can  be  obtained 
from  multiple  measurements  of  the  scattering  matrix  of  the  terrain  under 
observation.  Such  measurements  are  provided  directly  by  fully  coherent  po- 
larimetric  systems.  However,  with  coherent-on-receive  systems,  the  quan¬ 
tity  measured  is  the  Mueller  matrix.  Also,  theoretical  models  cannot  predict 
the  scattering  matrix  but  can  provide  results  for  the  Mueller  matrix.  Until 
recently,  no  process  was  available  by  which  the  probability  density  func¬ 
tions  of  the  co-polarized  and  cross-polarized  phases  could  be  derived  from 
the  Mueller  matrix.  References  [7]  and  [8]  contain  derivations  of  new  rela¬ 
tionships  that  allow  the  calculation  of  the  probability  density  functions  in 
terms  of  the  elements  of  the  Mueller  matrix.  The  new  relationships  provide 
a  critical  link  between  theory  and  experimental  observations  of  the  polar¬ 
ized  phase  statistics,  making  it  possible  to  use  theoretical  models  to  relate 
the  phase  statistics  to  physical  parameters  of  the  terrain. 
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2.3  Millimeter- Wave  Radar  Data  Base 

The  data  handbooks  [24,  25]  of  millimeter- wave  radar  scattering  from  ter¬ 
rain  contain  a  compilation  of  data  measured  at  35,  94,  140,  and  215  GHz, 
as  reported  in  the  open  literature.  The  unique  features  of  these  handbooks 
are: 

(a)  They  axe  the  only  handbooks  available  for  millimeter-  wave  radar  data 
of  terrain. 

(b)  All  data  is  presented  using  a  uniform  format  throughout. 

(c)  Only  calibrated  data  that  has  been  augmented  with  ” ground-  truth” 
information  is  included. 

3  LIST  OF  PUBLICATIONS 

The  following  list  contains  papers  published  in  scientific  journals  or  sym¬ 
posia  proceedings,  or  presented  at  technical  symposia  and  workshops.  All 
of  these  papers  were  generated  in  support  of  this  project. 

Journal  Papers 

1.  Ulaby,  F.T.,  T.H.  Haddock,  and  Y.  Kuga,  “Measurement  and  Modeling 
of  Millimeter- Wave  Scattering  From  Tree  Foliage,”  Radio  Science ,  Vol. 
25,  pp.  193-  203,  1990. 

2.  Haddock,  T.F.  and  F.T.  Ulaby,  “140  GHz  Scatterometer  System  and 
Measurements  of  Terrain,”  IEEE  Trans,  on  Geoscience  and  Remote 
Sensing ,  Vol.  28,  pp.  492-499,  July,  1990. 

3.  Kuga,  Y.,  F.T.  Ulaby,  T.F.  Haddock  and  R.  DeRoo,  “Millimeter- Wave 
Radar  Scattering  From  Snow:  Part  I-Radiative  Transfer  Model  with 
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Quasi-Crystalline  Approximation”  Radio  Science,. V ol.  26,  pp.  329- 
341,  March  1991. 

4.  Ulaby,  F.T.,  T.  Haddock,  R.  Austin  and  Y.  Kuga,  “Millimeter- Wave 
Radar  Scattering  From  Snow:  Part  II-Comparison  of  Theory  with  Ex¬ 
perimental  Observations”  Radio  Science,. V ol.  26,  pp.  343-351,  March 
1991. 

5.  Sarabandi,  K.  and  F.T.  Ulaby,  “High  Frequency  Scattering  from  Corru¬ 
gated  Stratified  Cylinders”,  IEEE  Trans.  Antennas  and  Propagation, 
Vol.  39,  pp.  512-520,  April,  1991. 

6.  Sarabandi,  K.,  F.T.  Ulaby,  and  T.B.A.  Senior,  “Millimeter  Wave  Scat¬ 
tering  Model  for  a  Leaf”,  Radio  Science,  V.  25,  N.  1,  pp.  9-18,  Jan.-Feb. 
1990. 

7.  Sarabandi,  K.,  “Derivation  of  Phase  Statistics  From  the  Mueller  Ma¬ 
trix,”  accepted  for  publication  in  Radio  Science,  1992. 

8.  Ulaby,  F.T.,  K.  Sarabandi,  and  A.  Nashashibi  “Statistical  Properties 
of  the  Mueller  Matrix  of  Distributed  Targets”,  Special  Issue  of  IEE 
Proceedings  Part  F:  Remote  Sensing  Radar,  1992. 

9.  Sarabandi,  K.,  Y.  Oh,  and  F.T.  Ulaby,  “Measurement  and  Calibra¬ 
tion  of  Differential  Mueller  Matrix  of  Distributed  Targets” ,  submitted 
for  publication  in  IEEE  Transactions  on  Antennas  and  Propagation, 
February,  1992 


Symposia  Papers 

10.  Kuga,  Y.,  Austin,  R.  T.,  Haddock,  T.  F.,  and  F.  T.  Ulaby,  “Calculation 
of  the  Diurnal  Backscattering  Characteristics  of  Snow  at  35  and  95 
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GHz,”  Progress  in  Electromagnetics  Research  Symposium ,  July  25-26, 
1989,  Cambridge,  Massachusetts. 

11.  Haddock,  T.  F.  and  F.  T.  Ulaby,  “140  GHz  Scatterometer  Measure¬ 
ments,”  1989  International  Geoscience  and  Remote  Sensing  Sympo¬ 
sium  (IGARSS  \ 89 ),  July  10-14,  1989,  Vancouver,  Canada. 

12.  Sarabandi,  K.,  F.  T.  Ulaby,  T.  B.  A.  Senior,  “Millimeter  Wave  Scatter¬ 
ing  Model  for  a  Leaf,”  1989  International  Geoscience  and  Remote  Sens¬ 
ing  Symposium  ( IGARSS  ’ 89 ),  July  10-14,  1989,  Vancouver,  Canada. 

13.  Ulaby,  F.  T.,  T.  F.  Haddock,  and  Y.  Kuga,  “Measurements  and  Mod¬ 
eling  of  Millimeter- Wave  Scattering  from  Tree  Canopies,”  1989  Inter¬ 
national  Geoscience  and  Remote  Sensing  Symposium  (IGARSS  ’89), 
July  10-14,  1989,  Vancouver,  Canada. 

14.  Ulaby,  F.  T.,  “Millimeter- Wave  Bistatic  Scattering  by  Terrain,”  AGARD 
Conference  Proceedings ,  Kopenhagen,  Denmark,  Sept.  9-13,  1989. 

15.  Kuga,  Y.,  R.  T.  Austin,  T.  F.  Haddock,  and  F.  T.  Ulaby,  “Calcu¬ 
lation  of  the  Diurnal  Backseat tering  Characteristics  of  Snow  at  35 
and  95  GHz,”  1989  Progress  in  Electromagnetic  Research  Symposium 
(PIERS),  July  25-26,  1989,  Massachusetts  Institute  of  Technology, 
Cambridge,  Massachusetts. 

16.  Ulaby,  F.  T.  ,  Y.  Kuga,  and  T.  Haddock,  “Measurement  and  Modeling 
Millimeter- Wave  Scattering  from  Snow,  URSI  Specialist  Meeting  on 
Signature  Problems  in  Microwave  Remote  Sensing,  16-18  May,  1990, 
Hyannis,  Massachusetts. 

17.  Ulaby,  F.,  T.  Haddock,  and  Y.  Kuga,  “Radar  Clutter  Database  at 
Millimeter- Wave  Frequencies”,  Workshop  on  Detection ,  Discrimina¬ 
tion  and  Classification  of  Targets  in  Clutter,  13-15  November,  1990. 


18.  Kuga,  Y.,  F.T.  Ulaby,  and  T.F.  Haddock,  “Millimeter- Wave  Radar 
Scattering  from  Snow”,  Workshop  on  Detection,  Discrimination  and 
Classification  of  Targets  in  Clutter ,  13-15  November,  1990. 

19.  Kuga,  Y.,  K.  Sarabandi,  A.  Nashashibi,  and  F.  Ulaby,  “Millimeter 
Wave  Polarimetric  Scatterometer  Systems:  Measurement  and  Calibra¬ 
tion  Techniques”,  AGARD,  6-10  May,  1991,  Montreal,  Canada. 

20.  Ulaby,  F.,  and  Y.  Kuga,  “Modeling  and  Measuring  Millimeter- Wave 
Scattering  from  Snow-Covered  Terrain”,  AGARD ,  6-10  May,  1991, 
Montreal,  Canada. 

21.  Nashashibi,  A.,  Y.  Kuga,  and  F.T.  Ulaby,  “Polarimetric  Observations 
of  Trees  at  35  and  94  GHZ”,  1991  North  American  Radio  Science 
Meeting  and  International  IEEE/AP-S  Symposium ,  June  24-28,  1991, 
The  University  of  Western  Ontario,  London,  Ontario. 

22.  Kuga,  Y.,  A.  Nashashibi  and  F.T.  Ulaby,  “Clutter  Measurements  by 
Millimeter-  wave  Radars”,  National  Telesystems  Conference ,  Atlanta, 
GA,  26-27  March,  1991. 

23.  Kuga,  Y.,  A.  Nashashibi  and  F.T.  Ulaby,  “Millimeter- wave  Polarimet¬ 
ric  Radar  Scattering  from  Snow”,  IGARSS  ’ 91 ,  3-7  June,  1991,  Fin¬ 
land. 

Technical  Reports 

24.  Haddock,  T.F.,  and  F.T.  Ulaby,  “Millimeter  -  Wave  Radar  Scattering 
from  Terrain:  Data  Handbook  I”,  Radiation  Laboratory  Technical  Re¬ 
port  026247-T-l,  The  University  of  Michigan,  Ann  Arbor,  Michigan, 
January,  1990. 
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25.  Haddock,  T.F.,  and  F.T.  Ulaby,  “Millimeter  -  Wave  Radar  Scattering 
from  Terrain:  Data  Handbook  IF,  Radiation  Laboratory  Technical 
Report  026247-T-2,  The  University  of  Michigan,  Ann  Arbor,  Michigan, 
September,  1990. 

4  PARTICIPATING  SCIENTIFIC  PERSONNEL 

The  following  people  participated  in  the  MMW  radar  program: 

Facultv  and  Research  Scientists 

_  ** 

Professor  Fawwaz  T.  Ulaby 
Professor  Yasuo  Kuga 
Dr.  Tom  Haddock 
Dr.  Kamal  Sarabandi 

Graduate  Students 

Dr.  Kamal  Sarabandi,  received  Ph.D  (1989) 

Mr.  Adib  Nashashibi,  expected  PhD  completion  in  1993 

Mr.  Roger  DeRoo,  received  MS  (1990),  expected  PhD  completion  in  1993 

Mr.  Richard  Austin,  expected  PhD  completion  in  1992 

Mr.  Steve  Ciccarelli,  expected  MS  completion  in  1992 

Mr.  Paul  Siqueira,  expected  PhD  completion  in  1994 

Research  Engineer 
Mr.  Ron  Hartikka 
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5  CONCLUSIONS 


The  work  conducted  under  this  program  has  led  to  significant  contributions 
in  terms  of  characterizing  the  radar  backscatter  from  terrain  at  millimeter 
wavelengths,  particularly  with  regard  to  the  new  field  of  radar  polarimetry. 
A  major  strength  of  the  program  has  been  the  strong  connection  between 
theoretical  modeling  and  experimental  observations.  A  key  area  of  future 
research  is  to  explore  the  relationships  between  the  statistics  of  the  co¬ 
polarized  and  cross-polarized  phase  angles  of  the  backscattered  signal  and 
the  physical  properties  of  the  terrain.  The  phase  statistics  provide  a  new 
vector  of  observational  parameters  which  may  convey  significant  informa¬ 
tion  about  the  observed  scene.  This  area  of  research  is  in  its  infancy  and 
deserves  serious  exploration. 
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APPENDIX  A:  REPRINTS  OF  SELECTED 

PAPERS 

This  appendix  contains  reprints  of  a  select  set  of  papers  generated  under 
this  program.  These  are: 

A.l  Ulaby,  F.T.,  T.H.  Haddock,  and  Y.  Kuga,  “Measurement  and  Model¬ 
ing  of  Millimeter- Wave  Scattering  From  Tree  Foliage,”  Radio  Science, 
Vol.  25,  pp.  193-  203,  1990. 

A. 2  Haddock,  T.F.  and  F.T.  Ulaby,  “140  GHz  Scatterometer  System  and 
Measurements  of  Terrain,”  IEEE  Trans,  on  Geoscience  and  Remote 
Sensing ,  Vol.  28,  pp.  492-499,  July,  1990. 

A. 3  Kuga,  Y.,  F.T.  Ulaby,  T.F.  Haddock  and  R.  DeRoo,  “Millimeter- Wave 
Radar  Scattering  From  Snow:  Part  I-Radiative  Transfer  Model  with 
Quasi-Crystalline  Approximation”  Radio  Science, ol.  26,  pp.  329- 
341,  March  1991. 

A. 4  Ulaby,  F.T.,  T.  Haddock,  R.  Austin  and  Y.  Kuga,  “Millimeter- Wave 
Radar  Scattering  From  Snow:  Part  II-Comparison  of  Theory  with  Ex¬ 
perimental  Observations”  Radio  Science, .V ol.  26,  pp.  343-351,  March 
1991. 

A. 5  Sarabandi,  K.,  F.T.  Ulaby,  and  T.B.A.  Senior,  “Millimeter  Wave  Scat¬ 
tering  Model  for  a  Leaf’,  Radio  Science,  V.  25,  N.  1,  pp.  9-18,  Jan.-Feb. 
1990. 

A. 6  Kuga,  Y.,  K.  Sarabandi,  A.  Nashashibi,  and  F.  Ulaby,  “Millimeter 
Wave  Polarimetric  Scatterometer  Systems:  Measurement  and  Calibra¬ 
tion  Techniques”,  AGARD,  6-10  May,  1991,  Montreal,  Canada. 
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Because  the  constituent  elements  of  a  tree  canopy,  namely  the  leaves,  needles,  branches,  and 
trunks,  have  complex  geometries  with  curvatures  and  surface  roughness  scales  that  are  comparable 
to  or  larger  than  the  wavelength  at  millimeter  wavelengths,  the  traditional  approach  used  to  compute 
the  phase  function  of  the  vegetation  volume  is  totally  impractical.  In  this  paper  we  propose  a 
relatively  simple  model  for  characterizing  the  phase  function  on  the  basis  of  direct  experimental  data. 
The  model  is  used  in  conjunction  with  a  solution  of  the  radiative  transfer  equation  to  predict  the 
backscattering  behavior  of  tree  canopies.  The  model  is  found  to  provide  very  good  agreement  with 
radar  observations  made  at  35.  94.  and  140  GHz. 


I.  INTRODUCTION 

The  leaves,  needles,  branches,  and  trunks  com¬ 
prising  a  vegetation  canopy  are  lossy  dielectric 
structures  with  complex  geometries.  Whereas  it 
may  be  acceptable  to  approximate  a  leaf  as  a  thin, 
flat  disc  at  centimeter  and  longer  wavelengths,  such 
a  treatment  is  invalid  at  millimeter  wavelengths 
because  the  leaf  curvature  and  its  thickness  are 
comparable  to  or  larger  than  A.  Similar  statements 
can  be  made  with  regard  to  the  size  and  surface 
roughness  of  branches  and  other  components  of  a 
vegetation  plant  or  tree.  Consequently,  it  is  very 
difficult,  if  not  impossible,  to  use  numerical  quadra¬ 
ture  techniques  for  solving  the  vector  radiative 
transfer  equation  [Ulaby  et  al.,  1986;  Tsang  et  al., 
1985]  to  compute  the  radar  backscattering  coeffi¬ 
cient  of  vegetation  at  millimeter  wavelengths.  The 
problems  encountered  are  that  (1)  it  is  not  possible 
to  compute  the  scattering  phase  function  of  the 
vegetation  volume  because  accurate  models  for  the 
scattering  matrices  of  the  scattering  elements 
(curved  leaves,  rough-surface  branches,  etc.)  are 
not  available  at  millimeter  wavelengths,  and  even  if 
such  models  were  available,  the  numerical  compu¬ 
tations  that  would  have  to  be  performed  to  obtain 
the  phase  function  (which  involve  integration  over 
size  and  orientation  parameters)  would  be  ex¬ 
tremely  expensive,  and  (2)  when  the  phase  function 
has  a  complicated  dependence  on  the  bistatic  scat- 
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tering  angles,  it  is  very  difficult  to  compute  the 
solution  of  the  radiative  transfer  problem  beyond 
the  first  order.  Hence  an  alternate  approach  is 
needed  for  computing  the  radar  backscatter  from 
vegetation  at  millimeter  wavelengths. 

In  this  paper  we  shall  propose  a  relatively  simple 
model  for  characterizing  the  phase  function  of  veg¬ 
etation  canopies  at  millimeter  wavelengths,  and 
then  use  it  in  a  second-order  solution  of  the  radia¬ 
tive  transfer  equation  to  compute  radar  backscat¬ 
tering  from  trees.  Model  results  are  compared  with 
radar  backscatter  measurements  for  tree  canopies 
at  35,  94,  and  140  GHz. 

2.  PHASE  MATRIX 

Except  for  the  main  trunk,  tree  foliage  consists 
primarily  of  randomly  distributed  leaves  (or  nee¬ 
dles)  and  branches,  most  of  which  are  much  larger 
than  the  wavelength  in  size  (at  millimeter  wave¬ 
lengths).  have  complex  shapes,  and  are  character¬ 
ized  by  quasi-uniform  orientation  distributions. 
These  properties  suggest  that  whereas  individual 
scattering  elements  may  exhibit  highly  complex  and 
polarization-dependent  scattering  patterns,  an  ele¬ 
mental  volume  dV  containing  many  of  these  ele¬ 
ments  is  likely  to  exhibit  propagation  and  scattering 
properties  that  are  weakly  polarization-dependent 
and  characterized  by  a  relatively  simple  scattering 
pattern.  This  prediction  is  supported  by  experimen¬ 
tal  observations  made  by  Ulaby  et  al.  [1988]  which 
show  that  bistatic  scattering  from  trees  exhibits 
comparable  results  for  horizontal  and  vertical  po¬ 
larizations. 
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The  random  nature  of  the  tree  foliage  supports 
the  use  of  radiative  transfer  theory  [Ishimaru,  1978] 
for  modeling  millimeter  wave  propagation  in  the 
canopy  [Ulaby  et  al. ,  1988;  Schwering  et  al. ,  1988]. 
In  the  radiative  transfer  model  the  formulation  is 
given  in  terms  of  the  phase  function  P(0S,  <£,;  0(,  <t>j) 
relating  the  specific  intensity  scattered  by  a  unit 
volume  of  the  scattering  medium  into  the  direction 
(0S,  to  the  specific  intensity  incident  upon  the 
unit  volume  from  the  direction  (0,,  <£,),  with  both 
sets  of  orientation  angles  being  defined  with  respect 
to  a  reference  coordinate  (x,  y,  i).  The  phase  matrix 
represents  the  average  Stokes  matrix  of  the  parti¬ 
cles  constituting  the  unit  volume.  To  relate  P  to  the 
properties  of  the  medium,  we  start  by  considering 
scattering  by  a  single  particle.  For  a  plane  wave 
with  electric  field  vector  E'  incident  upon  the  par¬ 
ticle  in  the  direction  k,  =  (0,,  <£,),  the  far-field  wave 
scattered  by  the  particle  in  the  direction  k5  =  (0S, 
<t>s)  is  a  spherical  wave  with  field  vector  E5 .  The 
vertical  and  horizontal  polarization  components  of 
E5  at  a  range  r  from  the  scatterer  are  related  to  the 
components  of  the  incident  field  through  the  scat¬ 
tering  matrix  S(0s,  <J>S;  0,,  <t>,)  of  the  particle  [Ulaby 
et  al. ,  1986], 


^(k,,  k,)  =  [4irr2  •  |£;|2/l£i|2]  =  4n\ Svh\2  (3) 

When  considering  an  elementary  volume  containing 
N  randomly  distributed  particles  per  unit  volume, 
we  can  characterize  bistatic  scattering  by  the  vol¬ 
ume  in  terms  of  the  bistatic  scattering  cross  section 
per  unit  volume  (or  bistatic  scattering  coefficient) 

«vh  -  % *> 

=  4irN<|S„/l|2>  (4) 

where  angle  brackets  denotes  ensemble  average. 

The  vector  radiative  transfer  equation  is  formu¬ 
lated  in  terms  of  the  specific  intensity  vector  I 
defined  through  the  modified  Stokes  parameters  /,., 
If,,  U,  and  V as  follows: 
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For  an  elementary  volume  of  length  ds  illuminated, 
in  the  general  case,  from  all  directions  k,  by  incident 
intensity  I'(k,),  the  intensity  scattered  in  the  direc¬ 
tion  k5  is  given  by 
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The  matrix  S  is  given  by  four  scattering  amplitudes, 
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$hv  $hh 


I(fe,)=/l  dS¥ik”kl)V{iL‘)da' 


(2) 


where  P  is  the  phase  matrix  given  by 
P(k„  fe*)  = 


^*11  P 12  P 13  P\* 
P 2i  P-n  Pn  Pi * 
P3\  P32  P33  P3 4 

P 41  Pi 42  P 43  P 44 


(6) 


For  a  specified  scattered/incident  polarization  com-  whose  elements  are  related  to  those  of  S  by 


r<|5n,|2) 

2<Re  (Sm,S;v)) 
2(Im  (SM,S:„)) 


(|S.*I2> 

(\S>*\2) 

2<Re  (5, *SJU)> 
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(Re  (5„.5S,)> 

(Re  (5Al,5&,)) 

(Re 

(Im  (5„,5^  +  S^SZ,,)) 


-(Im  (S,^)) 

-<Im  (S^Sjfc)) 

-(Im  -S,j,SZv)) 

(Re  (SwSZh  -  S^SU) 


(7) 


bination,  the  bistatic  scattering  cross  section  of  the 
particle  is  defined  in  terms  of  the  ratio  of  the 
scattered  to  incident  power  densities.  Fot  vh  polar¬ 
ization,  for  example. 


3.  PROPOSED  PHASE  MATRIX 

The  16  elements  of  P  can  be  readily  computed 
provided  we  know  (1)  the  number  density  N,  (2)  the 
probability  density  functions  for  the  sizes,  shapes. 
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and  orientations  of  the  particles,  and  (3)  the  dielec¬ 
tric  properties  of  the  particles,  and  additionally,  we 
have  available  appropriate  models  for  computing 
the  scattering  matrices  of  the  particles.  In  most 
cases,  this  information  is  not  available  for  terrain 
surfaces  and  volumes,  which  forces  investigators  to 
estimate  the  physical  parameters  of  the  canopy  and 
to  treat  the  canopy  constituents  as  spheres,  cylin¬ 
ders,  and  discs.  These  approximations  lead  to  errors, 
and  because  the  number  of  parameters  involved  is 
large,  it  is  difficult  to  assess  the  sources  of  error. 

3.1.  Phase  matrix  in  the  scattering  plane 

Instead  of  using  the  first  approach  described 
above  to  compute  the  elements  of  the  phase  matrix 
P.  we  propose  to  use  a  semiempirical  approach 
based  on  experimental  measurements.  Ulaby  et  al. 
[1988]  used  a  35-GHz  bistatic  radar  system  to 
examine  the  scattering  patterns  of  small  trees  under 
laboratory  conditions.  Two  types  of  experiments 
were  conducted:  (I)  transmission  measurements  to 
determine  the  extinction  coefficient  Kf  for  horizon¬ 
tal  and  vertical  polarizations  and  (2)  bistatic  scat¬ 
tering  measurements  in  the  plane  of  scattering  (de- 


Fig.  2.  Sketches  of  the  tree  architectures  and  photography  of  a 
Ficus  leaf  and  an  Arbor  Vitae  branch. 


V  =  270° 


6  a 

Fig.  I.  Configuration  used  for  measuring  bistatic  scattering 
from  tree  foliage.  The  tree  was  placed  on  a  rotating  platform,  the 
transmitter  was  in  a  fixed  location,  and  the  receiver  could  be  set 
at  any  angle  </». 


fined  to  be  the  plane  containing  the  incident  and 
scattered  directions  and  orthogonal  to  the  polariza¬ 
tion  planes  of  the  waves  [ Chandrasekhar .  I960])  to 
evaluate  the  angular  variations  of  the  like-  and 
cross-polarized  bistatic  scattering  cross  sections 
per  unit  volume  k!T(iM,  km,(<M.  and  Khh{tb), 

where  is  the  angle  shown  in  Figure  1 .  While  the 
transmitter  remained  in  one  location  with  the  beam 
pointing  at  the  crown  section  of  the  tree,  the  re¬ 
ceiver  was  moved  in  discrete  steps  around  a  circle 
in  the  horizontal  plane  with  the  tree  at  its  center.  At 
each  receiver  location,  defined  by  the  angle  <£,  the 
average  received  power  was  measured  and  then 
used  to  compute  «,„,(««,  kvHW),  kHv(\1 »).  and  Khh(<lt). 
The  averaging  process  was  realized  by  placing  the 
tree  on  a  rotating  platform  and  measuring  the  re¬ 
ceived  power  (for  a  given  receiver  location)  as  the 
tree  was  rotated  over  360°. 

Two  distinctly  different  types  of  trees  were  se¬ 
lected  for  examination:  Ficus  and  Arbor  Vitae  (Fig¬ 
ure  2).  The  Ficus  tree  had  small,  fiat,  simple  leaves 
approximately  10  cm2  in  area,  whereas  the  Arbor 
Vitae  tree  had  a  branching  trunk  arrangement  with 
branches  supporting  needles  approximately  1.5-3 
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FICUS  TREE 


Fig.  3.  Comparison  of  measured  bistatic  scattering  cross 
section  per  unit  volume  for  Ficus  tree  foliage  with  calculations 
based  on  the  model  functions  given  by  (14)  and  (15). 


ARBOR  VITAE  TREE 


Fig.  4.  Comparison  of  measured  bistatic  scattering  cross 
section  per  unit  volume  for  Arbor  Vitae  tree  foliage  with  calcu¬ 
lations  based  on  the  model  functions  given  by  (14)  and  (15). 


mm  in  length.  More  detailed  information  about 
these  test  trees  and  the  measurement  procedure  is 
given  in  Ulaby  et  al.  LI 988]. 

The  major  conclusions  derived  from  the  experimen¬ 
tal  observations  that  pertain  to  the  present  study  are: 

1.  For  both  types  of  trees,  the  like-polarized 
scattering  patterns,  ip)  and  Khh( ip),  were  approx¬ 

imately  the  same,  and  a  similar  result  was  observed 
for  Kvh( tf/)  and  Khv(tp).  Thus 


A 


XhhW~Klv{  I/O  =  K\(lp) 

(8) 

A 

KhvW  “  KvhW  =  X2W 

(9) 

2.  In  spite  of  the  fact  that  the  two  trees  were 
markedly  different  in  terms  of  the  shapes  and  sizes 
of  their  scattering  elements  (leaves,  needles, 
branches),  both  exhibited  similar  scattering  pat¬ 
terns.  Figures  3  and  4  shows  plots  of  the  measured 
values  of  K\(ip)  and  k2(i/>)  for  the  two  types  of  trees. 
Also  shown  are  plots  calculated  using  the  expres¬ 
sions  discussed  below. 

The  scattering  coefficient  k *  for  a  /i-polarized 
incident  intensity  is  given  by 

«S  =  4^  jjf  [«**&)  +  dtls  (10) 

and  a  similar  expression  can  be  defined  for  k's'.  In 
view  of  (8)  and  (9)  we  shall  set  =  k]'  =  kx.  If  we 
treat  foliage  as  an  isotropic  medium  and  if  we 


assume  azimuthal  symmetry  with  respect  to  the 
forward  scattering  direction  (ip  =  0),  the  like-  and 
cross-polarized  scattering  coefficients  can  be  ex¬ 
pressed  as 


«|(«A)  =  KiS'lW  (11) 

K2(4>)  =  KsSlW  (12) 

and  to  satisfy  (10),  the  sum  of<j|(i/0  andp2(*W  has  to 
satisfy  the  relation 


1 

2 


[9 1  («A)  +  ff2  W]  Sin  <p  =  1 


(13) 


In  view  of  the  shapes  of  the  measured  patterns 
(Figures  3  and  4),  g\(tp)  and  gi(ip)  can  each  be 
described  as  the  sum  of  a  relatively  weak  isotropic 
component  and  a  Gaussian-shaped,  strong  and  nar¬ 
row  forward-scattering  lobe  f(ip), 

9\W  =  L>i/iW  +  (1  —  at  1  )]C 


i 2  V 

/  <A\2 

exp 

.  U) . 

+  d  -«|) 

<14) 


92W  =  [« 2/2 ( •A)  +  (1  -  a2)](l  ~  C) 


+  (1  -  an) 


(I  -C) 


(15) 
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where  fit  and  /32  are  the  effective  beam  widths  of  the  pn  =  N(Rc  (S„,S^)) 

like-  and  cross-polarized  forward-scattering  lobes. 

These  expressions  have  the  following  properties:  =  N( Re  (|S1,.|e'*"|SlJe  ->*'*)) 


1 

2 


sin  i li  dip  =  C 


g2(ilf)  sin  ib  di/i  =  1  -  C 


(16) 


(17) 


L 


It 


afitli)  sin  i p  dtli 


(1  -  a)  sin  ill  dili 


=  «/( 1  -  a) 


(18) 


The  sum  of  properties  (16)  and  (17)  satisfies  (13),  the 
ratio  (1  -  C)/C  represents  the  ratio  of  total  scat¬ 
tered  cross-polarized  energy  to  total  scattered  like- 
polarized  energy,  and  the  ratio  (1  -  a)la  represents 
the  ratio  contained  in  the  isotropic  component  to 
the  energy  contained  in  the  main  lobe. 

The  “calculated”  plots  shown  in  the  Figures  3 
and  4  are  based  on  (11),  (12),  (14),  and  (IS),  with  the 
values  of  the  parameters  selected  to  provide  good 
agreement  between  the  measured  and  calculated 
plots. 

Now  let  us  return  to  the  phase  matrix  given  by 
(7).  The  element  P M  is  given  by 


/>„(*)  =  /V(|5„.|2) 


4ir 


1 


(19) 


Similarly,  it  is  easy  to  show  that 


=  NflSJlSrflcos  -  *,*)> 


=  NflS.jSjXcos  (*„.-  <»,*)>  (21) 

where  <f>tn,  is  the  phase  of  the  scattering  amplitude 
5ta,  (and  similar  definitions  apply  for  the  other 
scattering  amplitudes).  In  the  last  step  of  (21)  it  was 
assumed  that  the  magnitude  |5,J  and  the  phase 
difference  (<£,„,  -  <f>vh)  are  independent  random 
variables.  According  to  35-GHz  radar  measure¬ 
ments  of  the  backscattering  from  rocks  [Whitt  and 
Ulaby,  1988]  and  1.25-GHz  poiarimatric  data  ex¬ 
tracted  from  airborne  radar  images  of  forested  ar¬ 
eas,  the  phase  difference  (<f>w  -  <t>vh)  is  uniformly 
distributed  over  [0, 2v].  Hence  the  average  value  of 
cos  ~  4>vh)  *s  zero,  and  therefore  P)3  =  0. 
Similarly,  all  terms  in  P  involving  the  product  of  a 
like-polarized  scattering  amplitude  and  cross-polar¬ 
ized  scattering  amplitude  may  be  set  equal  to  zero. 

It  was  also  observed  (in  the  same  investigations) 
cited  above  that  the  phase  difference  -  d>hh) 
corresponding  to  the  product  of  the  like-polarized 
scattering  amplitudes  has  a  Gaussian-like  distribu¬ 
tion  centered  at  0°.  We  shall  therefore  adopt  the 
approximations  (cos  {d>m  -  4>hh ))  —  1  and  (sin  (<f>„. 
-  4>hh ))  =  0-  Furthermore,  in  view  of  (8)  we  shall 
assume  that  <|Sj  |S**I>  »  <|Sj2>.  Hence  for  the 
terms  involving  (Sm,SXh)  we  have 

N( Re  (51T5fc)>  =  /V<|S,„||S**|cos  <*,„  -  <t>hh )) 


=  JXcos  (*„,  -  *»*)> 

=  N(|S„,|2>  =  P  \\  (22) 


and 


N(\m  (S„.Sjf*))  =  0  (23) 


PziW)  —  Pu(<fr) 


and 


Pm(iIi)  =  =  (20) 

Next,  we  shall  make  certain  assumptions  to  sim¬ 
plify  the  remaining  terms.  Let  us  consider  the  term 
P 13  in  P, 


Upon  incorporating  the  preceding  results  in  (7),  we 
obtain  the  simplified  matrix 


PW  =  y- 
4ir 


9 1 

92 
0 
0 


92 

9\ 

0 

0 


0 

0 

9\  +92 
0 


0 

0 

0 

9\  ~92 


(24) 


with  g |  and  g2  as  given  by  (14)  and  (15),  respec¬ 
tively. 
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Fig.  5.  Coordinate  system  showing  the  incident  and  scattered 
angles*  (0, ,  d>! )  and  (05,  <£,),  and  their  relation  to  </*. 


3.2.  Phase  matrix  for  any  incident 
and  scattered  directions 

The  phase  matrix  P(tlt)  in  (24)  is  defined  in  terms 
of  the  scattering  angle  t{t  as  shown  in  Figure  1.  If  the 
incident  direction  is  (0,,  </>,)  and  the  scattered  direc¬ 
tion  is  (0S,  <t>s),  denoted  by  points  P\{0iy  <t>t)  and 


P <t>s\  0„  6,)  =  Un~y2)P(tlt)U-yt)  (25) 
where  L  is  a  linear  transformation  given  by 

cos’  y  sin2  y  2  sin  2y  o 

sin2  y  cos2y  -l  sin  2y  0 

y  -sin2y  sin  2y  cos  2y  0  ) 

0  0  o  1 

The  angles  yt  and  y2  are  defined  in  Figure  5. 
Introducing  the  abbreviations, 

(/,  /)  =  cos  yi  cos  y2  -  sin  y)  sin  y;  (27a) 

(/,  r)  =  -cos  yi  sin  y2  +  sin  yi  cos  y2  (21b) 

(r,l)~  sin  y!  cos  y;  +  cos  yj  sin  y2  (27c) 

(r,  r)  =  cos  yi  cos  y;  +  sin  yi  sin  y;  (21  d) 

we  can  write  the  phase  matrix  as 


P(0S,  <t>s;  0„  <*,)  =  -; 


gdMU)2 +g2mi,  r)2 
9\(<l>)(r ,  l)2  +  92(<li)(ry  r)2 
2lgd<l>)(l,  /)(r,  /)  +  g2(t}t)(r,  r)(l,  r)] 


gtmry  I)2  +  g2mr,  r)2 
g\mu)2+92mi,  r)2 
-  2[fl|(«/»)(r,  /){/,  /)  +  02(<frXr,  r)(l,  r)] 


-0i(</O(r,  /)(/,  /)  +  g2(>i>)(r,  r)(l,  r) 

5 UiMr,  /)(/,  /)  -  g2(<ii)(r,  r)(l ,  r) 

01  (*)((/.  I)2  ~  (r,  02]  +  0:(i/')[('\  r)2  ~  V*  ''I2! 


9\(<l>)  -  g2(>k) 


P2 (0S,  <t>s)  in  the  polar  coordinate  system  shown  in 
Figure  5,  the  plane  of  scattering  contains  the  trian¬ 
gle  OP\P2  and  the  scattering  angle  t^is  given  by  the 
angle  P{OP2.  The  radiative  transfer  equation,  on 
the  other  hand,  is  written  in  terms  of  the  polar 
angles  0  and  <t>.  We  need  to  obtain  a  new  phase 
matrix  in  terms  of  0  and  <t>  in  order  to  use  it  in  the 
radiative  transfer  equation.  The  details  of  the  trans¬ 
formation  from  P('lf)  to  P(0 5,  4>s\  0„  <i>, ).  which 
involves  linear  transformations  through  angles  yi 
and  ir  -  y2,  are  described  in  Chandrasekhar  [I960). 
The  transformed  phase  matrix  is  given  by 


Using  the  cosine  and  sine  laws  of  a  spherical 
triangle,  we  can  write  (l,  /),  (r,  r),  (r,  /),  and  (/,  r)  in 
terms  of  0S ,  <t>s,  0,  and  <f>,, 

sin2  (4>.  -  <t> .) 

(/,/)  = - r—i — : - sin  0,  sin  0, 

sin*  w 

•  (cos  tit  -  1)  -  cos  (4>,  -  <t>s)  (29a) 

sin2  (4>i  -  <t>s)  .  „ 

(r,  r)  = - r-j— ; - Sin  6S 

sur  tit 


sin  6 ,  (cos  \Jt  +  1)  -  cos  (<t>,  -  4>s) 
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r.  /)  =  | 


(sec  6;  +  sec  6S) 


sin  ( 4>i  -  cf>s) 


1  + 


sin  0S  sin  8;  cos  ip  cos  (<£,-  -  <ps) 
sin"  tp 


Sin  ((f);  —  <f>J  .. 

- — - [sin2  0,-  sec  0,  +  sin2  0S  sec  0j 

sirr  ip 


(29  c) 


/,  r)  =  | 


(/,  r)  -  -j  {sec  0,-  -  sec  6S) 


sin  (<t>j  -  4>s ) 


1  + 


sin  9S  sin  dj  cos  ip  cos  (<pj  —  <ps ) 


sin2  t/» 


sin  (<£,-  -  <£,)  r  .  2 


— i— — -  [sin*  0j  sec  0 
sin"  i/> 


,  -  sin2  05  sec  0r]  J 


(29d) 

where 

cos  tp  =  cos  0j  cos  dj  +  sin  0S  sin  8 ;  cos  (<£,■  -  <ps) 

(30) 


4.  RADIATIVE  TRANSFER  MODEL 

At  millimeter  wavelengths  the  penetration  depth 
of  foliage  rarely  exceeds  1  m.  Hence  it  would  be 
reasonable  to  neglect  the  backscatter  contribution 
of  the  underlying  ground  surface,  and  in  the  case  of 
most  vegetation  canopies  it  may  also  be  possible  to 
treat  the  canopy  as  semi-infinite  in  depth.  In  this 
section  we  seek  an  expression  for  the  vector  spe¬ 
cific  intensity  I*  scattered  from  a  forest  canopy 
characterized  by  a  phase  matrix  of  the  form  given 
by  (28).  To  this  end,  we  shall  develop  a  first-order 
solution  and  a  second-order  solution  of  the  radiative 
transfer  equation  and  then  compare  the  results  with 
the  exact  solution  (based  on  numerical  computa¬ 
tions  using  the  quadrature  method)  and  with  exper¬ 
imental  data. 

The  vegetation  canopy  is  modeled  as  a  continu¬ 
ous,  statistically  homogeneous,  horizontal  layer  of 
vertical  extent  d.  The  layer  has  diffuse  upper  and 
lower  boundaries  (Figure  6)  at  z  =  0  and  z  =  —  d, 
respectively,  and  it  is  illuminated  by  an  intensity 

I*  =  Io  5(cos  6  -  cos  60)S(<p  ~  <Po)  (31) 

incident  upon  the  upper  boundary  in  the  direction 
(it  —  6q ,  </>o).  Upon  solving  the  radiative  transfer 


Forest  Canopy 


I* 


Ground 

Fig.  6.  Geometry  of  the  scattering  problem. 


equation  to  obtain  an  expression  for  the  intensity 
<f>s)  scattered  in  any  direction  (0S,  <j>s),  we  can 
set  6S  =  0O  and  <ps  =  i t  +  4>0,  which  corresponds  to 
scattering  in  the  backward  direction,  to  compute 
the  backscattering  coefficient  from  the  equation 

o  ,  4tt  cos  80ISP(8 o,  it  +  (p0) 

<(0o)  = - P~  ,  . -  (32) 

Iq{TT  ~  00- 

where  p,  q  =  v  or  h  polarization. 

4.1.  Radiative  transfer  equations 

When  formulating  the  radiative  transfer  problem 
for  bounded  media,  the  standard  practice  is  to  split 
the  intensity  vector  into  upward-going  (I  +  (0i,  <f>s, 
z ))  and  downward-going  (I~(ir  —  9S,  <f>s,  z))  compo¬ 
nents,  noting  that  6S  varies  between  0  and  v/2 
[Ulaby  et  al.,  1986].  In  the  vegetation  layer  the 
intensity  I  +  (0J5  4>s,  z)  traveling  in  the  upward 
direction  (0S,  <f>s)  and  the  intensity  I“(tt  —  ds,  <ps,  z) 
traveling  in  the  downward  direction  (v  -  6S,  <f>s) 
must  satisfy  the  coupled  radiative  transfer  equations 

—  I  +  (ps,  <ps ,  z)  =  -—I  +  (/Xj,  (ps ,  z) 


+  F  +  (mj,  <t>s,  z) 


(33  a) 


d  Ke 

~7~  ^  (  —  Mr*  I  (  fix*  z) 

dz  ns 


+  F  (-p.s,<ps,z) 


(33  b) 
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where  Ke  is  the  extinction  coefficient  of  the  vegeta¬ 
tion  medium,  fxs  =  cos  0S,  and  -fxs  =  cos  (tt  -  9S). 
The  source  functions  F  +  (m,i  <f>s,  z)  and  F 
cf>s,  z )  account  for  directing  the  energy  incident 
upon  an  elemental  volume  from  all  directions  into 
the  direction  ( 0S ,  <f>s)  and  (v  -  6S,  respectively, 
and  are  given  by 

F  +  (m^  <t>s>  z) 

=  —  f  f  p (fj-s,  4>s\  Mm  +  (m;,  <i>i .  z)  dSli 

[Jo  Jo 

“t"  J  I  P(Mj»  $s<  Mm  $i')I  (  Mm  <£m  z)  dfl[- 

Jo  Jo 


(3  4a) 

F ~(-ns,  4>s .  z) 


1  r in  n 

-  I  P(  — Mil  4>s\  Ph  «#>,)!  (Mo  <f>i>  z)  dfl; 

Jo  Jo 


+  f  f  p (-Mi.  “Ml.  <*/)I  (~Mm  <t>i’  z)  d  ft; 

Jo  Jo 


I  (-p.f,  0J,  0)  =  I05(p,j  - -  (Ao)  (36a) 

I  +  (Mu  <t>s<  ~d)  -  0  (36*) 

and  because  F+  and  F“  are  themselves  integral 
functions  of  I+  and  I-,  we  have  to  use  numerical 
techniques  involving  segmentation  in  z  and  (9,  <j>)  in 
order  to  obtain  an  exact  solution  for  I*  =  I+(0*,  d>s, 
z  =  0).  While  this  may  be  useful,  particularly  for 
comparing  with  results  based  on  approximate  solu¬ 
tions,  the  numerical  technique  does  not  provide 
much  insight  with  regard  to  the  relative  importance 
of  various  scattering  contributions.  Hence  we  shall 
use  the  iterative  technique  to  develop  expressions 
for  the  first-order  and  second-order  solutions  of 
(35a)  and  (35*)  and  then  compare  their  results  with 
the  exact  results  of  the  numerical  solution.  The 
assumptions  underlying  the  iterative  technique  is 
that  the  medium  is  weakly  scattering,  i.e.,  the 
scattering  albedo  <o  =  Ks/Ke  <s:  1.  At  millimeter 
wavelength,  a>  —  0.6-0. 9  for  vegetation  [3],  and 
therefore  the  condition  is  not  satisfied.  Nonethe¬ 
less,  we  shall  now  proceed  with  the  iterative  tech¬ 
nique  and  then  evaluate  its  usefulness  in  a  later 
section. 

4.2.  First-order  solution 


(34*) 

where  dfl,  =  d^d^  =  s*n  and  P(Mj> 

M/>  <*,)  is  the  phase  matrix  of  the  vegetation  layer 
relating  the  intensity  incident  (upon  a  unit  volume  in 
the  medium)  in  the  direction  (0,  ,  <*,)  to  the  intensity 
scattered  in  the  direction  (8S,  <(>s).  The  phase  matrix 
is  defined  by  (28). 

The  solution  to  differential  equations  (33a)  and 
(33*)  can  formally  be  expressed  as 

I+(ms,  4>s,z)  =  e-K'iz  +  dv‘i‘lUns,  4>s,  ~d) 

+  J'  +  (Mj,  4>s,  z')  dz'  (35a) 

I  ~ (~Mj>  4>s,  z)  =  -(-fis,  4>s,  0) 

+  <t>5,  z')  dz'  (35*) 

Because  there  is  no  reflection  at  the  (diffuse) 
air-vegetation  boundaries  at  z  =  0  and  z  =  ~d,  the 
following  boundary  conditions  must  be  satisfied: 


We  start  with  the  zeroth-order  solutions,  which 
are  obtained  by  setting  P  =  0  in  (34a)  and  (34*), 
which  renders  F0+  =  F0“  =0  in  (35a)  and  (35*), 
where  the  zero  subscript  denotes  zero  order.  Using 
the  boundary  conditions  given  by  (34)  and  (35),  the 
zeroth-order  specific  intensities  are  given  by 

I0+(m,.  ^,z)  =  0  (37a) 

I0“  (“Mj>  4>s,  z) 

=  e*'^%8{ns  -  pom*,  ~  <£o)  (37*) 

The  zeroth-order  solution  corresponds  to  propa¬ 
gation  of  the  coherent  wave  through  the  medium 
with  scattering  ignored,  except  for  its  contribution 
to  extinction.  To  obtain  the  first-order  solution,  we 
first  need  to  insert  (37a)  and  (37*)  into  (34a)  and 
(34*)  to  compute  the  first-order  source  function  F  j1- 
and  Ff  and  then  insert  the  results  in  (35a)  and 
(35*).  This  process  leads  to 

e-K,Zfn, 

(Mji  4*  s  * 

Ps*  I 

<f>s,  -Mo^ollo 


(38) 
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If  ("Mi-  <*>f.  Z)  =  e«'^°l08(v.s  -  Mo Wd>s  -  <t> o) 

+ - [1  -e«3l]P(-M„  -Mo.  4>o)I0  (39) 

Ms*2 

where 

K|  =  Ke(l/fl0  +  lifts)  (40) 

A 

k2  =  Ke(l/fi0  -  l/ps)  (41) 

The  first-order  solution  for  the  backscattering  coef¬ 
ficient  can  be  obtained  by  setting  z  =  0,  ps  =  p0, 
and  4>s  =  7r  +  d>0  in  (38)  and  then  inserting  the  result 
in  (32).  These  steps  lead  to  the  expression 


<,(Oo)  = 


4v  cos  6q 
2  Ke 


[1 


_  ^  —2Kfd  sec  0o 


] 


*  [P(Mo.  tt  +  ^>o<  “Mo.  <£o)]u 


For  a  thick  canopy  such  that  ( 2Ked  sec  0O)  »  1 ,  the 
term  in  the  second  square  bracket  in  (42)  and  (43) 
reduces  to  1. 

4.3.  Second-order  solution 

The  second-order  solution  for  the  backscattered 
intensity  at  the  surface  I2+(0O,  tt  +  d>o.  0).  can  t>e 
obtained  by  (1)  replacing  (0S,  <j>s ,  z)  in  (38)  and  (39) 
with  ((?,-,  4>i,  z'),  (2)  inserting  the  resultant  expres¬ 
sions  in  (34a)  and  (34 b)  to  obtain  F f  (j is ,  <f>s,  z')  and 
Ff(— /u,,,  <f>s,  z'),  (3)  inserting  those  expressions  in 
(35a),  and  finally  (4)  replacing  (ps,  <f>s,  z)  with  (mo. 
t r  +  4> 0,  0).  This  process  leads  to 


1 2  (0o>  77  +  <£0,  0) 


F|+(mo.  tt  +  <t> o.  z’)  dz‘ 


(44) 


Ks  COS  0O  »  . 

=  [1  -e-2K'dMc9']g,W  (42) 

where  (P)  1 1  is  the  1 1  element  of  P  and  51(77)  is  given 
by  (14)  with  tp  =  tt. 

Similarly,  the  other  principal-polarization  back- 
scattering  coefficients  are  given  by 

o-°hh(0  0)  =  O^flo) 


«t(*o)  =  o^o) 


with 

F|+(M0i  7T  +  (p0,  z') 


=  —  [eK,z /,1"P(mo.  77  +  4>0;  -mo.  4>o)Io 
Mo 


nt  e-K'z/ti> 

- [*«,*  (Mo,  7r 

M,^3 


Kj  COS  0O 
2  Kf 


[1  -e-2't'dsec8%2(77) 


(43) 


Albedo  co 

Fig.  7.  Backscattering  coefficient  as  a  function  of  albedo  at 
normal  incidence,  computed  with  the  first-order,  second-order, 
and  numerical  solutions. 


+  4>o;  M/>  <£/)P(Mi.  -Mo.  4>o)Io  dft,- 

-2ir  ri  e*'*7**' 


ni 

M' 


k4 


[1  -eKU]P(Mo.  77 


+  0O;  “M/.  0;)P(-Mi.  4»/;  -mo.  0o)io  dfi/] 

(45) 


where 


*3 


(46) 


(47) 


Using  (14),  (15),  and  (28)  to  define  P,  the  integrals 
can  be  evaluated  numerically,  and  the  computed 
intensity  can  be  inserted  in  (32)  to  compute  cr®q 
for  any  p,  q,  =  v  or  h  polarizations. 
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Fig.  8.  Variation  of  backscattcring  coefficient  with  optica! 
thickness  t. 


(a)  35  GHz 


5.  MODEL  BEHAVIOR  AND  EXPERIMENTAL 
OBSERVATIONS 

Using  the  phase  matrix  given  by  (28)  with  the 
parameters  measured  for  the  Ficus  tree,  the  co¬ 
polarized  (co-pol)  and  cross-polarized  (x-pol)  back- 
scattering  coefficients  were  computed  for  a  variety 
of  canopy  conditions  in  accordance  with  (I)  the 
first-order  solution  of  section  4.2.  (a)  the  second- 
order  solution  of  section  4.3.  and  (3)  the  exact 
numerical  solution  using  the  quadrature  gradient 
technique  [Ulaby  et  al ..  1986].  Figure  7  shows  the 
variation  of  cr°  with  the  albedo  w  =  Ks/Ke .  for  a 
canopy  with  an  optical  thickness  t  =  Ked  =  I  Np. 
For  an  error  within  1  dB  of  the  numerical  solution 
the  first-order  solution  is  useful  up  to  &>  =  0.4  for  the 
co-pol  component  but  is  not  at  all  useful  for  the 
x-pol  component,  and  the  second-order  solution  is 
useful  up  to  a>  =  0.85  for  the  co-pol  component  but 
only  useful  up  to  w  =  0.5  for  the  x-pol  component. 
If  we  relax  the  error  margin  to  2  dB  for  the  x-pol 
component,  the  useful  range  of  w  may  be  extended 
up  to  0.85  for  the  second-order  solution. 

The  dependence  on  optical  thickness  is  illustrated 
in  Figure  8  for  all  three  solutions.  For  all  intents  and 
purposes,  cr°  is  independent  of  t  for  r  2:  1  Np.  This 
condition  is  almost  always  satisfied  for  tree  cano¬ 
pies  at  millimeter  wavelengths. 

Comparison  of  the  model  behavior  with  experi¬ 
mental  data  is  provided  in  Figure  9  which  shows 
measurements  of  (7°  as  a  function  of  incidence  angle 
for  a  canopy  of  Spruce  trees  at  35  GHz  and  a 
canopy  of  Bur  Oak  trees  at  94  and  140  GHz.  The 
model  calculations,  based  on  the  second-order  so¬ 
lution.  are  in  good  agreement  with  the  experimental 
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Fig.  9.  Comparison  of  theory  with  experimental  observa* 
tions  (<?)  35  GHz  for  a  canopy  of  Spruce  trees  and  at  ib)  94  GHz 
and  (c )  140  GHz  for  a  canopy  of  Bur  Oak  trees. 
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data  at  all  incidence  angles  and  frequencies  except 
for  the  70°  case  at  140  GHz  (Figure  9c).  No  obvious 
explanation  is  available  except  for  the  fact  that  oo  = 
0.95,  which  exceeds  the  region  of  validity  of  the 
second-order  solution. 

The  canopies  had  continuous  crown  sections,  the 
trees  were  about  10  m  in  height,  and  the  leaves  had 
a  moisture  content  of  53%  in  the  case  of  the  Spruce 
trees  and  27%  for  the  Bur  Oak  trees.  The  compu¬ 
tations  are  based  on  the  second-order  solution  using 
the  Ficus  phase-function  model  shown  in  Figure  3. 
The  only  free  parameter  used  in  attempting  to 
match  the  model  results  with  the  data  is  the  albedo 
(o,  which  was  chosen  to  be  equal  to  0.6  at  35  GHz, 
0.8  at  94  GHz,  and  0.95  at  140  GHz.  Similar  results 
were  obtained  in  attempting  to  match  the  model  to 
experimental  observations  for  (horizontally)  contin¬ 
uous  canopies  comprised  of  other  types  of  trees. 
This  observation  is  not  surprising  in  view  of  the 
strong  similarity  noted  earlier  between  the  bistatic 
scattering  patterns  shown  in  Figures  3  and  4  for  two 
trees  with  very  dissimilar  tree  architectures.  In 
other  words,  the  proposed  model  appears  to  apply 
to  a  wide  range  of  tree  types  of  continuous-crown 
canopies,  with  the  only  major  parameter  controling 
the  levels  of  the  co-pol  and  x-pol  backscattering 
responses  being  the  albedo  oo.  In  turn,  w  is  strongly 
dependent  on  the  wavelength  and  probably  depen¬ 
dent  on  leaf  moisture  content.  Further  study  is 
needed  to  establish  the  dependence  of  w  on  these 
two  parameters. 


6.  CONCLUSIONS 

Using  the  phase  matrix  model  proposed  in  this 
study,  radiative  transfer  theory  appears  to  provide 


excellent  agreement  with  experimental  observa¬ 
tions  of  the  backscatter  from  tree  canopies  at  35,  94. 
and  140  GHz.  The  only  free  parameter  used  in 
matching  the  model  to  data  is  the  scattering  albedo 
oo  which  appears  to  depend  on  only  two  parameters, 
the  wave  frequency  and  the  leaf  moisture  content. 
The  roles  of  shape  and  size  of  the  tree  leaves  or 
needles  and  the  tree  branch  architecture  appear  to 
be  secondary  in  importance.  Further  study  is 
needed  to  establish  the  exact  dependence  of  w  on 
moisture  content  and  frequency. 
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Abstract— The  goal  of  the  University  of  Michigan  millimeter-nave 
radar  program  is  to  characterize  terrain  scattering  at  35.  94,  and  140 
GHz.  The  140-GHz  channel  of  a  truck-mounted  scatterometer  system 
has  recently  been  added  to  give  the  full  desired  operating  capability. 
Two  injection-locked  45.33-GHz  Gunn  oscillators  use  triplers  to  supply 
the  up-  and  down-converters.  Full  polarization  capability  is  obtained 
through  the  use  of  rotatable  quarter-wave  plates.  Real-time  signal  pro¬ 
cessing  and  data  reduction  takes  place  in  an  HP  8510A  automatic  net¬ 
work  analyzer  on  the  truck-mounted  platform.  Sample  measurements 
of  millimeter-wave  radar  backscattering  from  vegetation  and  snow  are 
given. 

I.  Introduction 

MILLIMETER-WAVE  systems  offer  the  inherent  ad¬ 
vantages  of  high  resolution,  large  bandwidth,  and 
small  antenna  size.  In  recent  years,  significant  advances 
have  been  realized  in  the  development  of  millimeter-wave 
components  for  the  construction  and  operation  of  imaging 
airborne  radar  systems  at  the  atmospheric  window  fre¬ 
quencies  of  35,  94,  140,  and  215  GHz.  Hence,  there  is 
much  interest  in  measuring  terrain  scattering  at  these  fre¬ 
quencies  and  in  the  associated  development  of  theoretical 
and  empirical  scattering  models.  While  recent  work  has 
been  carried  out  at  millimeter  wavelengths  on  trees  [I] 
and  snow  [2],  [3],  such  data  are  still  sparse,  particularly 
at  140  GHz. 

The  University  of  Michigan  140-GHz  scatterometer 
system  is  the  latest  addition  to  the  network-analyzer-based 
millimeter-wave  scatterometer  system,  a  truck-mounted 
full-polarization  scatterometer  developed  in  support  of  a 
program  to  characterize  radar  scattering  from  terrain  at 
35,  94,  and  140  GHz.  The  basic  operation  of  the  scatter¬ 
ometer  system  has  been  described  in  detail  elsewhere  [4], 
[5].  Conversion  from  a  swept  2-4-GHz  intermediate  fre¬ 
quency  (IF)  to  the  millimeter-wave  RF  frequencies  is 
made  in  the  front  end,  allowing  flexible  real-time  signal 
processing  by  the  remotely  located  HP  85I0A  automatic 
network  analyzer.  An  outline  of  the  system  is  given  in 
Fig.  1.  An  HP  8350B  sweep  oscillator  is  swept  from  2  to 
4  GHz  by  the  HP  8510A  network  analyzer.  After  a  por¬ 
tion  of  this  IF  signal  is  taken  off  and  fed  to  the  a}  refer- 
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ence  port  of  the  HP  8511  frequency  converter,  it  is  up- 
converted  within  the  140-GHz  radar  and  transmitted  to 
the  target.  The  returned  signal  is  down-converted  to  the 
2-4-GHz  range  and  fed  to  the  b{  port  of  the  frequency 
converter.  Signal  processing  of  the  return  and  reference 
signals  takes  place  within  the  network  analyzer  and  is  sent 
on  the  HPIB  bus  to  the  HP  9920S  computer  where  final 
data  reduction  takes  place,  and  the  results  are  printed  out 
and  saved  on  disk.  The  HP  8510  computer-control  system 
allows  vector  error  correction  of  system  imperfections 
through  its  calibration  algorithms.  The  system  has  previ¬ 
ously  operated  in  this  mode  at  35  and  94  GHz,  and  the 
140-GHz  channel  is  its  latest  extension  in  frequency  ca¬ 
pability. 

II.  140-GHz  Scatterometer  Design 

A  block  diagram  of  the  140-GHz  front  end  is  shown  in 
Fig.  2.  The  transmit  portion  across  the  top  and  the  receive 
portion  across  the  bottom  are  driven  by  a  common  local 
oscillator  (LO)  chain.  The  LO  consists  of  a  45.33-GHz 
free-running  Gunn  oscillator,  two  circulator-coupled 
45.33-GHz  injection-locked  Gunn  oscillators  acting  as 
amplifiers  and  two  third-harmonic  frequency  multipliers. 
This  combination  provides  a  nominal  output  power  of  10 
dB  (1  mW)  from  each  multiplier  to  power  the  up-  and 
down-converters.  Other  combinations  of  fundamental  os¬ 
cillators,  amplifiers  or  frequency  multipliers  are  possible. 
This  particular  combination  provided  the  best  combina¬ 
tion  of  performance  and  cost.  Wave-polarization  is  con¬ 
trolled  by  a  fixed  quarter-wave  plate  followed  by  a  rotat¬ 
able  quarter-wave  plate.  A  90°  rotation  of  the  movable 
wave-plate  moves  the  electric  field  vector  through  90°  to 
give  either  vertical  or  horizontal  polarization.  The  polar¬ 
ized  signal  is  transmitted  though  a  conical  standard-gain 
horn  with  a  half-power  beam  width  of  1 1.8°.  The  received 
RF  signal  passes  though  a  3.0-in  diameter  lens-corrected 
horn  antenna  with  a  half-power  beamwidth  of  2.2°.  Hence 
the  antennas’  product  pattern  is  essentially  controlled  by 
the  receiver-antenna  pattern,  resulting  in  an  effective 
beamwidth  of  2, 1  ° .  Receive  polarization  is  determined  by 
movable  and  fixed  quarter- wave  plates,  in  the  same  man¬ 
ner  as  the  transmit  section.  The  RF  signal  is  down-con¬ 
verted  using  a  tripled  45.33-GHz  LO.  Since  the  conver¬ 
sion  processes  must  be  phase-coherent,  the  up-  and  down- 
converter  LO's  are  each  injection-locked  from  a  central 
dual-ended  Gunn  oscillator  running  at  45.33  GHz.  This 
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arrangement  gives  a  phase-coherent  LO  of  sufficient 
power  to  supply  both  up-  and  down-converters. 

III.  Calibration  and  Performance 

For  each  data  set.  measurement  of  a  sphere  of  known 
size  and  range  is  used  to  generate  the  401  VV  and  HH 
calibration  constants  for  each  of  the  40!  frequencies  in 
the  2-4-GHz  IF  band.  A  calibration  target  with  known 
cross-polarization  response  is  used  for  VH  and  HV  cali¬ 
bration.  For  an  incident  signal  consisting  of  either  pure 
vertical  (or  pure  horizontal)  polarization,  a  return  signal 
oriented  at  45°  to  vertical  is  generated  by  a  calibrator  con¬ 
sisting  of  a  rectangular  standard-gain  hom  followed  by  a 
38. 1-cm-long  section  of  WR6  waveguide  with  a  short  on 
the  end.  This  calibrator  is  placed  in  the  far  field  of  the 
140-GHz  radar,  and  pointed  toward  the  radar  with  the 
rectangular  aperture  of  the  hom  oriented  at  45°  to  hori¬ 
zontal.  While  a  portion  of  the  incident  radiation  is  re¬ 
flected  from  the  hom.  another  portion  passes  from  the  hom 
into  the  waveguide  and  propagates  in  the  TE)0  mode  with 
the  electric  field  in  line  with  the  short  axis  of  the  wave¬ 
guide.  which  is  oriented  at  45°  to  the  horizontal.  This 
signal,  comprising  equal  amplitude  vertical  and  horizon¬ 
tal  components,  is  reflected  by  the  short,  and  returns  to 
the  radar.  It  can  be  distinguished  from  the  return  from  the 
hom  aperture  by  its  longer  time  delay.  The  waveguide 
and  short  are  encased  in  a  metal  cylinder  to  prevent  return 
from  the  outside  of  the  guide  and  flange  at  the  range  of 
the  short.  Fig.  3  illustrates  the  W.  and  HH,  VH,  and  HV 
responses  of  the  cross-polarization  calibration  target. 
While  the  return  from  the  hom  is  complex,  the  return  from 
the  short  gives  a  known  cross-polarization  response.  At 
the  range  of  the  hom  in  Fig.  3.  the  like-  and  cross-polar¬ 
ized  responses  differ  by  approximately  5  dB.  At  the  range 
of  the  short,  where  for  a  perfect  radar  all  four  responses 
would  be  the  same,  all  returns  fall  within  a  ±  1  dB  range 
(within  experimental  uncertainties).  Measurement  of  this 
signal  is  used  to  generate  the  401  cross-polarization  cali¬ 
bration  constants. 


Fig.  2.  140-GHz  scaitcromeicr  front  end. 


Sphere  calibration  is  made  on  a  daily  basis,  but  the 
standard-gain  hom  cross-polarization  calibration  is  more 
cumbersome  and  is  made  less  frequently.  Cross-polari¬ 
zation  isolation  of  the  system  is  typically  about  15  dB. 
and  this  is  checked  at  each  use  of  the  system  by  making 
cross-polarization  measurements  of  the  sphere.  For  most 
natural  targets,  the  cross-polarized  return  at  140  GHz  is 
only  3-6  dB  below  the  like-polarized  return.  Hence  the 
cross-polarization  isolation  of  the  system  is  quite  ade¬ 
quate  at  140  GHz.  Noise  performance  of  the  system  is 
checked  after  each  calibration  by  making  measurements 
of  the  sky  at  typical  target  ranges.  Table  I  lists  the  mea¬ 
sured  system  performance  parameters. 

Fig.  4  illustrates  the  combined  effects  of  stability  and 
repeatability  of  calibration  of  the  140-GHz  system  over  a 
diurnal  cycle.  Repeatability  of  sphere  measufements  due 
to  pointing  only  is  typically  within  ±  0.5  dB.  Variations 
are  the  cumulative  result  of  system  gain  variations  and 
sphere  pointing  errors.  Installation  of  a  controlled  heater 
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Fig.  4.  System  calibration  constant  over  diurnal  cycle,  February  27,  28,  1989  (starting  10  a.m.). 


TABLE  1 

Parameters  of  Truck-Mounted  140-GHz  Scatterometer 


RF  frequencies 

138-140  GHz 

Transmit  power 

~4dB(l  mW) 

RF  bandwidth 

0-2.0  GHz 

Sweep  rate 

1  ms/freq.,  51,  101,  201,  401  freq/sweep 

Polarization 

VV%  VH%  HVt  HH 

Product-gain  beamwidth 

2.1° 

Incidence  angles 

o 

o 

1 

O 

o 

Platform  height 

2.7  m  minimum  to  18  m,  maximum 

Noise  equivalent  cr° 

- 20  dB 

Stability 

-0.2  dB/h 

Measurement  repeatability 

~1  dB 

Near-field  distance 

2.7  m 

Footprint 

minimum:  0.013  m2 
maximum:  16.8  m2 

Signal  processing 

HP  8510A/851 1A  based 

Output  products 

received  power  versus  range  (A R  =  c/2B  ) 
received  power  versus  frequency  (at  fixed  R  ) 
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Fig.  5.  Measured  backscatter  angular  response  of  pigweed  (Amaranthus  retroflexus)  grass  over  Stellaria  media  at  140  GHz. 


on  the  triple  LO  unit  was  required  to  achieve  the  excellent 
system  stability  observed  in  Fig.  4  ( ±  0.8  dB). 

For  most  terrain  measurements,  a  data  set  consists  of 
measurements  of  the  backscattering  coefficient  a0  as  a 
function  of  incidence  angle  for  W9  HH ,  and  VH  (or  HV) 
polarizations.  The  incidence  angle  is  set  by  an  elevation 
positioner  located  at  the  top  of  the  truck-mounted  boom. 
The  target  is  scanned  in  azimuth  to  obtain  spatially  inde¬ 
pendent  samples.  For  each  polarization  configuration  the 
number  of  independent  samples,  including  bandwidth  av¬ 
eraging,  is  at  least  50,  which  corresponds  to  a  measure¬ 
ment  precision  of  about  ±0.66  dB  [6].  Data  are  tabulated 
as  they  are  recorded  and  examined  in  real  time. 

IV.  Sample  Results 

Several  types  of  terrain  surfaces  and  covers  were  ob¬ 
served  by  the  140-GHz  scatterometer  in  1988  and  1989. 
Sample  results  are  shown  next  for  grasses,  trees,  and 
snow. 

A.  Backscatter  from  Grasses 

The  backscatter  plots  shown  in  Figs.  5  and  6  corre¬ 
spond  to  a  field  of  Amaranthus  retroflexus ,  a  spiny  weed 
about  50  cm  tall,  commonly  known  as  pigweed,  over 
ground  cover  of  Stellaria  media ,  a  low  ground-hugging 
weed,  commonly  known  as  chick  weed.  Fig.  5  shows  the 
140-GHz  backscatter  response  as  a  function  of  incidence 
angle  (measured  relative  to  normal  incidence)  for  W,  HH , 
and  HV.  Throughout  this  paper  the  “receive-transmit” 
convention  is  used.  As  expected  for  such  a  medium,  vol¬ 
ume  scattering  effects  predominate,  and  the  VV  and  HH 
returns  are  comparable  to  one  another  at  all  incidence  an¬ 
gles.  The  cross-polarized  return  is  approximately  6  dB 
lower  than  the  like-polarized  returns. 

Fig.  6  shows  the  HH  backscatter  response  from  the 
same  target  at  all  three  of  the  scatterometer  operating  fre¬ 
quencies:  35,  94,  and  140  GHz.  The  target  shows  a  weak 
sensitivity  to  frequency,  exhibiting  a  maximum  spread  of 
5  dB  between  the  three  curves. 


B.  Backscatter  from  Trees 

Fig.  7  shows  140-GHz  measurements  of  the  backscat¬ 
tering  coefficient,  plotted  as  a  function  of  incidence  angle, 
for  a  uniform  tree  canopy  of  Thuja  occidentalism  com¬ 
monly  known  as  white  cedar.  The  trees  were  approxi¬ 
mately  10  m  in  height,  and  the  average  water  content  of 
the  needles  was  measured  to  be  56.3%.  The  like-polari¬ 
zation  components  (HH  and  W)  are  essentially  identical 
in  level  and  exhibit  an  approximately  cos  6  dependence 
between  20°  and  70°.  The  HV  component,  on  the  other 
hand,  increases  with  increasing  incidence  angle,  and  its 
level  approaches  those  of  the  like-polarization  compo¬ 
nents  at  70°. 

In  a  separate  investigation,  the  backscattering  coeffi¬ 
cient  at  35  GHz  was  observed  as  a  function  of  time  over 
a  two-week  period  for  a  canopy  of  deciduous  trees  (bur 
oaks).  The  observation  period  covered  the  autumn  senes¬ 
cence  stage  during  which  the  moisture  content  of  the  trees 
decreased.  The  temporal  response  of  the  backscattering 
coefficient  (Fig.  8)  exhibited  a  3-dB  change  in  level  be 
tween  October  2  and  4  as  the  leaves  underwent  a  rapid 
change  in  moisture  content. 

C.  Backscatter  from  Snow 

In  February  and  March  of  1989,  the  University  of 
Michigan  millimeter-wave  system  was  used  to  measure 
the  backscatter  from  snow  at  a  site  near  Ann  Arbor,  MI. 
Figs.  9  and  10  illustrate  the  angular  variation  of  <j0  for 
wet  and  dry  snow.  Fig.  9  shows  the  response  of  dry  me¬ 
tamorphosed  snow  with  a  crystal  size  of  approximately 
2.2  mm,  and  Fig.  10  corresponds  to  fresh  wet  unmeta¬ 
morphosed  snow,  with  crystal  size  of  approximately  1.0 
mm  and  a  gravimetric  water  content  of  approximately 
1.9%. 

While  the  angular  dependence  of  the  two  plots  is  sim¬ 
ilar,  there  is  a  4-5-dB  level  shift  between  the  like-polar¬ 
ized  responses  of  the  dry  and  wet  snow  targets.  At  these 
frequencies,  snow  is  predominantly  a  volume-scattering 
medium,  and  the  presence  of  liquid  water  in  the  snow  me- 
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Fig.  8.  Two-week  temporal  response  of  35-GHz  backscatter  from  bur  oak  [Querces  macrocarpa)  trees  at  70°. 
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Fig  9.  Angular  response  of  bickscatter  from  dry  metamorphosed  snow  at  140  GHz. 


Fig.  10  Angular  response  of  backscatter  from  fresh  wet  unmetamorphosed  snow  at  140  GHz. 


dium  leads  to  increased  attenuation  and  decreased  albedo. 
For  the  dry  snow,  the  cross-polarized  response  is  lower 
than  the  like-polarized  response  by  about  4  dB,  while  for 
the  wet  snow  the  difference  in  level  is  about  6  dB. 

Fig.  1 1  shows  140-GHz  backscatter  measurements 
made  at  an  incidence  angle  of  40°  as  a  function  of  time 
over  a  12-h  interval  extending  from  noon  to  midnight  on 
February  27,  1989.  The  liquid  water  content  measured 
with  a  freezing  calorimeter  for  the  top  S-cm  snow  layer 
and  the  air  temperature  are  also  shown.  The  backscatter- 
ing  coefficient  is  observed  to  exhibit  a  3-dB  change  in  level 
at  around  14:30  in  response  to  the  decrease  in  tempera¬ 
ture  and  liquid  water  content.  At  140  GHz.  the  penetra¬ 
tion  depth  is  on  the  order  of  I  cm,  particularly  when  the 
snow  is  wet.  As  the  air  temperature  drops  below  0°C.  the 
snow  layer  starts  to  freeze  from  the  top  surface  down¬ 
ward.  Hence,  although  the  liquid  water  content  of  the  top 


5-cm  layer  may  still  be  greater  than  zero,  the  radar  re¬ 
sponds  only  to  the  top  1-2-cm  layer  and  therefore  exhibits 
a  time  response  that  leads  the  temporal  variation  exhibited 
by  the  measured  liquid  water  content.  This  dependence 
on  penetration  depth  is  illustrated  further  by  the  data  in 
Fig.  12  which  were  measured  by  the  33-GHz  channel. 
Because  of  the  greater  penetration  depth,  the  33-GHz  sys¬ 
tem  exhibits  a  much  more  gradual  change  in  level  be¬ 
tween  1400  and  2400  h.  Also,  the  magnitude  of  the  change 
in  level  is  12  dB  at  33  GHz,  compared  to  only  3  dB  at 
140  GHz.  At  94  GHz.  the  measured  diurnal  pattern  (not 
shown)  exhibited  a  response  similar  to  the  33-GHz  data, 
but  with  a  total  change  in  level  of  8  dB.  This  observed 
decrease  in  sensitivity  (of  the  backscattering  coefficient  to 
liquid  rate  content)  with  increasing  frequency  is  in  agree¬ 
ment  with  earlier  observations  reported  at  35  GHz  and 
lower  frequencies  [7]  .  [8]. 


498 


IEEE  TRANSACTIONS  ON  GEOSCIENCE  AND  REMOTE  SENSING.  VOL  28.  NO  4.  JULY  1900 


12  13  14  15  If  17  If  19  20  21  22  23  24 


Timo  (hours) 

(b) 


iwc  r%) 

Tomp.  (Ooq.  C) 


Fig-  11.  Diurnal  measured  temporal  variation:  (a)  140-GHz  radar  backscatter.  and  (b)  air  temperature  and  snow  gravimetric 

liquid  water  content  (top  5-cm  layer). 
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Fig  12  Diumal  measured  temporal  variation:  (a)  35-GHz  radar  backscatter.  and  (b)  air  temperature  and  snow  gravimetric 

liquid  water  content  (top  5-cm  layer). 
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V.  Conclusion 

This  paper  describes  the  operation  of  a  140-GHz  scat- 
terometer  system  with  a  measured  accuracy  of  1  dB.  Sam¬ 
ple  measurements  of  terrain  backscatter  at  35,  94,  and  140 
GHz  are  shown  for  grasses,  trees,  and  snow.  In  all  cases, 
the  angular  dependence  is  approximately  as  cos  0.  The 
two  like-polarized  components  (HH  and  W)  exhibit  es¬ 
sentially  identical  levels,  and  the  cross-polarized  response 
is  anywhere  from  1  to  6  dB  below  the  like-polarized  re¬ 
sponses,  depending  on  target  type.  At  140  GHz,  the  back 
scatter  from  snow  exhibits  a  dynamic  range  of  about  3  dB, 
in  response  to  change  in  liquid  water  content. 
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Millimeter-wave  (MMW)  remote  sensing  of  ground  snow  has  attracted  considerable  interest  in 
recent  years.  Because  the  size  of  the  snow  ice  particle  is  comparable  to  the  wavelength  in  the 
millimeter- wave  region,  we  can  no  longer  use  a  simple  Rayleigh  phase  function  or  the  small  particle 
approximation  usually  used  at  microwave  frequencies  for  calculating  the  extinction  coefficient.  In 
this  paper  we  present  a  model  for  MMW  scattering  from  snow  using  the  vector  radiative  transfer 
theory  and  a  Mie  phase  function.  Assuming  snow  to  consist  of  randomly  distributed  spherical 
particles  embedded  in  a  mixture  of  air  and  water,  the  vector  radiative  transfer  theory  is  solved  using 
the  discrete  ordinate  method.  The  values  of  the  extinction  coefficient  used  in  the  calculations  are 
based  on  a  combination  of  experimental  data  and  calculations  using  the  quasi-crystalline  approxi¬ 
mation.  The  backscattering  coefficient  is  calculated  for  different  liquid  water  contents  at  35,  95,  and 
140  GHz.  We  show  that  the  backscattering  coefficient  is  sensitive  to  liquid  water  content  at  all  three 
frequencies,  with  35  GHz  being  the  most  sensitive.  Except  for  normal  incidence,  the  effect  of  snow 
surface  roughness  is  negligibly  small  for  dry  snow,  and  it  is  somewhat  significant  for  wet  snow  at  35 
GHz.  but  not  at  the  higher  frequencies. 


1.  INTRODUCTION 

Although  several  papers  have  appeared  recently 
documenting  the  results  of  millimeter-wave  (MMW) 
radar  observations  of  snow-covered  terrain  [Stiles 
and  Ulaby ,  1980:  Ulaby  and  Stiles ,  1980;  Hayes  et 
al.,  1984;  Hallikainen,  1985 ,  Baars  and  Essen,  1988; 
Williams  et  al.,  1988;  Ulaby  et  al.,  1988;  Currie  et 
al.,  1988;  Narayanan  and  McIntosh,  1990],  the 
interaction  mechanisms  responsible  for  the  ob¬ 
served  radar  response  are  not  well  understood  at 
the  present  time.  At  millimeter  wavelengths,  snow 
is  a  highly  lossy  medium,  particularly  when  wet; 
consequently,  the  penetration  depth  is  only  of  the 
order  of  a  few  centimeters  [Ulaby  et  al.,  1986,  p. 
1608;  Hallikainen,  1985;  Hallikainen  et  al.,  1986]. 
For  dry  snow  the  attenuation  is  dominated  by 
scattering  because  the  ice  particles  are  comparable 
to  the  wavelength  in  size,  and  for  wet  snow  both 
absorption  and  scattering  are  important. 

The  physical  parameters  that  exhibit  the  stron¬ 
gest  importance  on  the  radar  backscatter  from  snow 
are  snow  surface  roughness,  crystal  size,  snow 
depth,  and  the  liquid-water  profile  with  depth.  This 
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paper,  which  is  part  1  of  a  two-paper  sequence, 
provides  a  radiative  transfer  model  for  characteriz¬ 
ing  MMW  scattering  from  snow,  using  experimental 
data  in  conjunction  with  the  quasi-crystalline  ap¬ 
proximation  [Tsang  et  al.,  1985,  pp.  461-469]  to 
compute  the  extinction  coefficient  of  the  snow  me¬ 
dium.  Part  2  [Ulaby  et  al.,  this  issue]  describes  the 
results  of  experiments  conducted  at  35,  95,  and  140 
GHz,  and  includes  comparison  between  theory  and 
experiment  for  certain  cases. 


2.  SNOW  MODEL 

In  our  snow  model  we  assume  ground  snow  to 
consist  of  spherical  ice  particles  embedded  in  a 
background  medium.  Liquid  water,  when  present, 
is  included  as  part  of  the  background  medium.  The 
size  of  the  water  inclusion  is  usually  much  smaller 
than  the  wavelength  for  millimeter-wave  remote 
sensing.  Therefore  it  is  reasonable  to  assume  that 
the  water  is  uniformly  distributed  in  the  snow  and 
the  dielectric  constant  of  water  can  be  included  as  a 
part  of  the  background.  Thus  the  wet  snow  medium 
is  modeled  in  terms  of  ir :  crystals  in  a  “wet  air” 
background.  It  is  also  possible  to  consider  the  snow 
as  lossy  particles  with  a  thin  film  of  water  sur¬ 
rounded  by  air,  in  which  case  the  dielectric  constant 
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TABLE  I.  Dielectric  and  Propagation  Properties  of  a  Snow  Medium  With  Ice  Volume  Fraction/  =  0.4,  Average  Ice  Particle 

Diameter  of  1  mm,  and  Number  Density  N  =  7.64  x  10® 


m(,  % 

^(jjj  *  ni 

K{  QCA) 

Re[/i(QCA)), 

m  1 

Kr(QCA),  m  1 

Real 

Imaginary 

Real 

Imaginary 

35  GHz  elce  =  3.15  +  /0.003  <r,  = 

=  1.2  x  10"® 

<T0  =  5  X  to-10 

*,<EFA)  = 

8.24  m'1 

0 

1 

0 

0 

930 

0.269 

1.268 

0.538 

I 

1.0427 

0.0129 

5.56 

942 

0.267 

1.285 

0.534 

2 

1.104 

0.0532 

22.2 

958 

0.264 

1.307 

0.528 

3 

1.1776 

0.0869 

35.2 

977 

0.262 

1.332 

0.524 

4 

1.2405 

0.1268 

50.0 

992 

0.259 

1.353 

0.518 

5 

1.3046 

0.1723 

66.2 

1007 

0.257 

1.373 

0.514 

95  GHz 

€1CC  =  3.15  +  *0.0085  <rt  = 

5.89  x  10-7 

<ra  =  4.28  x  10 

*f(EFA) 

=  402  m"! 

0 

1 

0 

0 

2562 

11.35 

1.2884 

22.69 

1 

1.0296 

0.0138 

16.2 

2582 

10.91 

1.298 

21.82 

2 

1.0626 

0.0397 

45.9 

2604 

10.43 

1.3095 

20.86 

3 

1.1055 

0.0697 

78.7 

2632 

9.84 

1.323 

19.68 

4 

1.141 

0.0961 

107.2 

2654 

9.39 

1.335 

18.78 

5 

1.1765 

0.1244 

136.6 

2676 

8.95 

1.346 

17.9 

140  GHz 

£jce  *  3.15  +  *0.012  a,  = 

1.95  x  10“6 

a*  =  2.33  x  10~8 

k,(EFA) 

=  1331  m" 

0 

1 

0 

0 

3876 

56.58 

1.32 

113.15 

1 

1.0254 

0.0115 

20 

3891 

54.95 

1.326 

109.9 

2 

1.0534 

0.0307 

52.67 

3909 

53.94 

1.331 

107.9 

3 

1.0883 

0.0540 

91.15 

3933 

53.36 

1.339 

106.7 

4 

1.1181 

0.0741 

123.38 

3955 

53.2 

1.347 

106.4 

5 

1.1482 

0.0953 

156.54 

3977 

53.1 

1.354 

106.3 

Here  m„:  liquid  water  content  by  volume,  c»:  background  dielectric  constant,  Kag:  background  absorption  coefficient.  K( QCA): 
propagation  constant  in  the  snow  layer  obtained  by  QCA,  Re(n(QCA)]:  real  part  of  the  effective  index  of  refraction  obtained  by  QCA, 
k,(QCA):  extinction  coefficient  of  ice  particles  obtained  by  QCA,  x,(EFA):  extinction  coefficient  of  ice  particles  obtained  by  EFA. 
a,  and  cra:  total  and  absorption  cross  sections  of  ice  particles  obtained  by  EFA. 


of  the  water  inclusions  is  part  of  the  dielectric 
constant  of  the  lossy  particles.  This  approach  is 
inappropriate,  however,  because  liquid  water  in 
snow  usually  occupies  the  spaces  between  adjacent 
ice  crystals  rather  than  coat  the  crystals  [Ulaby  et 
al. ,  1986,  p.  2067], 

From  ground-truth  data  we  know  that  the  ice 
particles  have  an  average  diameter  of  the  order  of 
0. 1-2  mm  and  their  shapes  are  round  but  nonspher- 
ical.  In  our  model  we  treat  the  snow  ice  particles  as 
spheres  with  a  normal-size  distribution  with  an 
average  diameter  of  I  mm  and  a  standard  deviation 
of  0.2  mm.  The  values  of  the  dielectric  constant  of 
the  ice  particles  needed  in  this  study  were  obtained 
from  Mdtzler  and  Wegmiiller  [1987]  and  are  listed  in 
Table  1.  Also  from  the  ground-truth  data  measured 
in  conjunction  with  the  experimental  observations 
reported  in  part  2,  the  volume  concentration  of 
snow  is  approximately  40%  and  the  snow  depth  is 
0.45  m.  These  values  are  used  in  our  model  calcu¬ 
lations. 


2.1.  Background  absorption  by  water 
inclusions 

The  imaginary  part  of  the  background  dielectric 
constant  is  directly  related  to  the  background  ab¬ 
sorption  coefficient.  We  use  a  dielectric  mixing 
formula  to  calculate  the  background  dielectric  con¬ 
stant,  and  we  assume  that  scattering  by  the  water 
particles  to  be  much  smaller  than  absorption.  The 
absorption  coefficient  of  the  background  me¬ 
dium  is  given  by 

=  2  Itn  {Ml  -/)}  (1) 

Kag  =  2*„(1  “/)  Im  {V^T}  (2) 

where  k0  is  the  free  space  wavenumber,  c*  is  the 
background  dielectric  constant,  and /is  the  volume 
fraction  of  ice  particles.  To  obtain  the  dielectric 
constant  of  the  background,  which  is  a  mixture  of 
air  and  water,  we  use  the  Van  Santen  mixing 
formula  [Ulaby  et  al .,  1986,  p.  1609]: 
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where  e„  is  the  dielectric  constant  of  water,  mv  is 
the  volumetric  snow  wetness,  and  Au  is  the  depo¬ 
larization  factor  which  depends  on  the  shape  of  the 
water  droplet. 

If  the  water  droplet  is  spherical.  Au  is  constant 
and  given  by  Au  =  1/3.  but  the  water  particles  in 
snow  are  usually  nonspherical  and  they  change 
shape  with  the  amount  of  liquid  water  present 
[Hallikainen  et  al .,  1986].  If  the  liquid  water  content 
is  low,  known  as  the  pendular  regime,  the  values  of 
Au  are  close  to  those  for  a  needle.  On  the  other 
hand,  if  the  liquid  water  level  is  high,  known  as  the 
funicular  regime,  the  values  of  Au  are  close  to  those 
for  a  disk.  Hallikainen  et  al.  [  1986]  provide  plots  for 
Aa,Ab,  and  Ac  as  a  function  of  mv ,  which  show  that 
the  transition  from  the  pendular  to  funicular  regime 
occurs  at  around  mv  =  2.5.  These  plots  are  used  in 
all  calculations  of  eb  in  this  and  the  following  papers 
(see  Table  l). 

2.2.  Scattering  characteristics  of  ice  particles 

The  scattering  characteristics  of  ice  particles  are 
calculated  using  the  Mie  solution  [ Ishimaru ,  1978, 
p.  27],  The  background  dielectric  constant  is  as¬ 
sumed  to  be  the  real  part  of  eb  listed  in  Table  1.  The 
imaginary  part  of  eb  is  not  used  because  of  the 
difficulty  in  calculating  the  Mie  solution  when  the 
imaginary  part  of  the  normalized  dielectric  constant 
is  negative.  The  average  total  and  absorption  cross 
sections  are  shown  in  Table  1. 

2.3.  Extinction  coefficient  of  dry  snow 

When  snow  is  dry,  the  attenuation  at  millimeter 
wavelengths  is  due  mainly  to  scattering  by  the  snow 
particles.  In  a  sparsely  distributed  medium  in  which 
the  correlation  between  particles  can  be  neglected, 
the  effective  field  approximation  (EFA)  can  be 
applied  and  the  extinction  coefficient  is  linearly 
proportional  to  the  concentration  of  particles 
]7sflrt(?  et  al..  1985.  p.  458].  Ground  snow,  however, 
has  an  ice  volume  fraction  of  10  to  40%,  and  the 
dielectric  constant  of  ice  is  much  larger  than  that  of 
the  background  medium.  Hence  the  independent- 
scatterers  assumption  is  inappropriate  for  snow. 


which  means  that  correlation  between  adjacent  ice 
particles  should  be  considered  [Tsang  et  al..  1985. 
chaps.  5  and  6].  Two  additional  important  phenom¬ 
ena  must  be  considered  when  treating  propagation 
in  a  dense  medium  like  snow,  namely  the  backscat- 
tering-enhancement  effect  and  the  decrease  of  the 
extinction  coefficient  when  the  density  is  high  [Ish¬ 
imaru  and  Kuga.  1982;  Kuga  and  Ishimaru.  1984: 
Tsang  and  Ishimaru.  1984], 

Backscattering  enhancement  is  caused  by  the 
constructive  interference  of  two  waves  propagating 
in  opposite  directions  and  is  important  only  in  the 
backscattering  direction  when  the  phase  difference 
of  the  two  waves  is  zero.  The  backscattering  en¬ 
hancement  effect  has  been  observed  for  both  ran¬ 
domly  distributed  discrete  particles  and  very  rough 
surfaces.  The  importance  of  the  backscattering  en¬ 
hancement  effect  that  has  been  recognized  in  optics 
and  solid  state  physics,  but  its  significance  in  mi¬ 
crowave  and  millimeter-wave  remote  sensing  has 
not  yet  been  evaluated.  Since  the  angular  width  of 
the  backscattering  enhancement  pattern  is  much 
less  than  1°  for  discrete  particles,  the  backscattering 
enhancement  effect  is  not  observable  if  the  detec¬ 
tor's  field  of  view  is  large.  In  most  microwave  and 
millimeter-wave  systems,  the  receiving  cone  of  the 
antenna  is  much  larger  than  1°  and  the  observation 
configuration  is  not  truly  monostatic.  Hence  the 
backscattering-enhancement  effect  is  probably  not 
significant  and  may  be  ignored,  at  least  to  first 
order,  although  we  suspect  that  it  may  be  respon¬ 
sible  for  the  difference  in  level  between  the  calcu¬ 
lated  and  observed  backscatter  at  95  GHz  and  140 
(antenna  beamwidth  about  2°)  reported  in  part  2  of 
this  sequence  [Ulaby  et  al.,  this  issue]. 

When  the  volume  fraction  of  particles  is  more 
than  1%,  the  extinction  coefficient  is  no  longer 
linearly  proportional  to  the  number  density.  The 
deviation  from  the  linear  relationship  applicable  at 
low  densities  is  related  to  the  size  parameter,  the 
dielectric  constant  and  volume  fraction  of  the  par¬ 
ticles.  Extensive  experimental  and  theoretical  stud¬ 
ies  on  the  extinction  coefficient  in  a  dense  medium 
have  been  conducted  in  recent  years  [ Ishimaru  and 
Kuga.  1982;  Tsang  et  al.,  1985].  The  theoretical 
models  include  Twersky's  model,  the  perturbation 
solution  with  hole-correction,  the  quasi-crystalline 
approximation  (QCA)  with  Percus-Yevick  pair  cor¬ 
relation  function,  and  the  quasi-crystalline  approx¬ 
imation  with  coherent  potential  (QCA-CP).  Twer- 
sky’s  model  is  simple  but  it  is  applicable  only  for 
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Water  Content  111^%) 

Fig.  1.  Optical  distance  r  versus  liquid  water  content  mv  at 
35  and  95  GHz.  Snow  depth  is  0.45  m,  the  ice  volume  fraction  is 
0.4,  and  average  ice  particle  diameter  is  I  mm.  EFA  is  the 
effective  field  approximation  and  QCA  is  the  quasi-crystalline 
approximation. 


small  particles.  The  formula  based  on  the  hole- 
correction  is  valid  if  the  volume  fraction  is  much 
less  than  10%.  For  higher  concentrations,  QCA  and 
QCA-CP  with  the  Percus-Yevick  pair  correlation 
function  have  been  shown  to  be  effective  [Tsang  et 
al .,  1985,  p.  479). 

In  the  millimeter-wave  region  where  the  size 
parameter  is  close  to  1  we  cannot  use  a  small- 
particle  approximation.  We  need  to  solve  the  QCA 
numerically.  We  calculated  the  extinction  coeffi¬ 
cient  of  snow  using  the  QCA  with  the  Percus- 
Yevick  pair  distribution  function  at  35,  95,  and  140 
GHz.  The  results  are  listed  in  Table  1.  When  the 
volume  fraction  /  is  0.4,  the  extinction  coefficients 
given  by  QCA  are  only  6.5%  of  those  calculated 
according  to  the  EFA  at  35  GHz  and  similar  per¬ 
centages  apply  at  95  and  140  GHz.  The  optical 
distance,  which  is  defined  as  r  =  (k,  +  Kag)d,  is 
shown  in  Figure  I  as  a  function  of  liquid  water 
content  for  both  EFA  and  QCA. 

According  to  recent  extinction  measurements 
conducted  for  dry  snow  at  35  and  95  GHz  [Halli- 
kainen .  1985).  the  extinction  coefficient  was  found 
to  exhibit  a  strong  dependence  on  snow  type.  The 
reported  values  of  the  extinction  coefficient  covered 
the  range  between  0.96  and  15.4  m_l  at  35  GHz  and 
between  1.9  and  30.7  m-1  at  95  GHz,  with  esti¬ 


mated  median  values  of  about  3.7  m~’  at  35  GHz 
and  19  m_l  at  95  GHz.  The  median  value  at  95  GHz 
is  close  to  that  computed  using  the  QCA  method, 
but  the  median  extinction  coefficient  at  35  GHz  is 
much  higher  than  that  computed  according  to  the 
QCA  method. 

If  the  optical  distance  is  greater  than  5.  the 
backscattering  coefficient  of  dry  snow  becomes 
essentially  independent  of  the  optical  distance  (in 
the  millimeter-wave  region).  For  a  snow  thickness 
of  0.45  m  and  extinction  coefficient  of  19  m-1  the 
optical  distance  is  already  more  than  8  at  95  GHz. 
Therefore  the  exact  value  of  extinction  coefficient  is 
not  important  at  95  and  140  GHz.  However,  an 
accurate  estimate  of  the  extinction  coefficient  is 
important  at  35  GHz  because  r  is  smaller.  In  our 
model  calculations  with  a  volume  fraction  /  =  0.4 
we  used  extinction  coefficients  of  3.7  m-1  at  35 
GHz  and  19  m  ” 1  at  95  GHz,  which  seem  to  offer  the 
best  fit  to  the  experimental  data.  For  snow  with 
other  values  of  /  the  extinction  coefficient  is  scaled 
in  proportion  to  the  change  in  the  QCA  model.  This 
method  is  chosen  because  the  extinction  coefficient 
is  a  nonlinear  function  of  f  and  the  density  depen¬ 
dence  curve  given  by  QCA  has  been  confirmed 
experimentally  [Ishimaru  and  Kuga ,  1982).  With 
this  correction,  for  example,  the  extinction  coeffi¬ 
cient  at  35  GHz  becomes  5.18,  4.45,  and  3.7  m~' 
when  the  snow  density  /  is  20.  30,  and  40%. 
respectively.  Because  of  the  lack  of  measured  ex¬ 
tinction-coefficient  data  at  140  GHz,  we  will  use  the 
value  calculated  by  the  QCA  model. 

3.  RADIATIVE  TRANSFER  THEORY 

Microwave  remote  sensing  of  ground  snow  has 
been  studied  by  many  researchers  in  the  past  [VI- 
aby  et  al.,  1986,  chap.  13;  Tsang  et  al.,  1985,  chap. 
4).  The  model  is  usually  based  on  the  radiative 
transfer  theory,  and  the  Rayleigh  phase  function  is 
used  for  modeling  the  ice  particles.  This  is  a  good 
approximation  for  microwave  remote  sensing  be¬ 
cause  the  ice-particle  size  is  much  smaller  than  the 
wavelength.  For  millimeter-wave  remote  sensing, 
where  the  ice  particle  size  is  comparable  to  the 
wavelength,  the  Mie  phase  function  must  be  used 
instead. 

In  a  dense  medium  like  snow  we  need  to  modify 
the  conventional  radiative  transfer  theory  to  take 
into  account  the  correlation  between  particles.  A 
dense-medium  radiative  transfer  theory  was  re- 
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cently  developed  using  the  Dyson  equation  with 
QCA-CP  and  the  Bethe-Salpeter  equation  under  the 
ladder  approximation  of  correlated  scatterers 
[Tsang  and  Ishimaru,  1987;  Wen  et  al.,  1990],  The 
form  of  the  dense-medium  radiative  transfer  theory 
is  the  same  as  the  conventional  radiative  transfer 
theory  and.  therefore  the  same  numerical  tech¬ 
niques  can  be  used  for  its  solution  [Cheung  and 
Ishimaru,  1982:  Ishimaru  et  al.,  1982;  Wen  et  al., 
1990],  The  difference  between  the  conventional 
radiative  transfer  theory  and  the  dense-medium 
radiative  transfer  theory  is  the  extinction  rate, 
which  can  be  obtained  by  the  QCA-CP  and  the  new 
form  of  the  albedo.  The  phase  function  is  still  the 
single  particle  phase  function,  which  can  be  approx¬ 
imated  by  the  Rayleigh  phase  function  when  the 
particles  are  small  in  size  compared  to  the  wave¬ 
length. 

In  this  paper  we  shall  use  a  modified  version  of 
the  conventional  radiative  transfer  model,  but  the 
formulation  is  not  as  rigid  as  the  one  given  by  Wen 
et  al.  [1990].  We  will  use  the  extinction  coefficient 
values  given  at  the  end  of  the  previous  section 
rather  than  those  based  on  the  effective  field  ap¬ 
proximation.  This  seems  to  give  the  best  fit  to  the 
experimental  data.  The  effective  dielectric  constant 
of  the  snow  layer  will  be  calculated  using  the  QCA 
and  the  real  part  will  be  used  for  calculating  the 
reflectivity  and  transmissivity  of  the  air-snow 
boundary.  The  phase  function  will  be  calculated 
using  the  Mie  theory.  The  background  absorption 
will  be  calculated  using  a  mixing  formula  and  in¬ 
cluded  in  the  total  extinction  coefficient.  This  will 
reduce  the  effective  value  of  the  albedo  when  the 
liquid  water  content  increases  in  the  background. 

3.1.  Problem  formulation 

We  consider  a  plane-parallel  medium  containing 
spherical  particles  as  shown  in  Figure  2.  A  linearly 
polarized  wave,  which  can  be  either  vertically  or 
horizontally  polarized,  is  obliquely  incident  with 
incident  angles  0o  and  <b0.  The  dielectric  constants 
of  media  1.  2.  and  3  may  be  different.  In  our  model 
medium  I  is  air.  medium  2  is  snow,  and  medium  3  is 
the  ground.  We  assume  that  e,  and  e3  do  not  change 
with  temperature,  but  ez  (the  background  dielectric 
constant  of  medium  2)  varies  in  response  to  changes 
in  liquid  water  content  of  the  snow  layer. 

The  formulation  of  the  radiative  transfer  theory 
for  an  oblique-incidence  case  has  been  derived 


previously  in  terms  of  the  Fourier-series  expansion 
for  the  Stokes  vectors  [Ishinaru  et  al.,  1982],  If  the 
incident  wave  is  normally  incident  and  linearly 
polarized,  only  two  terms  in  the  Fourier  series  are 
necessary,  but  the  oblique-incidence  case  requires 
all  the  components  of  the  Fourier  series.  In  this 
paper  we  will  briefly  describe  the  formulation.  The 
details  can  be  found  elsewhere  [Cheung  and  Ishi¬ 
maru,  1982;  Ishimaru  et  al.,  1982;  Ma  et  al..  1990]. 

For  convenience  we  will  use  the  modified  Stokes 
parameters  (/],  I2.  U,  V), 

/ 1  =  (£|Ef)  (4a) 

h  ~  (E^ET)  (4 b) 
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[5]  =  [Sy]  =  4x4  Mueller  matrix. 

[/,]  =  4  x  1  incident  specific  intensity  matrix 
(defined  later  by  (14)  and  (22)). 


In  a  sparse  medium  with  number  density  N  and  k , 
=  N{<rt)  the  above  equation  becomes 

x,  (<r,) 

—  = -  ( 10) 

k,  <tr,>  +  (KagIN) 


The  Mueller  matrix  (which  also  is  called  the  scat¬ 
tering  matrix)  is  expressed  in  terms  of  the  scattering 
amplitudes  /,|,/12,  and/22  [Cheung  and  Ishimaru , 
1982]: 

/(l/.il)2  <l/i:l)2 

m-_Ll  <l/2jl>2 

m_(<r,>  I  2  Re  </„/?,)  2  Re 

\2Im  (/„/?,)  2  1m  </l2/fe> 


This  form  is  similar  to  the  albedo  of  a  single  particle 
which  is  defined  as 


<T,  CTj  +  <7  a 


where  cr5  and  cra  are  the  particle's  scattering  and 
absorption  cross  sections.  If  we  define  the  effective 
albedo  as 


R  e</„/fc>  -lm(fnftz) 

Re</zi/?:>  -Im  iftxfW 

Re  {fwfh  +/«/?!>  -I mifuftz-fnfh) 

Im  (f\ |  f$2  +/t2/?i)  Re  (fufiz  -/i:/|i> 

where  <  )  denotes  ensemble  averaging  over  the 
particle  size  distribution  and  (cr,)  is  the  average 
extinction  cross  section  of  the  particles.  In  terms  of 
the  optical  distance  r  and  the  extinction  Ke ,  the 
equation  of  transfer  becomes 

d[/(r,  j)]  k,  r 

M  — 7 —  =  -[/(T,  i)}  +  -  ^'[S][/(T.  a*)] 

J4,r 


+  (7) 


with 


=  cos  0 

(8  a) 

da>'  =  dfi'  d<t>' 

(8b) 

dr  =  ds 

(8c) 

k’=K,  +  Kag 

(8 d) 

"5- 


_ 

(K„g/N) 


(12) 


we  observe  that  it  decreases  as  the  background 
absorption  increases.  For  example,  the  effective 
albedo  of  dry  snow  at  95  GHz  is  W'a  -  0.99.  If  the 
background  is  wet  and  the  snow  wetness  is  5%,  the 
effective  albedo  becomes  W’a  =  0.69.  A  similar 
calculation  at  35  GHz  shows  that  W’0  decreases 
from  0.95  for  dry  snow  to  W'0  =  0.09  for  a  snow 
wetness  of  5%.  In  our  model  we  assume  KelKe  to  be 
independent  of  the  density  and  given  by  (9).  For 
calculating  the  optical  distance  and  the  reflection 
coefficients,  however,  we  used  QCA  to  obtain  the 
effective  propagation  constant  of  the  snow.  We 
chose  this  approach  because  if  QCA  is  used  for 
calculating  k,/k^,  the  backscattering  coefficient  be¬ 
comes  very  small  for  slightly  wet  snow  and  it  does 
not  agree  with  the  experimental  data.  Our  approach 
seems  to  give  a  reasonable  fit  to  experimental  data, 
as  discussed  in  Part  2.  One  problem  with  our 
approach  is  that  when  r  becomes  small,  the  inten¬ 
sity  of  dry  snow  may  become  less  than  that  for  wet 
snow.  For  millimeter-wave  remote  sensing,  how¬ 
ever,  ris  large  and  this  deficiency  is  not  critical. 


Absorption  by  the  background  medium  (wet  air)  is 
now  included  in  k 't  in  (?).  If  the  background  absorp¬ 
tion  is  zero,  'D  reduces  to  (2)  in  the  work  by 
Ishimaru  et  al.  1982].  However,  if  Kag  is  not  zero, 
Kag  acts  like  an  additional  loss  term  in  the  particle's 
extinction  cross  section. 


3.2.  Incident  specific  intensity 
for  vertical  polarization 

Because  of  reflection  at  the  boundaries,  we  have 
upward  and  downward  traveling  incident  waves. 

[/,]  =  [/,*]  +  (/,']  (13) 


1l 


Kf  +  Ka„ 


where  [/,+  ]  is  for  p.  >  0  and  [/,  ]  is  for  p.  <  0.  For  a 
vertically  polarized  incident  wave  /„[  1  0  0  0]r  inci- 
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dent  along  the  direction  (mo>  0).  (7,+  ]  is  given  in 
terms  of  the  reduced  incident  intensity  I, .7  by 


[S][I-]dw- 


[/,*]=  [SIi;  ]  d» 


loCp  j 


[5]J  q  J8(m'  ~  n0)&{4>‘)e  M„  dn'  d<t>' 


=  I0C 


F 

p\  f31 

\f4, 


-  1oCpR2>  J  [5](  q 


8(n‘  +  n„)8{it>')  d/j.'  d<t>'  exp  (  - 


F2\ 

=  IoCpRv  F 


(\fn\-) 


Re  </ii/!i> 


<l/:il2> 


2  1m  </„/?,) 


3.3.  Incident  specific  intensity 
for  horizontal  polarization 

For  a  horizontally  polarized  incident  wave  given 
by  7o(0  1  0  0)r,  [7/+]  is  given  by 


H0  =  cos  0o,  4>a  =  0;  F||,  F2l.  Fji,  F4,  are 
evaluated  at  n  =  p.0\ 


[If  ]  -  IoC, 


( F 

r22 

pi  F32 


expl--  (23, 


1  -  « 23*21  exp  I - 

\  Mo 


n cos  6 

Tn=~ - f\'h\2 

n  |  cos  #i 


*21  =  kjil2 


<23  =  lr23l 


<l/.2|2> 


2  Re  (fnfh) 
(cr,) 


O/22 1 2> 


2  Im  </,:/*> 


The  reflected  specific  intensity  for  horizontal  polar¬ 
ization  is  given  by 


where  the  superscript  v  denotes  vertical  polariza¬ 
tion,  nj  is  the  index  of  refraction  of  medium  j,  and  rjj 
and  t}i  are  the  Fresnel  reflection  and  transmission  =  r  r  1 ? 

coefficients  at  boundary  ij. 


nf  cos  $;  -  n,  cos  $. 

rti  =  - - : - T  (20) 

r\j  cos  0,  +  fij  cos  6j 

2/1,  cos  6j 

fj  = -  (21) 

nj  cos  6,  +  n,  cos  0j 

The  reflected  incident  specific  intensity  [/,"]  is 
given  by 


exp  - 


2  T0  ~  T 


The  quantities  Cp  and  R2 3  in  (23)  and  (24)  are  given 
by  the  same  expressions  given  previously  for  the 
vertical  polarization  case,  except  for  replacing  the 
superscript  v  with  h  everywhere  and  using 

n,  cos  6.  -  nj  cos  6, 

r*  - - - - -  (26, 

n ,  cos  6,  +  rij  cos  0y 
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2 ni  cos  Oi 

t-.  - - 

n,  cos  d,  +  rtj  cos  6j 


3.4.  Outline  of  the  solution 


where 


t/L  = 


Wo  =  [5,]0a 


The  upper  and  lower  boundary  conditions  for  the  [f , , 

incoherent  specific  intensity  [/]  inside  medium  2  are  [F  =  p2, 
given  by 


[F -)o  = 


[I(r  =  0.  ~n)  =  [^:i(-m)]U<t  =  0’  ^)] 


[/(r  =  t0.  -/i)]  =  =  r0,  m)] 


(28)  For  m  >  0 


-1  <  m  <0 


1/FI  /*l  K  f 

FJT~=  -[']«  +  -  -  -f 

dr  «,  J_,  K; 


0  <  **  <  1  •  CP[F  +  ]me  -  CpR2i[F~]me  ~ 


where  (F2|(m)]  and  are  the  reflectivity 

matrices  given  by  the  generic  form 


Re  (rtf*) 
Im  (rtf  *) 


-Im  (rtf/, 
Re  (rtf/) 


where 


It  is  possible  to  eliminate  the  d>  dependence  from  (7) 
and  obtain  a  new  equation  in  terms  of  r  and  fi 
[Cheung  and  Ishimaru .  1982;  Ishimaru  etal .,  1982]. 
First  we  expand  [/],  [5],  and  [F]  in  Fourier  series  in 
<t>.  For  the  plane  wave  case  it  has  been  shown  that 
ail  Fourier  components  of  the  equation  of  transfer 
are  independent  of  each  other,  and  for  the  vertically 
and  horizontally  polarized  incident  waves,  the  first 
two  terms  of  the  stokes  vector  are  even  functions  of 
<i>  and  the  last  two  terms  are  odd  functions  of  <t>.  The 
incoherent  specific  intensity  is  given  by  the  sum  of 
all  Fourier  terms, 


[/]  =  2  [/]m  cos  m<t>  +  2  V)m  sin  m<t>  (31) 


We  can  obtain  the  equation  of  transfer  for  each 
Fourier  component.  For  the  Fourier  component  m 
-  0  and  a  vertically  polarized  incident  wave, 

d[/]0  k,  r\  k, 

M  -  -(/]„  +  -  [L]0[l}0  du  +  -f 

dr  Kf  J_,  Kf 

•  Cp[F*]0e~"’»°'  *  CMF-)0e-{2T"  -  (32) 


[F  +  ]m  = 


[F-]m  = 


[L]m  = 


[S  ,]" 
-[Sit 


The  equation  of  transfer  given  by  (32)  and  (34) 
together  with  the  boundary  conditions  given  by  (28) 
and  (29),  constitute  the  complete  mathematical  for¬ 
mulation  of  the  problem.  An  analytical  solution  of 
the  integrodififerential  equation  given  by  (32)  and 
(34)  is  not  available,  but  it  is  possible  to  solve  the 
equation  of  transfer  numerically  by  several  tech¬ 
niques.  In  this  study  we  use  the  discrete  ordinate 
method.  The  details  of  this  technique  are  given  by 
Ishimaru  (1978,  chap.  11].  Once  [/]„  and  [/]m  have 
been  found,  the  total  incoherent  specific  intensity 
emerging  from  the  layer  at  a  point  above  the  air- 
snow  boundary  can  be  computed  from 


[/]  =  [f2.] 


cos  m<t> 
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x 

2 


m  -  I 


sin  m<t> 


r  =  0 


(36) 


where  [72|]  is  the  transmissivity  matrix  of  the  upper 
boundary.  If  the  incidence  angle  is  larger  than  the 
critical  angle,  [Ty]  =  0;  otherwise,  it  is  given  by 

[T,j] 


rij  COS  $j 
n3  cos  0, 


0  0 

1^1 2  0 

0  Re  (/£/$*) 
0  Im  (#*•) 


:  \ 

-Im  (#$*)  I 

Re  <#£*>  / 


(37) 


where  t, j  is  given  by  (21)  for  vertical  polarization 
and  by  (27)  for  horizontal  polarization. 

Once  the  specific  intensity  has  been  computed  for 
a  vertically  polarized  incident  wave,  the  backseat- 
tering  coefficients  of  the  snow  volume  for  inj  and  hv 
polarizations  can  be  obtained  from  the  first  and 
second  elements  of  [/], 


=  4 ir  cos  0OI\ 
CTj,',’  =  Air  cos  Bgh 


(38) 


A  similar  solution  is  used  for  computing  cr** 

h 


and 


4.  SURFACE  SCATTERING  AT  BOUNDARIES 

In  the  millimeter-wave  region,  most  natural  sur¬ 
faces,  such  as  ground  snow,  are  rough  compared  to 
the  wavelength.  Surface  scattering  at  the  snow- 
ground  interface  is  not  important  for  millimeter- 
wave  remote  sensing  because  of  the  large  amount  of 
attenuation  in  the  snow  layer,  however,  the  scatter¬ 
ing  at  the  air-snow  interface  becomes  important 
when  the  snow  is  wet  and  the  volume  fraction  of  ice 
particles  is  high. 

If  the  top  surface  is  rough,  the  incident  wave 
transmitted  into  the  snow  layer  is  no  longer  a  plane 
wave.  Hence  we  need  to  consider  the  interaction 
between  surface  scattering  and  volume  scattering. 
This  interaction  has  been  examined  in  the  literature 
[Ulaby  et  al.,  1986,  chap  13;  Tsang  et  a/.,  1985,  p. 
203]  for  a  Rayleigh  layer  bounded  by  a  Kirchhoff 
rough-surface  interface.  On  the  basis  of  these  stud¬ 
ies  we  can  ignore  the  influence  of  the  rough-surface 


on  the  wave  transmitted  into  the  snow  layer  be¬ 
cause,  as  we  have  shown  in  Table  1,  the  effective 
index  of  refraction  of  the  snow  layer  is  around  1.3. 
However,  we  have  to  account  for  the  direct  back- 
scatter  from  the  rough  surface,  which  we  can  add 
incoherently  to  the  volume  backscattering  contribu¬ 
tion  derived  in  the  preceding  section. 

Two  widely  used  techniques  for  analyzing  sur¬ 
face  scattering  are  the  Kirchhoff  approximation  and 
the  small  perturbation  method  [ Tsang  et  al.,  1985. 
chap  2].  If  a  surface  has  a  large  rms  slope,  we  can 
obtain  a  simple  expression  using  the  Kirchhoff 
stationary-phase  approximation.  Assuming  that  the 
surface  characteristics  can  be  modeled  by  a  Gaus¬ 
sian  correlation  function,  we  can  write  the  copolar¬ 
ized  backscattering  coefficient  of  the  snow  surface 
as 


]r(0)l2  exp  [-(tan2  0o/2 m2)] 
2 m2  cos4  0O 


(39) 


where  d0  is  the  incidence  angle,  m  is  the  rms  slope, 
and  r(0)  is  the  Fresnel  reflection  coefficient  evalu¬ 
ated  at  normal  incidence.  The  Kirchhoff  approxima¬ 
tion  does  not  produce  depolarization  in  the  back- 
scatter  direction.  Therefore  the  total-snow 
backscattering  coefficient  is 

(Tu  =  (40a) 

=  <rj?  +  (406) 


—.hv  _  —.vh  _  vh  _  hv 
** rs  <Tfs  ® sv  &sv 


(40c) 


5.  RESULTS 

In  this  section  we  shall  discuss  the  variation  of 
the  backscattering  coefficient  with  snow  wetness 
(liquid  water  content),  snow  thickness,  and  angle  of 
incidence,  and  we  shall  evaluate  the  contribution  of 
surface  scattering  to  the  total  backscattering  coef¬ 
ficient.  All  calculations  were  performed  for  a  snow 
layer  with  an  ice  volume  fraction  of  0.4,  containing 
ice  particles  characterized  by  a  normal  size  distri¬ 
bution  with  a  mean  diameter  of  1  mm  and  a  standard 
deviation  of  0.2  mm.  The  snow  surface  is  assumed 
to  have  a  rms  slope  m  =  0.5. 

Evaluation  of  the  Fourier  components  of  (36)  led 
to  the  conclusion  that  it  is  necessary  to  include  only 
the  first  four  components  in  the  computation  of  the 
scattered  intensity  1  because  the  contributions  of 
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Liquid  Water  Content  mv  (%) 

Fig.  3.  Backscattcring  coefficients  <rvs  and  crts  versus  mv  at 
35  GHz.  Symbols  are  the  calculated  values  based  on  the  radia¬ 
tive  transfer  theory.  Solid  lines  are  third-order  polynomial  curve 
fits  to  calculated  results.  The  incidence  angle  is  60  =  40°  and  the 
incident  wave  is  vertically  polarized.  The  rms  slope  of  the  snow 
surface  is  m  =  0.5. 


the  higher-order  components  are  negligibly  small. 
Moreover,  except  when  the  snow  layer  is  dry  and 
its  thickness  is  only  a  few  centimeters,  the  effects  of 
reflections  by  the  underlying  ground  surface  may  be 
ignored. 

5.1.  Surface  scattering  contribution 

The  plots  shown  in  Figure  3  depict  the  35-GHz 
variation  of  the  snow-volume  backscattering  coef¬ 
ficient  crn,  with  liquid  water  content  for  a  0.45-m 
thick  snow  layer  at  an  incidence  angle  da  —  40°.  An 
additional  curve  is  shown  for  the  copolarized  case 
representing  the  total-snow  backscattering  coeffi¬ 
cient  crM,  which  includes  cr„,  and  the  surface  con¬ 
tribution  crsl.  We  observe  that  <7^  decreases  rapidly 
with  increasing  liquid  water  content  due  to  the 
corresponding  increase  in  background  absorption. 
The  rate  of  decrease  of  the  cross-polarized  compo¬ 
nent  is  much  higher  than  the  rate  for  the  copolarized 
component,  indicating  less  multiple  scattering  as 
the  medium  becomes  highly  absorptive.  Because 
we  use  spherically  shaped  ice  particles  in  the 
model,  single  scattering  does  not  produce  depolar¬ 
ization:  the  depolarized  return  is  caused  exclusively 
by  multiple  scattering. 

For  dry  snow  (m,.  =  0)  the  contribution  of  surface 


scattering  ass  is  negligibly  small  in  comparison  with 
asv.  However,  oss  increases  with  increasing  mv  and 
becomes  comparable  in  magnitude  with  osv  at  m,,  = 
5%.  These  observations  pertain  to  35  GHz.  Similar 
calculations  made  at  95  and  140  GHz  indicate  that 
ass  c  osv  even  at  mv  =  5%.  Thus,  surface  scatter¬ 
ing  need  not  be  included  at  these  higher  frequen¬ 
cies. 

The  variation  of  the  backscattering  coefficient 
with  liquid  water  content  at  95  and  140  GHz  is 
shown  in  Figure  4.  Compared  with  the  35-GHz 
results,  asv  decreases  with  liquid  water  content  at  a 
much  slower  rate,  and  the  cross-polarized  compo¬ 
nent  remains  significant  in  magnitude  even  at  high 
liquid  water  contents,  indicating  strong  volume 
scattering  within  the  snow  layer.  Although  the 
background  absorption  coefficient  Kag  is  approxi¬ 
mately  twice  as  large  at  95  GHz  than  its  value  at  35 
GHz  (see  Table  1),  the  effect  of  the  background 
absorption  on  the  total  backscattered  signal  is  less 
important  at  95  GHz  than  at  35  GHz  because  the 
scattering  cross  section  of  the  ice  particles  is  much 
larger  at  95  GHz,  as  a  result  of  which  multiple 
scattering  becomes  dominant  in  comparison  to  ab¬ 
sorption.  As  was  noted  earlier,  at  95  GHz  the 
effective  albedo  of  the  snow  volume  decreases  from 
0.99  for  dry  snow  to  0.69  for  wet  snow  with  mv  = 
5%,  compared  to  a  corresponding  decrease  from 
0.95  to  0.09  at  35  GHz.  At  140  GHz  the  effective 
albedo  decreases  from  0.99  for  dry  snow  to  0.90  for 
mv  =  5%. 

5.2.  Angular  variation 

The  angular  variation  shown  in  Figure  5  is  very 
close  to  a  cos  ^-dependence  for  both  the  copolar¬ 
ized  and  cross-polarized  scattering  coefficients.  The 
copolarized  response  represents  both  hh  and  vv 
polarizations  because  according  to  the  model  calcu¬ 
lations,  there  is  little  difference  between  the  two 
responses.  This  cos  0-like  variation  was  also  ob¬ 
served  at  95  and  140  GHz. 

5.3.  Variation  with  snow  thickness 

Figure  6  shows  the  response  of  the  backscattering 
coefficient  to  snow  thickness  at  35  and  95  GHz  for 
60  =  40°.  Two  wetness  conditions  are  shown.  In  all 
cases  the  minimum  thickness  considered  in  the 
model  calculations  was  5  cm.  For  dry  snow,  crv5 
increases  rapidly  with  increasing  snow  thickness, 
particularly  for  cross  polarization,  until  it  ap- 
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Liquid  Water  Content  mv  (%) 


Liquid  Water  Content  mv  (%) 


Fig.  4.  Backscattering  coefficient  arvs  versus  mv  at  95  and 
140  GHz.  Symbols  are  the  calculated  values  based  on  the 
radiative  transfer  theory.  Solid  lines  are  third-order  polynomial 
curve  fits  to  calculated  results.  The  incident  angle  is  $0  *  40°  and 
the  incident  wave  is  horizontally  polarized. 


proaches  a  saturation  level  beyond  which  the  in¬ 
crease  becomes  very  gradual. 

5.4.  Simulation  of  diurnal  response 

Let  us  consider  a  diurnal  cycle  during  which  the 
liquid  water  content  of  a  0.45-m  thick  snow  layer 
exhibits  the  variation  shown  in  Figure  la .  which  is 
Gaussian  shaped  with  a  peak  value  of  2%.  In  fact. 


Incident  Angle  (degrees) 

Fig.  5.  Backscattering  coefficient  <tvs  versus  incidence  angle  at 
35  GHz  for  a  0.45-m-thick  layer. 

the  wetness  need  not  be  uniform  in  depth  because 
the  radar  response  is  essentially  controlled  by  the 
wetness  of  the  very  surface  layer  (except  when  the 
surface  layer  is  dry  and  the  snow  contains  a  sub¬ 
surface  layer  with  nonzero  wetness). 


Snow  Thickness  (cm) 


Fig.  6.  Backscattering  coefficient  <rvs  versus  snow  thick¬ 
ness.  The  incident  angle  is  B0  *  40*  and  the  incident  wave  is 
vertically  polarized. 
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Fig.  7.  (a)  Model  for  the  liquid  water  content  variation  in  time.  The  coefficient  m,.  *  exp  ( — (/  - 

I3.85i:/cr:  where  =  2%  is  the  peak  value.  I  is  time  in  hours  and  <r  =  2  is  the  standard  deviation,  (h).  (c). 

and  id)  Calculated  diurnal  variations  of  the  backscattering  coefficient  a,.,  at  3$.  95.  and  140  GHz. 


Figures  lb ,  7c.  and  Id  show  the  corresponding  6.  CONCLUSIONS 

diurnal  radar  response  at  35.  95.  and  140  GHz.  All  Using  radiative  transfer  theory,  we  developed  a 
cases  include  both  the  copolarized  and  cross-polar-  reasonably  and  computationally  efficient,  uncompli- 
ized  responses.  As  expected,  the  radar  diurnal  cated,  model  for  relating  the  backscattering  coeffi- 
responses  are  approximately  mirror  images  of  the  cient  of  snow  at  millimeter  wavelengths  to  the 
liquid-water  diurnal  variation,  with  35  GHz  exhib-  physical  properties  of  the  snow  layer.  According  to 
iting  the  greatest  dynamic  range  and  140  GHz  the  model,  snow-surface  roughness  is  unimportant 
exhibiting  the  smallest.  when  the  snow  is  dry,  but  when  the  snow  is  wet 
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snow-surface  roughness  is  important  at  35  GHz  but 
not  at  the  higher  frequencies  examined  in  this  paper 
(95  and  140  GHz). 

Within  an  error  of  about  1  dB  the  backscattering 
coefficient  crls  varies  with  incidence  angle  as  cos  80 
for  both  the  copolarized  and  cross-polarized  config¬ 
urations  at  all  millimeter-wave  frequencies  and  liq¬ 
uid  water  contents  considered  in  this  study.  The 
effective  penetration  depth  of  the  snow  medium  is 
strongly  dependent  on  frequency  and  the  liquid 
water  content.  At  35  GHz  the  effective  penetration 
depth  decreases  from  about  30  cm  for  dry  snow 
down  to  about  1.5  cm  for  mv  =  5%.  At  95  and  140 
GHz  the  penetration  depth  of  dry  snow  is  5  and  I 
cm,  respectively.  The  backscattering  coefficient  de¬ 
creases  with  increasing  liquid  water  content,  with 
35  GHz  exhibiting  the  strongest  sensitivity  to  liquid 
water  variations  and  140  GHz  exhibiting  the  weak¬ 
est. 
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Using  a  truck-mounted  platform,  backscatter  measurements  were  made  at  35.  95,  and  140  GHz  for 
a  variety  of  snow  conditions  to  evaluate  the  radar  response  to  incidence  angle,  surface  roughness, 
and  liquid  water  content.  Good  agreement  was  obtained  between  the  experimental  observations  and 
theoretical  calculations  based  on  the  numerical  solution  of  the  radiative  transfer  equation  presented 
in  the  preceding  paper.  A  notable  exception  is  when  the  snowpack  is  in  the  refreezing  phase  of  the 
diurnal  cycle,  during  which  the  snowpack  is  characterized  by  a  dry  surface  boundary  with  wet  layers 
underneath.  To  accommodate  this  type  of  condition,  a  hybrid  first-order  numerical  solution  is 
proposed.  The  hybrid  approach  provides  excellent  agreement  between  theory  and  experiment. 


I.  INTRODUCTION 


In  Part  I  [Kuga  et  al.,  this  issue]  of  this  two-part 
sequence  we  proposed  a  radiative  transfer  model 
for  characterizing  radar  scattering  from  snow  at 
millimeter  wavelengths.  The  purpose  of  the  present 
paper  is  to  (I)  provide  a  summary  of  the  observed 
behavior  of  the  radar  backscattering  coefficient  cr° 
at  35,  95,  and  140  GHz,  including  its  dependence  on 
incidence  angle,  surface  roughness,  and  liquid  wa¬ 
ter  content,  (2)  compare  the  measured  data  to 
theoretical  calculations,  where  possible,  and  (3) 
propose  a  hybrid  first-order  numerical  model  for 
explaining  the  diurnal  variation  of  cr° . 

The  measurements  reported  in  this  paper  were 
acquired  by  the  University  of  Michigan's  truck- 
mounted  millimeter- wave  scatterometer  [Ulaby  et 
al.,  1988o;  Haddock  and  Ulaby,  1990].  A  summary 
of  the  system's  specifications  is  given  below. 


frequencies 
IF  bandwidth 
transmit  power 


sweep  rate 

polarization 
incidence  angles 
platform  height 


35.  94.  140  GHz; 

0  to  2.0  GHz; 

35  GHz,  +3  dBm; 

94  GHz.  0  dBm; 

140  GHz,  -4  dBm; 

1  m-s/frequency,  51,  101,  201,  401; 
frequencies/sweep ; 

HH.  HV,  VV,  VH; 

0*  to  70°; 

3  m  minimum,  to  18  m  maximum; 
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noise  equivalent  a-0 


crosspol  isolation 


phase  stability 


near-field  distance 


beamwidth 


antenna  diameter 


signal  processing 
output  products 


35  GHz,  -22  dB; 

94  GHz,  -28  dB; 

140  GHz,  -21  dB; 

35  GHz,  23  dB; 

94  GHz,  20  dB; 

140  GHz,  10  dB; 

35  GHz,  ~!  deg/hour; 

94  GHz,  ~1  deg/min; 

140  GHz,  —10  to  50  deg/s; 

35  GHz,  2.7  m; 

94  GHz,  7.3  m; 

140  GHz,  2.7  m; 

35  GHz,  R.  4.2°  T.  4.2°; 

94  GHz,  R:  1.4°  T:  2.8°; 

140  GHz.  R:  2.2°  T:  11.8°; 

35  GHz,  R:  6  in  T:  6  in; 

94  GHz.  R:  6  in  T:  3  in; 

140  GHz,  R:  3  in  T :  0.36  in; 
HP  8510A/8511A  based; 
received  power  versus  range; 
received  power  versus 
frequency  (at  fixed  RY, 
phase  and  amplitude  for  each 
frequency. 


Ground-truth  observations  were  made  for  the  fol¬ 
lowing  parameters:  (1)  air  temperature,  (2)  near- 
surface  snow  temperature,  as  well  as  the  tempera¬ 
ture  at  deeper  locations  in  the  snow  layer,  (3)  snow 
density  ps  (g/cm3),  (4)  height  profile  of  the  snow 
surface  using  a  graded  metal  plate  inserted  edgewise 
into  the  snow,  from  which  the  rms  height  s  is  calcu¬ 
lated,  (5)  depth  of  the  snow  layer.  (6)  volumetric  liquid 
water  content  of  the  surface  5-cm  layer. 
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measured  using  the  freezing  calorimeter  technique, 
and  (7)  microscope  photographs  of  thin  snow  sam¬ 
ples,  from  which  the  ice  crystal  size  is  estimated. 

The  standard  measurement  procedure  consisted 
of  measuring  the  backscattered  power  for  hh.  hv. 
vh.  and  w  polarizations  at  each  of  30  or  more 
spatial  locations.  The  backscattering  coefficient  was 
obtained  by  averaging  the  measurements  for  the 
different  spatial  locations.  In  addition  to  spatial 
averaging,  frequency  averaging  was  used  to  further 
improve  measurement  precision  [Ulaby  et  at.. 
19886].  The  estimated  uncertainty  associated  with 
the  values  reported  in  this  paper  is  ±0.5  dB.  Anal¬ 
ysis  of  the  data  shows  that  the  hv  and  vh  measure¬ 
ments  are  essentially  identical  (within  a  fraction  of  I 
dB),  which  is  expected  from  the  reciprocity  rela¬ 
tion.  In  almost  all  cases  the  copolarized  responses, 
< T^h  and  cr®  „  were  within  1-2  dB  of  each  other. 

2.  ANGULAR  RESPONSE 

The  data  shown  in  Figure  1  were  measured  for  a 
12-cm-thick  layer  of  dry,  freshly  fallen,  unmetamor¬ 
phosed  snow  composed  of  ice  crystals  with  diame¬ 
ters  on  the  order  of  1  mm.  The  measured  rms  height 
was  1.4  mm.  Only  hv  and  w-polarized  data  are 
shown  because  the  difference  between  and  <r°, 
is  I  dB  or  less  across  the  entire  angular  range  at  all 
three  frequencies.  The  curves  shown  in  the  figure 
were  calculated  according  to  the  theoretical  model 
described  in  the  preceding  paper.  For  w  polariza¬ 
tion,  theory  and  experiment  are  in  good  agreement 
at  35  and  95  GHz.  At  140  GHz,  however,  the  level 
predicted  by  theory  is  lower  than  the  experimental 
observations  for  w  polarization  by  about  4  dB.  We 
attribute  the  difference  to  the  backscattering  en¬ 
hancement  effect,  which  the  model  does  not  take 
into  account. 

A  typical  example  of  the  angular  dependence  for 
wet  snow  is  shown  in  Figure  2.  The  theoretical 
curves  were  computed  assuming  a  mean  ice  crystal 
diameter  of  1  mm  and  a  rms  slope  of  0.07.  Reason¬ 
able  agreement  between  theory  and  experiment  is 
obtained  except  for  w  polarization  of  140  GHz;  we 
again  attribute  the  difference  to  the  backscattering 
enhancement  effect,  although  no  strong  evidence 
exists  to  support  this  contention. 

3  EFFECT  OF  SNOW  SURFACE  ROUGHNESS 

According  to  the  model  results  presented  in  the 
previous  paper,  snow  surface  roughness  should 


Fig.  1.  Measured  and  calculated  backscattering  coefficient 
for  a  dry  snowpack  with  the  following  parameters:  depth.  12  cm; 
snow  density,  0.2  g/cm};  mean  crystal  diameter.  1  mm:  and  rms 
surface  slope.  0.07. 


have  a  minor  effect  on  the  level  of  tr°  (except  at 
normal  incidence)  when  the  snow  is  dry.  For  wet 
snow,  however.  cr°  should  increase  by  as  much  as  5 
dB  at  35  GHz  if  the  surface  is  made  rough  relative  to 
the  wavelength,  but  the  increase  should  not  be 
significant  at  95  GHz  or  higher  frequencies.  There  is 
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Fig.  2.  Measured  and  calculated  backscattering  coefficient 
for  a  wet  snowpack  with  the  following  parameters:  depth.  27  cm: 
snow  density.  0.4  g/cm1:  mean  crystal  diameter.  I  mm:  and  rms 
surface  slope.  0.07. 

plenty  of  experimental  evidence  to  support  the 
model  expectations  with  regard  to  dry  snow,  both  at 
35  GHz  I  Stiles  and  Ulaby,  1980]  and  95  GHz 
[Williams  et  al. .  1988].  Additional  support  is  pro¬ 
vided  by  the  data  shown  in  Figure  3  which  shows 
plots  of  <7°  versus  the  incidence  angle  0O  for  dry 


snow  at  95  GHz.  The  data  includes  measurements 
for  a  25-cm-deep  snowpack  made  before  and  after 
artificially  roughening  the  surface.  The  measured 
rms  height  was  found  to  be  0.7  cm  for  the  undis¬ 
turbed  surface  and  2.0  cm  for  the  roughened  sur¬ 
face.  Considering  that  A  =  0.3  cm  at  95  GHz.  both 
surface  conditions  are  electromagnetically  rough, 
so  it  is  not  surprising  that  the  angular  curves  shown 
in  Figure  3  for  the  two  surface  conditions  are  within 
1  dB  of  each  other  for  both  in>  and  hv  polarizations. 

A  more  detailed  examination  of  the  effect  of 
surface  roughness  on  the  radar  response  from  snow 
was  conducted  by  observing  the  radar  backscatter 
as  a  function  of  time  for  three  sections  of  a  71-cm- 
deep  snowpack  with  different  surface  roughnesses. 
Twelve  sets  of  observations  were  made,  all  at  0O  = 
40°.  for  each  surface  roughness,  starting  at  0630  LT 
and  ending  at  2200  LT;  each  set  consisted  of  mea¬ 
surements  at  35  and  95  GHz  for  all  linear  polariza¬ 
tion  combinations.  The  measured  rms  heights  of  the 
three  surfaces  were  jj  =  0.49  cm,  =  0.89  cm,  and 
s3  =  1.98  cm,  and  the  liquid  water  content  of  the 
5-cm  surface  layer  exhibited  a  Gaussianlike  varia¬ 
tion  with  a  peak  value  of  4.8%  (Figure  4). 

Figure  4  shows  a  sample  of  the  measured  data, 
specifically  the  cross-polarized  diurnal  responses  of 
all  three  roughnesses  at  both  35  and  95  GHz. 
According  to  the  35  GHz  data,  increasing  the  rms 
height  from  0.49  to  0.88  cm  (which  corresponds  to 
increasing  si  A  from  0.57  (slightly  rough)  to  1.02 
(rough))  causes  a0  to  increase  by  1-3  dB,  but 
increasing  the  roughness  further  to  s3  =  1.98  cm 
does  not  seem  to  have  much  of  an  impact  on  o-0.  At 
95  GHz,  even  the  least  rough  surface  (with  s\l A  = 
1.5)  is  electromagnetically  very  rough.  Hence  a0 
exhibits  approximately  the  same  diurnal  pattern  for 
all  three  surface  roughnesses.  A  similar  behavior 
was  observed  for  w  polarization. 

4.  DIURNAL  RESPONSE 

The  plot  shown  in  Figure  5  depicts  the  temporal 
variation  of  the  volumetric  liquid  water  content  of 
the  top  5-cm  snow  layer  over  a  period  of  14  hours, 
starting  at  0800  LT.  The  snowpack  was  dry  until 
1 100  LT.  then  mv  increased  rapidly  to  a  peak  value 
of  1%  at  1400  LT,  and  then  returned  to  the  dry 
condition  by  1900  LT.  This  description  applies  to 
only  the  top  5-cm  layer,  treated  as  a  single  layer;  it 
does  not  provide  any  information  on  lower  layers  or 
on  the  depth  profile  of  m,  within  the  5-cm  layer. 
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Fig.  3.  Measured  backscattering  coefficient  of  a  dry  snowpack  at  95  GHz  under  two  different  surface 

roughness  conditions. 


The  significance  of  this  statement  will  become  ap¬ 
parent  later  when  we  compare  measured  values  of 
mv  with  those  predicted  by  theory. 

Ideally,  it  would  be  desirable  to  measure  the 
entire  depth  profile  of  mv  with  a  vertical  resolution 
on  the  order  of  millimeters.  In  practice,  however,  it 
is  very  difficult  to  make  accurate  measurements  of 
mv  for  samples  thinner  than  S  cm.  Furthermore, 
because  it  takes  about  30  min  to  process  each  snow 
sample  (using  the  freezing  calorimeter  technique) 
and  because  mv  varies  rapidly  with  time,  it  is 
impractical  to  sample  more  than  one  depth  layer, 
unless  additional  trained  manpower  is  used  (two 
people,  working  full  time,  were  required  to  generate 
the  data  in  Figure  5). 

The  radar  response  to  the  observed  variation  in 
liquid  water  content  is  shown  in  Figure  6  for  35, 95, 
and  140  GHz,  measured  at  an  incidence  angle  of 
40°.  We  observe  that: 

1.  The  shapes  of  the  diurnal  patterns  of  <7°  are 
qualitatively  mirror  images  of  the  mv  pattern  for  all 
polarizations  and  frequencies. 

2.  The  dynamic  range  of  the  variation  is  greater  at 
35  GHz  and  smallest  at  140  GHz. 

3.  The  dynamic  range  of  the  variation  is  slightly 
greater  for  cross  polarization  than  for  like  polariza¬ 
tion. 

All  three  observations  are  consistent  with  the 
predictions  of  the  theoretical  model  given  in  the 


preceding  paper.  However,  as  we  will  see  later, 
quantitative  agreement  between  theory  and  experi¬ 
ment  for  the  melting  part  of  the  cycle  is  quite 
different  from  the  agreement  for  the  freezing  part  of 
the  cycle.  To  examine  this  question,  we  computed 
<7°  as  a  function  of  mv  using  the  numerical  solution 
outlined  in  the  preceding  paper  for  the  snowpack 
parameters  measured  in  the  field  (depth,  density, 
crystal  size,  etc.).  It  was  assumed  that  the  value  of 
mv  measured  for  the  top  5-cm  layer  is  valid  for  the 
entire  depth  of  the  snow  layer.  The  results  are 
shown  in  Figure  7  for  hh  polarization  at  0O  =40°. 
(Also  shown  in  Figure  7  is  the  computed  response 
of  <7°  based  on  the  first-order  solution  of  the  radia¬ 
tive  transfer  equation,  which  will  be  discussed  later 
in  section  4).  Allowing  for  a  possible  level  shift 
between  data  and  theory,  we  used  the  measured  cr° 
data  shown  in  Figure  6  (relative  to  <r°  of  dry  snow) 
to  compute  mv  according  to  the  numerical  solution 
curves  given  in  Figure  7.  The  plots  shown  in  Figure 
8  depict  the  diurnal  variation  of  mv  as  measured  in 
the  field  and  as  predicted  by  the  <7°  data  and  the 
model  calculations  at  35,  95,  and  140  GHz.  We 
observe  good  overall  agreement  between  the  mea¬ 
sured  values  of  mv  and  those  predicted  by  the  35- 
and  95-GHz  data  for  the  period  between  0800  and 
1300  LT,  and  although  the  values  predicted  on  the 
basis  of  the  140-GHz  data  are  somewhat  lower  in 
level  than  those  predicted  by  the  lower-frequency 
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Fig.  4.  Observed  diurnal  variation  of  mv,  the  liquid  water 
content  of  the  top  5-cm  layer,  and  <r°  at  35  and  95  GHz  for  cross 
polarization.  The  snowpack  was  28  cm  deep,  had  a  mean  crystal 
diameter  of  0.6  mm.  and  a  density  of  0.4  g/cm3.  The  measured 
rms  heights  of  the  surfaces  were:  i|  =  0.49  cm  (smooth),  r;  = 
0.89  cm  (slightly  rough),  and  sj  =  1.98  cm  (very  rough). 


channels,  the  temporal  variation  is  similar  to  that 
exhibited  by  the  other  plots. 

A  totally  different  situation  is  observed  in  Figure 
8  for  the  period  from  1300  to  1900  LT,  during  which 
the  snowpack  undergoes  a  refreezing  process.  In 
fact,  when  the  measured  mv  is  at  its  peak  value  of 
1%  (at  approximately  1400  LT),  all  three  radar- 
predicted  values  of  m„  are  in  the  0.8-1. 5%  range. 
Clearly,  the  models,  as  used  thus  far.  are  inappro¬ 
priate  for  modeling  the  backscatter  from  wet  snow 
under  these  "refreezing*’  conditions.  This  lack  of 
correspondence  between  theory  and  measurement 
is  caused  by  the  fundamental  assumption  that  mv 


Fig.  5.  Measured  air  temperature,  snow  surface  tempera¬ 
ture.  and  liquid  water  content  of  the  top  5-cm  layer  on  March  1 . 


1990. 


has  a  uniform  depth  profile  between  the  snow  and 
ground  surfaces.  During  the  melting  phase  of  the 
diurnal  cycle  the  melting  process  starts  at  the  air- 
snow  interface  because  the  source  of  thermal  en¬ 
ergy  is  the  warmer  air  mass  above  the  snow  pack  or 
direct  solar  radiation  (under  cloud-free  conditions). 
For  the  experiment  under  discussion  the  air  temper¬ 
ature  rose  from  -8°C  at  0800  LT  to  a  high  of  +6°C 
at  1145  LT,  and  then  decreased  down  to  -3°C  at 
2000  LT.  As  the  snow  surface  layer  starts  to  melt, 
its  thermal  conductivity  increases,  thereby  allowing 
the  transfer  of  thermal  energy  down  to  lower  layers. 
During  the  melting  phase,  mv  is  highest  at  the 
surface  and  decreases  in  an  exponentiallike  manner 
with  depth.  At  millimeter  wavelengths  the  penetra¬ 
tion  depth  in  dry  snow  is  of  the  order  of  1  m  at  35 
GHz,  5  cm  at  95  GHz,  and  1  cm  at  140  GHz;  and  for 
mv  =  3%,  the  penetration  depth  is  only  3  cm  at  35 
GHz  and  less  than  1  cm  at  the  higher  frequencies. 
Hence,  when  the  top  snow  layer  is  wet,  it  alone 
governs  the  radar  backscatter,  and  the  lower  layers 
exercise  an  inconsequential  influence  on  <r°. 

A  different  situation  occurs  during  the  refreezing 
phase.  As  the  air  temperature  starts  to  decrease,  we 
have  a  reversal  in  energy  balance  with  the  snow¬ 
pack  becoming  the  source  of  excess  energy.  Freez¬ 
ing  starts  at  the  surface  and  proceeds  slowly  down¬ 
ward  to  lower  layers.  This  results  in  an  inverted 
profile  with  mv  =  0  at  the  surface,  increasing  slowly 
to  a  maximum  value  at  some  depth  below  the 
surface,  and  then  decreasing  again  as  the  depth  is 
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Fig.  6.  Backscaf'ring  variation  measured  at  35.  95.  and  140 
GHz  on  March  1.  1990.  The  incidence  angle  was  40°.  The 
snowpack  density  was  0.32  g/cm}.  depth  was  12  cm.  and  other 
parameters  are  given  in  Figure  5. 


s 


Volumetric  Liquid  Water  Content  my  (  % ) 


Fig.  7.  Calculated  backscattering  coefficient  for  a  deep 
snowpack  versus  liquid  water  content  for  hh  polarization  at  an 
incidence  angle  of  40°.  The  calculations  are  based  on  the  first- 
order  and  the  numerical  solutions  of  the  radiative  transfer 
equation  for  a  snowpack  with  the  following  properties:  mean 
crystal  diameter.  I  mm:  and  snow  density.  0.32  g/cm}. 


increased  further.  Figure  9  shows  a  family  of  pro¬ 
files  generated  on  the  basis  of  this  simple  logic, 
constrained  by  the  condition  that  the  value  of  m,. 
averaged  over  the  5-cm  layer  has  to  be  equal  to  the 
value  measured  by  the  freezing  calorimeter.  With 
the  uppermost  surface  layer  being  dry  in  the  invert¬ 
ed-profile  case,  the  radar  response  will  be  governed 
not  only  by  the  dry  surface  layer,  but  by  the  wet 
layer  immediately  underneath  it  as  well.  It  is  not 
surprising,  therefore,  that  the  model-predicted  val¬ 


ues  of  mv  are  quite  different  from  the  5-cn.  average 
measured  in  the  field.  To  solve  this  problem,  we 
propose  to  use  a  hybrid  first-order  numerical  model 
for  characterizing  the  radar  response  of  snow. 

4.1.  Hybrid  first-order — numerical  model 

Let  us  consider  the  advantages  and  limitations  of 
the  first-order  solution  and  the  numerical  solution  of 
the  radiative  transfer  equation  as  we  apply  them  to 
the  snow  problem.  The  firs’  -'rder  solution  takes 
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Fig.  8.  Diurnal  variation  of  liquid  water  content.  The  plot 
labeled  "measured"  refers  to  the  value  of  mv  of  the  top  5-cm 
layer  as  measured  by  the  freezing  calorimeter.  The  other  three 
plots  refer  to  m.  as  predicted  by  the  model  using  the  radar 
observations  at  the  noted  frequencies  (with  hh  polarization  and 
0O  =  40°). 

into  account  only  single  scattering  in  the  snow 
volume,  whereas  the  numerical  solution  takes  into 
account  all  orders  of  multiple  scattering.  Conse¬ 
quently,  the  first-order  solution  may  underestimate 
the  magnitude  of  o-0.  particularly  when  the  scatter¬ 
ing  albedo  is  large.  The  scattering  albedo  decreases 
with  increasing  liquid  water  content  because  of  the 
increased  absorption  by  the  background  material 
(“wet”  air)  and  because  the  ratio  of  the  index  of 
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Fig  9  Proposed  depth  profiles  for  m. .  The  pack  is  assumed  to 
be  totally  dry  both  prior  to  1 100  LT  and  after  2200  LT. 
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refraction  of  the  ice  crystal  to  that  of  the  back¬ 
ground  decreases  with  increasing  snow  wetness.  To 
compare  the  first-order  solution  with  the  numerical 
solution,  we  refer  to  the  plots  in  Figure  7  which 
show  crffh  versus  mv  for  an  incidence  angle  of  40c. 
For  dry  snow  (m,,  =  0)  the  first-order  solution 
underestimates  cr°  by  about  3  dB  at  35  GHz.  8  dB  at 
95  GHz,  and  11  dB  at  140  GHz.  The  differences 
between  the  two  solutions  decrease  with  increasing 
mv.  becoming  insignificant  at  35  GHz  for  m.  >  29c. 
but  not  so  at  the  higher  frequencies.  Another  major 
advantage  of  the  numerical  solution  is  that  it  pro¬ 
vides  a  result  for  the  cross-polarized  scattering 
coefficient  <r°v,  whereas  the  first-order  solution  does 
not  predict  cross  polarization  because  the  scattering 
particles  are  assumed  to  be  spherical  in  shape. 

Having  stated  the  advantages  of  the  numerical 
solution  over  the  first-order  solution,  let  us  now 
consider  the  converse.  The  numerical  solution  re¬ 
quires  a  great  deal  of  computation,  in  spite  of  the 
fact  that  the  snow  medium  is  assumed  to  have 
uniform  properties  throughout.  If  we  are  to  formu¬ 
late  the  solution  for  a  multilayer  structure  in  order 
to  accommodate  a  nonuniform  depth  profile  for  m... 
the  complexity  of  the  numerical  approach  would 
make  the  solution  computationally  impractical.  On 
the  other  hand,  the  first-order  solution  is  perfectly 
amenable  to  computing  the  backscatter  from  a 
medium  with  nonuniform  properties  in  the  depth 
dimension.  Hence,  we  propose  to  use  a  hybrid 
model  that  takes  advantage  of  the  more  accurate 
feature  of  the  numerical  solution  and  the  easier 
structure  of  the  first-order  solution.  The  procedure 
is  as  follows. 

For  a  snowpack  with  specified  density  and  crystal 
size  distribution,  we  compute  o-0  for  the  semi- 
infinitely  deep  case  (at  the  angle  of  incidence, 
frequency,  and  polarization  of  interest)  as  a  func¬ 
tion  of  mv  (assuming  a  uniform  profile  with  depth) 
using  the  numerical  solution  outlined  in  the  preced¬ 
ing  paper.  We  shall  denote  this  backscattering  co¬ 
efficient  cr°(60.  /,  pq),  where  p  and  q  denote  the 
polarization  configuration  of  the  receive  and  trans¬ 
mit  antennas,  respectively.  Figure  10  shows  the 
results  of  such  a  calculation  at  35.  95.  and  140  GHz 
for  0q  =  40°  and  pq  =  hh.  The  form  of  the  variation 
can  be  fitted  to  the  functional  form: 

10  log  [cr^im, )]  =  a0  +  aim,1 '  +  a*m,  ( I ) 
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tional  techniques  is  in  tj,,  (or  its  “equivalent"  in  the 
case  of  the  numerical  solution).  We  propose  to  use 
the  results  of  the  numerical  technique  to  obtain  a 
more  exact  estimate  of  17,.  and  then  use  it  in  (2)  to 
compute  cr0  for  any  mv(z)  profile.  We  do  this  by 
evaluating  (2)  for  a  semi-infinite  medium  with  uni¬ 
form  liquid  water  content  mt.: 

0  17, Am,.)  cos  0o 

=  Tp(00 .  mv)Tq(6 0.  mx.) -  (3) 

2k,(/t!1.) 

and  then  equating  the  result  to  (1).  This  process 
leads  to 


Liquid  Water  Content  (  % ) 

Fig.  10.  Numerical  fits  using  the  form  given  by  (l ). 

where  a0 ,  a, ,  and  a2  are  constants  (for  a  given  set 
of  radar  parameters  (0O,  /,  PQ))- 
For  a  snow  medium  with  wetness  profile  mv(z) 
observed  at  an  incidence  angle  0O,  the  first-order 
solution  of  the  radiative  transfer  equation  leads  to 
the  following  expression  for  the  pq- polarized  back- 
scattering  coefficient. 

«7?  =  Tp[0o.wtl,(O)]r(,[0o.m,.{O)] 

x  J  tj,.  exp  £-2  J  Kf  sec  0'  dz  dz'  (2) 

where  Tp  and  Tq  are  the  transmissivities  of  the 
air-snow  boundary  for  backscattered  polarization  p 
and  incident  polarization  q,  m„(0)  is  the  snow 
liquid-water-content  at  z  -  0  (air-snow  surface),  d  is 
the  depth  of  the  snow  layer.  6'  is  the  refraction 
angle  in  the  snow  medium,  Ke  is  the  extinction 
coefficient,  and  Tjr  is  the  volume  backscattering 
coefficient  in  m^/m3.  Both  K'e  and  i)v  are  functions 
of  z  if  mv(z)  is  nonuniform.  The  form  of  (2)  is  very 
convenient  for  computing  the  backscattering  coef¬ 
ficient  when  mv(z)  is  not  uniform  with  depth,  but  as 
was  mentioned  earlier,  the  first-order  solution  un¬ 
derestimates  the  level  of  cr°  for  like  polarization, 
and  it  does  not  predict  a  cross-polarized  component 
because  tj,.  =  0  if  p  *  q  for  spherical  particles.  The 
transmissivity  coefficient  TpTq  is  identical  for  both 
the  numerical  and  first-order  solution,  and  the  same 
approach  used  for  computing  the  extinction  coeffi¬ 
cient  K'e  can  be  used  in  both  cases. 

The  key  difference  between  the  two  computa- 


2 <T°(mt.)«;(mr)  sec  0O 
Tp(d0,  mt,)7,(0o.  mv) 


2*c^(m,.)  sec  0O 

Tp(60,  mJTJOo’  mv) 


antilog 


Qq  +aimv*+a2mv 


Thus  tj ,,(m„),  as  given  by  the  above  expression, 
represents  the  “effective”  backscattering  coeffi¬ 
cient  per  unit  volume  of  the  snow  medium  at  the 
specific  frequency  /  and  polarization  combination 
pq  used  in  computing  <r°. 

To  summarize,  the  proposed  procedure  consists 
of  the  following  steps: 

1.  For  a  given  set  of  snow  parameters  (density 
and  particle  size  distribution),  the  QCA  method 
(adjusted  with  respect  to  experimental  data  as  dis¬ 
cussed  in  section  2.3  of  the  preceding  paper)  is  used 
to  compute  K'e(mv),  and  the  numerical  solution  is 
used  to  compute  tr°(m,,)  versus  m,,  (for  a  medium 
with  uniform  wetness  profile)  at  the  specified  wave 
parameters  (0O,  /,  pq). 

2.  For  convenience,  <7°(m,,)  is  fitted  to  a  function 
of  the  form  given  by  (1). 

3.  Equation  (5)  is  used  to  compute  77,, (m,,). 

4.  Equation  (2)  is  then  used  to  compute  o-0  for  any 
specified  profile  mv(z)- 


4.2.  Results 

This  procedure  was  applied  to  the  wetness  pro¬ 
files  shown  in  Figure  9,  and  the  results  are  shown  in 
Figure  11  for  hh  polarization  at  35,  95.  and  140 
GHz.  The  plots  compare  the  measured  diurnal 
variation  of  cr°  (relative  to  its  value  for  dry-snow 
conditions)  with  the  variation  computed  using  the 
hybrid  solution  outlined  above.  Very  good  overall 
agreement  is  observed  at  all  three  frequencies,  and 
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similar  agreement  was  obtained  for  the  other  polar¬ 
izations  and  for  other  diurnal  experiments. 

5.  CONCLUSIONS 

The  radiative  transfer  model  with  the  quasi-crys¬ 
talline  approximation  presented  in  the  preceding 
paper,  and  the  hybrid  first-order  numerical  solution 
proposed  in  this  paper,  together  provide  an  excel¬ 
lent  tool  for  examining  the  radar  response  of  snow 
at  millimeter  wavelengths.  This  conclusion  is  sup¬ 
ported  by  comparisons  of  theoretical  predictions 
with  experimental  observations  made  at  35,  95,  and 
140  GHz. 
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At  millimeter  wave  frequencies  a  typical  leaf  is  a  significant  fraction  of  a  wavelength  in  thickness, 
and  its  nonuniform  dielectric  profile  now  affects  the  scattering.  To  provide  a  simple  and  efficient 
method  for  predicting  the  scattering,  two  types  of  physical  optics  approximations  are  examined.  The 
first  approximates  the  volume  polarization  current  by  the  current  which  would  exist  in  an  infinite 
dielectric  slab  with  the  same  profile,  while  the  second  (and  simpler)  one  employs  the  surface  current 
which,  on  the  infinite  slab,  produces  the  known  reflected  field.  It  is  shown  that  the  first  method  is 
superior,  and  provided  the  actual  dielectric  profile  is  used,  it  predicts  the  scattered  field  to  an 
accuracy  which  is  adequate  for  most  practical  purposes. 


1.  INTRODUCTION 

Leaves  are  a  key  feature  of  any  vegetation  can¬ 
opy,  and  in  order  to  model  the  scattering  from 
vegetation-covered  land,  it  is  necessary  to  develop 
an  efficient  and  effective  technique  for  predicting 
the  scattering  from  a  single  leaf.  At  microwave 
frequencies  where  a  typical  leaf  is  electrically  thin 
with  lateral  dimensions  at  least  comparable  to  the 
free  space  wavelength  A0,  several  methods  have 
been  proposed  le.g.,  Le  Vine  et  al.,  1985;  Willis  et 
at.,  1988]  all  based  on  the  physical  optics  approxi¬ 
mation  applied  to  a  uniform  dielectric  slab.  In 
particular,  if  the  leaf  thickness  is  no  more  than 
about  A0/50,  physical  optics  in  conjunction  with  a 
resistive  sheet  model  predicts  the  scattering  at  most 
angles  of  incidence  [Senior  et  al.,  1987]  and  can  also 
handle  curved  leaves  [Sarabandi  et  at.,  1988]. 

On  the  other  hand,  at  millimeter  wavelengths  the 
thickness  can  be  a  significant  fraction  of  a  wave¬ 
length,  and  it  is  also  necessary  to  take  into  account 
the  internal  structure  of  a  leaf.  At  least  two  different 
types  of  cell  can  be  distinguished,  and  their  differing 
water  content  affects  the  dielectric  constant,  lead¬ 
ing  to  a  nonuniform  dielectric  profile.  To  compute 
the  scattering  at  these  higher  frequencies,  two  dif¬ 
ferent  physical  optics  approximations  are  exam¬ 
ined.  The  first  of  these  employs  the  polarization 
current  which  would  exist  in  an  infinite  slab  con¬ 
sisting  of  one,  two,  or  more  layers  simulating  the 
dielectric  profile  of  the  leaf,  and  this  is  referred  to  as 
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the  volume  integral  physical  optics  (VIPO)  approx¬ 
imation.  When  there  are  many  layers,  a  convenient 
method  of  implementation  is  described  in  the  ap¬ 
pendix.  In  the  second  (and  simpler)  approach  the 
equivalent  surface  currents  are  approximated  by 
the  electric  and  magnetic  surface  currents  of  an 
infinite  slab.  New  expressions  for  the  physical  op¬ 
tics  currents,  which  are  more  convenient  than  the 
standard  ones  [Beckmann,  1965],  are  introduced, 
and  this  technique  is  referred  to  as  the  surface 
current  physical  optics  (SCPO)  approximation. 

For  an  electrically  thin  leaf  or  plate  the  two 
approximations  are  indistinguishable,  but  as  the 
thickness  (or  frequency)  increases,  the  predicted 
scattering  differs  in  most  directions,  and  by  com¬ 
parison  with  the  results  of  a  moment  method  solu¬ 
tion  of  the  volume  integral  equation  it  is  shown  that 
VIPO  is  superior.  In  addition,  for  a  two-layer  ma¬ 
terial  it  is  no  longer  adequate  to  treat  the  plate  as  a 
homogeneous  one  having  an  average  dielectric  con¬ 
stant.  Provided  the  actual  dielectric  profile  of  a  leaf 
is  simulated,  it  appears  that  VIPO  can  predict  the 
scattering  behavior  of  a  leaf  to  an  accuracy  that  is 
sufficient  for  most  practical  purposes  at  millimeter 
wavelengths. 

2.  STRUCTURE  OF  A  LEAF 

The  structure  of  a  typical  vegetation  leaf  is  shown 
in  Figure  1 .  The  type  and  number  density  of  cells 
may  vary  as  a  function  of  depth  into  the  leaf  which, 
in  turn,  results  in  a  nonuniform  dielectric  profile. 
The  effect  of  this  nonuniformity  becomes  observ¬ 
able  at  higher  frequencies  where  the  thickness  of 
the  leaf  is  comparable  to  the  wavelength. 
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Fig.  1.  The  structure  of  a  typical  vegetation  leaf. 


Leaves  contain  two  types  of  photosynthetic  cells: 
palisade  parenchyma,  consisting  of  column-shaped 
cells  in  which  most  photosynthesis  takes  place,  and 
spongy  parenchyma,  consisting  of  irregularly 
shaped  cells  with  large  spaces  between  them  [Curtis 
and  Barnes,  1985].  Because  a  large  part  of  the 
vegetation  material  is  water,  its  dielectric  constant 
is  strongly  influenced  by  the  dielectric  constant  of 
water  and  the  water  content.  For  most  leaves  the 
water  content  is  higher  in  the  upper  layer  (palisade 
region)  than  in  the  undersurface  (spongy  region). 
The  sensitivity  of  the  dielectric  constant  to  water 
content  is  much  greater  in  the  lower  part  of  the 
millimeter  wave  spectrum  than  in  the  upper,  but 
this  is  more  than  counterbalanced  by  the  thickness 
to  wavelength  ratio.  The  net  result  is  that  the 
sensitivity  to  dielectric  variations  is  greater  at  the 
higher  frequencies. 

To  examine  the  effect  of  the  nonuniform  dielec¬ 
tric  profile  on  the  scattering  properties  of  the  leaf  at 
millimeter  wavelengths,  we  computed  the  normal 
incidence  reflection  coefficient  T0  of  a  two-layer 
dielectric  slab  and  compared  it  with  the  reflection 
coefficient  of  a  uniform  dielectric  slab  whose  dielec¬ 
tric  constant  is  the  average.  The  computation  was 
performed  for  a  leaf  thickness  of  0.5  mm,  and  the 
water  content  ratio  of  the  two  layers  was  chosen  to 
be  4  to  1,  representing  a  marked  variation  between 
the  upper  and  lower  surfaces  of  the  leaf.  From  the 
data  in  Table  1  it  is  seen  that  when  the  two-layer 
slab  is  approximated  by  a  uniform  slab  the  error  in 


the  reflection  coefficient  increases  with  increasing 
frequency  and  is  as  large  as  4  dB  at  140  GHz. 

3.  PHYSICAL  OPTICS  APPROXIMATIONS 

At  microwave  frequencies  where  a  typical  leaf  is 
no  more  than  about  Aq/50  in  thickness  with  lateral 
dimensions  comparable  to  or  larger  than  the  wave¬ 
length,  the  scattering  properties  can  be  accurately 
predicted  using  the  physical  optics  approximation 
applied  to  a  resistive  sheet  model  of  a  leaf  [Sara- 
bandi  et  al.,  1988].  In  effect,  the  leaf  is  modeled  as 
an  infinitesimally  thin  layer,  but  as  the  frequency 
increases,  it  is  necessary  to  take  the  leaf  thickness 
into  account.  There  are  now  two  types  of  physical 
optics  approximation  that  can  be  employed.  One  is 
the  surface  current  (SCPO)  approach  in  which  an 
infinite  dielectric  slab  is  replaced  by  an  equivalent 
sheet  current  that  produces  a  plane  wave  identical 
to  the  reflected  wave  of  the  slab.  This  current  is 
then  used  as  an  approximation  to  the  equivalent 
surface  current  over  the  upper  surface  of  a  finite 
dielectric  plate.  Alternatively,  the  induced  (volume) 
polarization  current  in  the  plate  can  be  approxi¬ 
mated  by  the  current  in  the  infinite  dielectric  slab, 
and  we  shall  refer  to  this  as  the  volume  integral 
physical  optics  (VIPO)  method.  It  is  more  accurate 
than  the  SCPO  method,  although  the  latter  is  more 
convenient  to  use  for  evaluating  the  scattered  field. 

To  illustrate  the  two  procedures,  consider  a  di¬ 
electric  plate  consisting  of  a  homogeneous  dielec¬ 
tric  of  thickness  d\  and  relative  permittivity  f|  atop 
a  second  material  of  thickness  d2  -  d\  and  relative 
permittivity  e2.  The  plate  occupies  the  region  -all 
s  x  s  all,  -b/2  <  y  s  bl2,  and  -d2  <  c  ^  0,  as 
shown  in  Figure  2,  and  is  illuminated  by  an  E- 
polarized  plane  wave  whose  electric  vector  is 

£i  _  j^iAoUsinflo  ~  zcosflo)  ( 

where  k 0  is  the  propagation  constant  in  the  free 
space  medium  above  and  below  the  plate.  When  the 


TABLE  I.  Voltage  Reflection  Coefficient  for  a  Two-Layer  and  Average  Dielectric  Slab 
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The  corresponding  results  for  a  single  layer  of 
thickness  dl  and  relative  dielectric  constant  gj  can 
be  obtained  by  putting  d2  =  dx  and  k2z  =  ku, 
implying B2  =  Bx  and A2  =  Ax. 

Given  a  volume  distribution  of  electric  current  J 
in  free  space,  the  corresponding  Hertz  vector  is 


n(r)  = 


iZ  o 
4vk.Q 


gikft\r  -  r‘\ 

Tr-r'l 


dv 


(6) 


Fig.  2.  The  geometry  of  the  scattering  of  a  plane  wave  from  a 
two-Jayer  dielectric  slab. 

plate  is  treated  as  an  infinitely  extended  slab,  the 
electric  field  can  be  written  as 

Ey  =  (e  “*°'z  +  Ye^y***™9**  0  <  z 

Ey  =  +Aie^)e^eax  _dl  <  0 

(2) 

Ey  =  ( B2e  ~ik2’z  +  A1eikl‘z)eik°s  ia9°*  -d2<z<-dx 

Ey  =  B3e  -**o^e*osm fl<*  z  <  _ d2 

where 


k0z  =  k0  cos  60  kjz  =  kQ(ej  -  sin2  0O)V2 

for  j  =  1,  2.  If  Ri  and  R2  are  the  reflection 
coefficients  at  the  upper  and  lower  surfaces  where 

„  koz-ku  n  k0z  -  k2z 
koz  +  kiz  kfjz  +  k2z 


C±  =l±-j^{l+R2e2ikl’{d2-d>y} 
k]z 

•{1  -R2e2ik^2-d,)y\ 

application  of  the  boundary  condition  at  the  three 
interfaces  gives 


C  +  (l  +  /?i) 
C+  +C-R,e2* 


Ax  =  ^e2^d'Bx  (3) 


e'i{ki2  -  ki2)di  2 


B3  =ei{k2’-ko')d2(l-R2)B2 


where  Zq(=  1/Yq)  is  the  free  space  impedance,  and 
the  resulting  field  is 

E(r)  =  V  x  V  x  n(r) 


H(r)  =  -ik0Y0V  x  n(r) 

In  the  far  zone  of  the  current  distribution. 


E(r)  =  -k$r  x  t  x  II(r)  (8) 

In  the  dielectric  slab  the  volume  current  J  is  the 
polarization  current 

J  =  -ik0Y0(ej  -  l)Eyy  (9) 

where  Ey  has  the  value  appropriate  to  each  layer 
(j  =  1,  2),  and  when  this  is  inserted  into  (6)  and  the 
integration  carried  out  over  the  volume  occupied  by 
the  plate,  we  obtain  the  VIPO  approximation.  For 
scattering  in  the  direction  8S  indicated  in  Figure  2 
the  expression  for  the  Hertz  vector  is 


n vlpo  =  y 


e‘k°r  k0ab  sin  X 
X  1 


where 


cos  Bs) 


{j  _  e-i(kiz  -koco&6s)d\ 

i(ku  -  k0  cos  6 

J  _  +*ocos«,)rf|  "I 

_  _ _  JJ  l 

Uku  +  A0  cos  0S)  1 J 

’ 6  ~i(kiz  -  kocos0i)d 

+ (e*  ■  ‘a — 


i  _  e~i(ki2  -kocoses)d2 


k0  cos  es) 


C  +  R\  +C.e2lkud 
C~'C~R\en^ra 


e*(/f22  +^ocos^>/!  _  €i{kiz  +kQcases)d2 

i(k2z  +  k0  cos  6S) 
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kQa 

X  =  ~  (s,n  +  SIn  eo) 


The  far  zone  scattered  field  can  then  be  obtained 
from  (7)  and  written  as 


where  £,  =  ( ik^  ej)~^  Z^bHyld  x  and  Ez  = 

~(ik0 £j)  ~  ZQdHy/d  x  have  the  values  appropriate  to 
each  layer  (j  =  1,  2).  The  Hertz  vector  can  be 
computed  using  (6),  and  for  the  scattered  field  HJ 
the  far-field  amplitude  is  found  to  be 


Es  =  j^Se(0s,0o)  (13) 

where  0O)  is  the  far-field  amplitude,  and  for 

the  VIPO  approximation  the  result  is 

klab  sin  X 

Sf?°{es,6*)=1-~-irF  (14) 

In  terms  of  the  far-field  amplitude  the  bistatic  scat¬ 
tering  cross  section  is 


<r(es,  0O)  =  —  IS(0„  0o)l2  (15) 

7 T 

The  more  conventional  SCPO  approximation  can 
be  obtained  by  noting  that  the  electric  current  sheet 


J  =  -2Po  cos  0oTe 


i*o  sin  0o*  s 


vido  kfcb  sin  X 

S%lpo(0s,  0o)  =  S~^r  —  (cos  0SF\  -  sin  0SF{) 


where 

F,  \ku{e i  -  I) 

k^e  i  i{ku  &o  cos  0S) 


l^ei»u+ha»»,)d,  1  k  ,  ,) 

_| _ D’  4.  ZZ—Z 

i(k\z  +  k0  cos  0S)  k0e2 


-ioco3«,)di  _  e~i(k\z  -  k0cose,)dj 


i(k2z  -  ko  COS  0S) 
e i(h ;  +  *ocose,)di  _  e‘(ku  +  kocos6,)di 

i(k2z  +  k0  cos  0S) 


produces  a  plane  wave  identical  to  the  field  re¬ 
flected  from  the  dielectric  slab.  As  evident  from  the 
impulse  function  S(z )  in  (16),  the  current  is  located 
at  the  upper  surface  of  the  slab,  and  when  (16)  is 
inserted  into  (7)  we  find 

e‘k,,r  -i  sin  X 

jjSCPO  ~  ^  __  -  —  cos  0or ab  — —  (17) 

and  the  far-field  amplitude  is  then 

pz'urt  t&o  sin 

SlCVO(0s  >  $o)  =  $  cos  e°rab  —  (18) 


Fj  =  sin  0< 


I  __  £-i(*iz  -*ocos0jMi  J  _  ei{k\z  +*ocos0J)dl 

_ _ __  A  r  _  _  n' 

i(k\z  -  *o  cos  0S)  i(k\z  +  k0  cos  0S )  1 


__  j\  {  -KAu-*oa»i,)di  _e-i(kit-kocos0j)d2 
_j_  L-i _ _  _ : _ ^ ' 

£2  i{k2z-k0  cos  0S) 


^l(*2z  +  *oCOS0,)t/i  _  giijkh  +  *oCOS0j)<^2 


i(k2z  +  k0  cos  0S ) 


In  the  specular  (0,  =  —  0O)  and  backscattering  (6S  = 
0O)  directions  it  can  be  verified  that  (14)  and  (18)  are 
identical,  but  in  the  other  directions  the  two  approx¬ 
imations  differ. 

In  the  case  of  H  polarization  for  which 

|| f  =  j^ffco(jrsin0o  “  zcos0o) 

the  analysis  is  similar.  With  Hy  represented  as 
shown  in  (2)  the  various  coefficients  (now  indicated 
by  primes)  differ  from  those  for  E  polarization  in 
having  klz  replaced  by  k\z/si  and  k2z  replaced  by 
k2zle2  everywhere  except  in  the  exponents.  The 
induced  polarization  current  then  has  two  compo¬ 
nents  and  is  given  by 

J  =  -ik0Z0(ej  -  l)(Exk  +  Ezi)  (20) 


The  SCPO  approximation  can  also  be  obtained  by 
noting  that  a  magnetic  current  sheet  of  the  form 

J*  =  -2Zq  cos  0orVfc°smfl°'r£te)?  (24) 

generates  a  plane  wave  identical  to  the  reflected 
wave.  When  we  use  this  as  the  equivalent  surface 
current  on  the  dielectric  plate,  the  magnetic  far-field 
amplitude  becomes 

ppnfi  -  ikl  sin  X 

s £cpo(0„  do)  =  S  cos  e0r ’ab  —  (25) 

As  in  the  case  of  E  polarization,  the  two  approxi¬ 
mations  are  identical  in  the  specular  direction,  but 
(21)  and  (25)  differ  in  all  other  directions,  including 
backscattering  (0S  =  0O)  unless  60  =  0. 
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Fig.  3.  Amplitude  of  the  ratio  of  the  bistatic  far-held  ampli¬ 
tude  of  V1PO  to  SCPO  for  £  polarization  of  a  dielectric  plate  with 
d2  *  V4  and  £|  *  ej  *  3  +  iO.l  at  6b  =  30°. 


Fig.  4.  Phase  of  the  ratio  of  the  bistatic  far-held  amplitude  of 
VlPO  to  SCPO  for  £  polarization  of  a  dielectric  plate  with  d2  = 
Afl/4  and  ej  *  e2  *  3  +  (0. 1  at  6b  «  30°. 


4.  NUMERICAL  RESULTS 

To  illustrate  the  difference  between  the  VlPO  and 
SCPO  approximation,  we  consider  a  homogeneous 
(single  layer)  plate  of  thickness  d2  =  A0/4  with  e2  - 
£i  =  3  +  <0.1.  For  an  £-polarized  plane  wave 
incident  at  30°  the  amplitude  and  phase  of  SViP0/ 
Sscro  ^  gjven  jn  Figures  3  and  4,  and  these  show 
that  the  difference  increases  away  from  the  specular 
and  backscattering  directions.  At  a  fixed  scattering 
angle  the  difference  increases  with  the  electrical 
thickness  of  the  plate  up  to  the  first  resonance  and 
then  decreases.  To  test  their  accuracy,  the  two 
approximations  have  been  compared  with  the  re¬ 
sults  of  a  moment  method  solution  of  the  volume 
integral  equation.  The  particular  code  used  is  a 
two-dimensional  one  which  was  extended  to  three 
dimensions  by  assuming  that  the  induced  currents 
are  independent  of  the  y  coordinate.  Since  the 
dielectric  constant  of  most  vegetation  materials  is 
high,  it  is  necessary  to  have  the  cell  sizes  very 
small,  and  one  consequence  of  this  is  the  need  to 
compute  the  matrix  elements  extremely  accurately, 
especially  for  H  polarization  [Sarabandi,  1989].  For 
a  2A0  square  plate  formed  from  the  above  men¬ 
tioned  layer  and  illuminated  by  an  £-polarized  plane 
wave  at  normal  incidence,  the  two  approximations 
are  compared  with  the  moment  method  solution  in 
Figure  5,  and  the  superiority  of  VlPO  is  clear. 

In  the  case  of  a  thin  plate  the  two  approximations 
are  indistinguishable.  This  is  illustrated  in  Figure  6 
showing  the  VlPO  expression  (14)  and  the  moment 
method  solution  for  a  2Aq  square  plate  of  thickness 


d2  -  Ao/50  for  £  polarization.  The  plate  is  a  homo¬ 
geneous  one  having  e  =  13  +  <12  corresponding  to 
the  average  permittivity  at  35  GHz  in  Table  1 .  The 
SCPO  expression  (18)  yields  the  same  results,  as 
does  a  two-layer  model  having  the  permittivities 
listed  in  Table  1.  The  analogous  data  for  H  polar¬ 
ization  are  given  in  Figure  7,  and  over  a  wide  range 
of  scattering  angles  the  approximate  and  moment 
method  solutions  are  in  excellent  agreement. 

As  the  frequency  and  hence  the  electrical  thick¬ 
ness  of  the  plate  increase,  the  superiority  of  the 
VlPO  approximation  becomes  apparent  and,  in 
addition,  it  becomes  necessary  to  take  the  layering 


Fig.  5.  The  bistatic  cross  section  of  a  2Aq  *  2Aq  plate  for  £ 
polarization  with  d2  =  A<>/4  and  €\  *  c2  *  3  +  rO.l  at  normal 
incidence:  moment  method  solution  is  denoted  by  the  solid  line, 
VlPO  by  the  short-dashed  line,  and  SCPO  by  the  long-dashed 
line. 
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Fig.  6.  The  bisutic  cross-section  area  of  a  2A«  x  2Ao  plate 
for  £  polarization  with  </2  “  A<,/50  and  eavf  *  13  +  il2  at  normal 
incidence:  moment  method  solution  is  denoted  by  the  solid  line. 
VIPO  or  SCPO  is  denoted  by  the  dashed  line. 

of  the  plate  into  account.  In  Figures  8  and  9  the 
simulated  frequency  is  140  GHz,  but  to  keep  the 
moment  method  calculations  tractable,  the  plate 
has  been  reduced  in  size  to  1.4Ao  by  2Ao-  The 
curves  shown  are  for  a  two-layer  plate  having  d2  = 
2d\  -  0.5  mm  with  £|  =  5  +  i4  and  e2  =  2  +  il,  and 
for  a  single  layer  having  the  average  permittivity 
£aVg  =  3.5  +  /2.5  (see  Table  1).  Since  the  accuracy 
of  the  physical  optics  approximation  increases  with 
the  plate  size,  the  agreement  between  the  two-layer 
VIPO  approximation  and  the  moment  method  solu¬ 
tion  is  remarkably  good,  and  significantly  better 
than  if  a  single  layer  had  been  used. 


Fig.  7.  The  bisutic  cross  section  of  a  2Aq  x  2Ao  plate  for  H 
polarization  with  <f2  -  Aq/50  and  *av|  »  13  +  <12  at  normal 
incidence:  moment  method  solution  is  denoted  by  the  solid  tine. 
VIPO  or  SCPO  is  denoted  by  the  dashed  line. 


Fig.  8.  The  bistatic  cross  section  of  a  1 .4Ao  x  2Aq  plate  for  £ 
polarization  with  d2  “  2 dt  m  0.5  mm  and /»  140  GHz  at  normal 
incidence:  The  solid  tine  denotes  moment  method  solution  with 
*i  *  5  +  »4,  e2  «  2  +  <1 ;  the  short-dashed  tine  denotes  VIPO  with 
<i=5  +  i4,  e2  *  2  +  il;  and  the  long-dashed  line  denotes  VIPO 
with  e2  *  <|  ■  3.5  +  <2.5. 

5.  CONCLUSIONS 

A  typical  leaf  has  at  least  two  dielectric  layers 
whose  cells  have  differing  water  content,  and  this 
produces  a  nonuniform  dielectric  profile  which  can 
now  affect  the  scattering.  At  microwave  frequen¬ 
cies  where  the  leaf  is  no  more  than  (about)  V50  in 
thickness,  the  nonuniformity  is  not  important,  and 
as  shown  by  Senior  et  al.  [1987]  the  leaf  can  be 
modeled  as  a  resistive  sheet  using  an  average  value 
for  the  permittivity.  If  the  physical  optics  approxi- 


Fig.  9.  The  bisutic  cross-section  ares  of  a  I  .4Aq  x  2Aq  plate 
for  H  polarization  with  dj  ■  2dt  m  0.5  mm  and /  *  140  GHz  at 
normal  incidence:  The  solid  tine  denotes  moment  method  solu¬ 
tion  with  C|  *  5  +  »4.  «2  »  2  +  il;  the  short-dashed  tine  denotes 
VIPO  with  *i  «  5  +  i4,  *2  *  2  +  fl ;  and  the  long-dashed  line 
denotes  VIPO  with  r2  ■  <|  »  3.5  +  »2.5. 
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mation  is  then  applied,  the  resulting  scattering  is 
attributed  to  a  surface  current,  and  this  method  is 
equivalent  to  the  SCPO  approximation.  At  higher 
frequencies,  however,  the  thickness  and  structure 
of  a  leaf  are  more  significant.  At  100  GHz  and 
above,  a  leaf  is  a  considerable  fraction  of  a  wave¬ 
length  in  thickness,  and  in  spite  of  the  reduced 
sensitivity  to  water  content,  the  nonuniformity  af¬ 
fects  the  scattering. 

For  a  two-layer  model  of  a  leaf  the  SCPO  approx¬ 
imation  has  been  compared  with  the  volume  inte¬ 
gral  (VIPO)  approximation.  When  the  leaf  is  thin, 
the  two  approximations  are  identical  and  in  good 
agreement  with  data  obtained  from  a  moment 
method  solution  of  the  integral  equation,  but  as  the 
electrical  thickness  increases,  the  two  approxima¬ 
tions  diverge  in  all  directions  except  the  specular 
and  (for  E  polarization)  backscattering  ones.  Al¬ 
though  the  VIPO  approximation  is  more  compli¬ 
cated,  its  accuracy  is  greater,  and  the  agreement 
with  the  moment  method  data  is  better  using  a 
two-layer  model  than  when  a  single  layer  of  average 
permittivity  is  employed. 

For  most  practical  purposes  it  would  appear  that 
VIPO  in  conjunction  with  an  accurate  dielectric  pro¬ 
file  of  a  leaf  provides  an  adequate  approximation  to 
the  scattering  at  millimeter  wavelengths.  As  our 
knowledge  of  the  profile  increases,  it  may  be  desirable 
to  use  a  multilayer  model  which  could  even  simulate  a 
continuous,  nonuniform  profile,  and  a  convenient 
way  of  doing  this  is  described  in  the  appendix.  We 
also  note  that  at  frequencies  for  which  the  leaf  thick¬ 
ness  is  comparable  to  Am/2,  where  Am  is  the  (average) 
wavelength  in  the  leaf,  the  scattering  is  greatly  re¬ 
duced  at  some  angle  of  incidence,  and  because  the 
permittivity  is  complex,  there  is  actually  a  range  of 
angles  for  which  this  is  true.  Since  the  reduction  is 
accompanied  by  an  increase  in  the  field  transmitted 
through  the  leaf,  this  could  provide  a  means  for 
penetration  through  a  vegetation  canopy. 

APPENDIX 

Al.  Combined  sheets  model 

When  using  the  VIPO  approximation,  an  efficient 
way  to  take  into  account  the  effect  of  any  nonuni¬ 
formity  in  the  dielectric  profile  is  to  model  the  leaf 
as  a  stack  of  N  combined  current  sheets.  Each  sheet 
simulates  a  very  thin  dielectric  layer  whose  thick¬ 
ness  is  less  than  A/ 1 5  where  A  is  the  wavelength  in 


the  material.  A  combined  sheet  consists  of  coinci¬ 
dent  resistive  and  modified  conductive  sheets  that 
support  electric  and  magnetic  currents,  respec¬ 
tively,  with  the  conductive  sheet  accounting  for  the 
electric  currents  flowing  perpendicular  to  the  di¬ 
electric  layer.  The  mth  layer  sheets  are  character¬ 
ized  by  a  complex  resistivity  and  conductivity  Rm 
and  R*m,  respectively,  where 


M-U*  -  1) 


R 


*  — 

m  ~~ 


•Ypem 

*o Aw(£m- ,) 


(Al) 


Here  en  and  Am  are  the  relative  dielectric  constant 
and  thickness  of  the  mth  layer,  and  t=2£=,  Am  is 
the  total  thickness  of  the  dielectric  slab. 

The  boundary  conditions  at  the  mth  combined 
sheet  are  as  follows  [Senior  and  Volakis,  1987]: 


Ax{Ax[E++E-]}=-2RMJm  (A2) 

Jm  =  ftx[H+-H-]  (A3) 

where  Jm  is  the  total  electric  current  supported  by 
the  resistive  sheet,  and 

A  x  {A  x  [H *  +  H-]} 

“^A^£[E*  +  E"3=  -2RZJS,  (A4) 

J*  =-Ax[E+-E‘]  (A5) 

where  JJJ,  is  the  total  magnetic  current  supported  by 
the  conductive  sheet.  The  superscripts  plus  and 
minus  refer  to  the  upper  (+)  and  lower  (-)  sides  of 
the  sheet,  and  A  is  the  unit  vector  outward  normal  to 
the  upper  side. 

A2.  The  induced  currents  of  N  planar 
sheets 

Consider  a  stack  of  N  infinite  planar  combined 
sheets  all  parallel  to  the  xy  plane  of  a  Cartesian 
coordinate  system  (x,  y,  z )  as  depicted  in  Figure  Al . 
The  top  sheet  is  in  the  z  =  0  plane,  and  the  mth 
sheet  is  located  at  z  =  -dm,  where  dx  =  0.  The 
space  between  the  mth  and  (m  +  l)th  sheets  is 
referred  to  as  region  m,  and  we  note  that  region 
0  (z  >  0)  and  region  N  (z  <  - d„ )  are  semi-infinite 
free  space.  A  plane  wave  whose  plane  of  incidence 
is  parallel  to  the  xz  plane  Impinges  on  the  stack  of 
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Fig.  A 1 .  Layer  of  N  combined  sheets  simulating  infinite  dielec¬ 
tric  slab. 


sheets  from  above.  From  the  symmetry  of  the  prob¬ 
lem,  all  the  field  vectors  are  independent  of  y  (i.e., 
d/dy  =  0),  as  a  result  of  which  the  field  components  in 
each  region  can  be  separated  into  E-  and  //-polarized 
waves  which  are  the  dual  of  each  other. 

The  analysis  is  similar  to  that  employed  in  study¬ 
ing  reflection  and  transmission  by  a  layered  half¬ 
space  [Kong,  1985],  In  the  case  of  E polarization  the 
incident  field  is  given  by  (1),  and  the  field  compo¬ 
nents  in  region  m  can  be  expressed  as 


£  _  | g—ikozosOnz  ^i^ocosflozj^^osinflo* 

Hmx  =  T0  cos  e0[C‘me-ik,>cose°z  -  c>,'*»cos8»V*°sinfloJt 


(A6) 


Hmi  =  Y0  sin  0o[C>-''*»cosfi»z  +  crmeik»c0*e°z]eik°sin9'‘x 

The  coefficients  Clm  and  Crm  are  the  amplitudes  of 
the  waves  traveling  in  the  —  z  and  +z  directions, 
respectively,  in  region  m.  In  region  0,  Cq  =  1  and 
Cq  =  T£  (the  total  reflection  coefficient),  and  in 
region  N,  CrN  =  0  and  C'N  =  TE  (the  total  transmis¬ 
sion  coefficient).  Hence,  using  the  boundary  condi¬ 
tions  (A2MA5),  there  are  2 N  unknowns  and  2N 
equations  that  can  be  solved  simultaneously.  On 
substitution  of  (A6)  into  (A4)  the  left-hand  side 
vanishes,  showing  =  0.  As  expected,  the  con¬ 
ductive  sheet  is  not  excited  with  this  polarization 
since  there  is  no  current  in  the  z  direction  in  the 
dielectric  slab,  and  in  the  absence  of  a  magnetic 
source  the  tangential  component  of  the  electric  field 
must  be  continuous  as  given  by  (A5).  On  inserting 
the  expressions  (A6)  into  (A5)  and  (A2)  and  defining 
the  reflection  coefficient  in  region  m  as 

cr 

pE  —  m  o  ”2Wrocos0o</m  +  i 

m  cL 


the  following  relations  are  obtained: 

_  -!  +  (2 Y0  cos  0„  Rm  _ ,  )e2tt«cos*i>(‘<- -  >  - 
-  1  1  +  2Yq  COS  0o  Rm  + 

(A7) 

rt  _ _ 1  +  -  I _  ci 

'-m  [  +  e2i*«cos#o(d„.i  -d„)p£  <-m-  1 

The  induced  electric  current  in  the  mth  sheet  can  be 
found  from  (A3)  and  expressed  as  (excluding  the 
phase  factor  el*osinfl°JC) 

jf  =  f2Y0  cos  0oe'*oeosSo</" 

'  i  +  rg 

J  4.  e2/*oCOS«cWn+!  -dm)p£ 


Vr  V  1  +  rf-i  A 

^  I  j  e2iAoCOs9oWi+ 1  —  <//)p£  J  (A8) 

The  total  reflection  coefficient  in  region  0  (T E(&)  = 
T<f)  can  be  evaluated  from  the  recursive  relation 
(A7)  by  noting  that  r$  =  0  (the  region  N  is 
semi-infinite).  The  total  transmission  coefficient  can 
also  be  obtained  from  (A7)  as  follows: 


ri  N 

tew-t?-  n 

Lq 


m  =  I 


j  ^2/ioCosfloW™.  i 


(A9) 


Unlike  the  E-polarized  case  where  the  magnetic 
current  is  zero,  an  //-polarized  wave  excites  a 
magnetic  current  in  the  y  direction,  and  the  tangen¬ 
tial  electric  and  magnetic  fields  are  both  discontin¬ 
uous  across  the  combined  sheets.  For  H  polariza¬ 
tion  the  tangential  field  vectors  in  region  m  can  be 
obtained  by  applying  the  duality  relationships  to 
(A6).  In  this  case  the  amplitudes  of  the  waves 
traveling  in  — z  and  +z  directions  are  denoted  by  B‘m 
and  Brm,  respectively.  By  applying  the  boundary 
conditions  (A2MA5)  at  the  mth  sheet  and  denoting 
the  reflection  coefficient  in  region  m  by 

Br 

ptf  =  _ Hi  „“2t*oCOS0O</m+  I 

1  m  D  i 

**m 

after  some  algebraic  manipulation  we  obtain 
VH 

1  m  -  1 


(QmPm 


(A  10) 
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(gw-  i)+(i+6m)r"-,  _  ■ 

m  -(l+em)+(l-Qm)r"e2‘ioCOs9oW"'*,-t'")  m_  1 


where  the  parameters  (2m  and  are 


sin2  0o 

2R*Z0  cos  0O 


2Rm  sec  0q 


Since  —  0  the  recursive  relation  (All)  gives  the 
total  reflection  coefficient  in  region  0,  and  using 
(A  12)  the  total  transmission  coefficient  is  given  by 


Fig*  A2.  The  geometry  of  the  scattering  of  a  plane  wave  from 
a  finite  TV-layer  combined  sheet. 


_ (g„-i) +(i+g„)r"-1 


The  induced  electric  and  magnetic  currents  can  be 
expressed  in  terms  of  the  reflection  coefficients  as 

=  __££J*oCOS0odm-l 


A3,  Scattering  by  a  rectangular  stack 

Consider  a  portion  of  the  ^-layered  stack  of 
combined  sheets  in  the  form  of  a  rectangle  occupy¬ 
ing  the  region  -a/2  ^  x  £  a/2,  -b/2  <  y  <  b/2  as 
depicted  in  Figure  A2.  In  the  far  zone  the  approxi¬ 
mation 


I?  —  I  **  r  +  sin  6,x'  +  cos  6sdm 


|^(1  +  r"_  ,)  +  (!  +  r"e2tto“sfloW"*1  “d")) 

(gm~  1)  +  (1  +  QmWm-\ _ ‘ 

a  +  Qm)  +  a  -  Qm)rykM{d-' 


mn‘  (Qi-DHi+QiWf-i 

'  (l + Qi) + ( i _ Qt)T^e 


J "*  =  JZo  cos  0O 


(i  -  r"_  ,)  +  (i  —  r"e2iA:(iCOsM‘/"'+i  _rf-)) 


(gm  -  1)  +  a  +  g„)r"_ , 

(T+  em)  +  (i  -  gm)r^"t»cos8o(rf",‘-<'" 


leads  to 


e  “r  iZ0  fa/2  fW2  /  £  , 


jlioCOSfljt/m 


^  or  /To  p>/2/  " 

mn-^LLiiwv 


.  e/Aosin0,x'  ^ 

Using  the  physical  optics  approximation,  the  cur¬ 
rents  obtained  for  the  infinite  sheets  are  substituted 
into  (A  16)  to  find  the  scattered  fields.  For  E  and  H 
polarizations  the  far-field  amplitudes  are 


s£(0J,  0o)  =  f-£-klabZ0l  2  S^Y~ 

\  m  =  1 


|-T  _ v> Cl  1  V1  1 

(i+gt)+a-g/)r,"e 


(gf-i)+(i+g/)rfL, _ 

4_f  i  _/r),^r?^^cos®o^+ 1  ~dl)  ^ 


S/j(0,  #o)  -  9  T~  kQab 


where 


J  *(*)  =  j"*(x)  = 


2  (cos  0,j£  +  Y0J%‘)eikMd"  —  (A18) 

m  =  1 
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where,  as  before,  X  =  {k0a/2)  (sin  6S  +  sin  60).  In 
the  backscattering  (0S  -  0O)  and  specular  (0,  = 
-0O)  directions  the  summation  term  in  (A  17)  re¬ 
duces  to  a  telescopic  series  resulting  in 

N 

2  j£e*tcos'*J”  =  2 To  cos  <*o  C'0  =  2T0  cos  0or£(0o) 

m  m  1 

(A19) 

and  the  backscattering  cross  section  is  then 

<T£(0 o.  ®o)  =  47T  —j~T~  COS2  0oir£(0o)l2 


tion  can  also  be  obtained  from  the  far-field  ampli¬ 
tude  using 

A2 

atxt  =  —  Im  [S(0O,  Oo  +  ^)]  (A23) 

7 r 

from  which  we  obtain 

<rfxt  =  2  ab  cos  0O  Re  [1  -  7*£(0)] 

(A24) 

<7 Jf1*  =  2fli  cos  0O  Re  [1  —  7W(0O)] 
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sin2  (k0a  sin  0O)^ 

(jfcofl  sin  Sq)* 

Also,  for  H  polarization, 
i  (cos  6oJZ  + 

m  -  1 

N 

=  -2  cos  «o  2  (b;  . ,  - 

m  ■  1 


(A20) 


=  -2  cos  $oBq  =  -2  cos  0O  r«(0o)  (A21) 

which  leads  to 


«r«(0o»  0o)=4ir 


cos2  0oir//(0o)l2 


sin2  ( ka  sin  0O) 
(ka  sin  0o)2 


(A22) 


These  results  are  identical  with  the  ones  obtained 
by  the  SCPO  technique.  The  extinction  cross  sec¬ 
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SUMMARY 

The  target  and  system  phase-stability  during  the  time  to 
measure  the  scattering  matrix  is  a  major  problem  for  millimeter 
wave  polarimetric  radars.  This  is  particularly  true  for  network 
analyzer-baaed  systems.  To  circumvent  this  phase^etability 
problem,  we  have  developed  new  fully  polarimetric  radars  at  35 
and  94  GHz.  The  system  is  based  on  a  relatively  inexpensive 
network  analyzer  and  La  capable  of  operating  in  either  the 
coherent  or  the  incoherent  polarimetric  measurement  mode.  In 
the  coherent  mode,  the  scattering  matrix  can  be  measured  within 
2  ms.  In  the  incoherent  mode,  the  average  Mueller  matrix  is 
measured  directly  by  transmitting  four  different  polarizations  and 
measuring  the  Stokes  vector  of  the  backscattered  signal.  To 
compare  the  performance  of  the  true  measurement  modes,  the 
average  Mueller  matrix  and  the  statistics  of  the  phase  difference 
of  the  two  co-polarized  signals  were  measured  for  a  rhododendron 
tree  and  for  a  metallic  tree  The  average  Mueller  matrices 
obtained  from  the  coherent  and  incoherent  polarimetric 
measurement  modes  were  similar,  The  target  motion  during  the 
data  acquisition  period  did  not  change  the  average  Mueller 
matrix  in  the  incoherent  measurement  mode.  The  probability 
density  function  of  the  phase  difference  of  the  two  co-polarized 
signals  computed  from  the  average  Mueller  matrix  is  essentially 
the  same  as  the  one  measured  with  the  coherent  polarimetric 
measurement  mode. 

1  INTRODUCTION 

Increasing  interest  has  been  expressed  in  recent  years  for 
understanding  the  statistical  properties  of  data  obtained  with 
fully  polarimetric  radars  for  remote  sensing  applications  [Ulaby 
and  Elochi,  1990].  At  centimeter  wavelepgths,  polarimetric  data 
has  been  found  to  be  useful  for  land-use  classification  [Van  Zyl  et. 
al.,  1987]  and  for  measuring  the  biophysical  properties  of  forest 
canopies  [McDonald  et.  al.,  1990].  For  the  MMW  region,  however, 
it  is  still  not  clear  what  type  of  information  can  be  extracted  from 
polarimetric  radar,  over  and  above  the  magnitude  information 
provided  by  conventional  radar  systems.  Unlike  the  microwave 
region,  the  complexity  and  the  cost  of  building  a  fully 
polarimetric  radar  at  millimeter-wave  frequencies  is  still  very 
expensive,  and  progress  has  been  rather  slow,  which  is  due,  in 
part,  to  the  limited  availability  of  experimental  data. 

At  microwave  frequencies  the  traditional  approach  used  for 
measuring  the  polarimetric  radar  response  of  a  given  target  is 
based  on  the  direct  measurement  of  the  target’s  scattering  matrix, 
S.  For  distributed  targets,  such  as  terrain  surfaces,  multiple 
measurements  of  5  are  made,  corresponding  to  statistically 
independent  samples,  each  measurement  is  used  to  compute  its 
corresponding  Mueller  matrix  £,  and  then  an  ensemble  average  is 
performed  to  obtain  an  estimate  of  the  average  Mueller  matrix, 

<  C  >.  Whereas  the  scattering  matrix  measurement  technique  is 
appropriate  at  microwave  frequencies,  it  is  difficult  to  implement 
at  millimeter  wavelengths  because  it  requires  that  both  the 
system  and  target  phases  remain  stable  during  the  time  it  takes 
to  measure  S.  This  is  particularly  true  for  network 
analyzer-based  polarimetric  radars  [Ulaby  et.  d.,  1990]. 

To  circumvent  this  phase-stability  problem,  we  have  developed 
new  fully  polarimetric  radars  at  35  and  94  GHz.  The  system  is 
based  on  relatively  inexpensive  network  analyzer  and  is  capable  of 
operating  in  either  the  coherent  or  the  incoherent  polarimetric 
measurement  mode.  In  the  coherent  mode,  the  scattering  matrix 
can  be  measured  withiir  2  ms.  In  the  incoherent  mode  the  average 
Mueller  matrix  is  measured  directly  by  transmitting  four  different 
polarizations  and  recording  the  horizontally  polarized  and 
vertically  polarized1  components  of  the  backscattered  field.  This 
paper  includes  a  detailed  analysis  of  the  two  measurement  modes, 
and  provides  comparisons  of  data  measured  using  the  two  modes 
for  a  rhododendron  tree  and  an  artificially  made  metallic  tree. 


2  NWA  BASED  POLARIMETRIC 
RADARS 

The  fully  polarimetric  radar  configuration  based  on  the  vector 
network  analyzer  (NWA)  is  easy  to  construct  and  is  widely  used 
for  remote  sensing  investigations  [Ulaby  et.  al.,  1990].  These 
systems  usually  are  operated  in  the  swept  frequency  mode  over  a 
given  bandwidth.  The  minimum  sweep  time,  which  depends  on 
the  number  of  frequency  points  and  the  type  of  NWA,  is  typically 
between  100  to  400  ms.  The  decorrelation  time  of  the  MMW 
wave  scattered  from  trees,  on  the  other  hand,  can  be  shorter  than 
10  ms  [Narayanan  et.  al.,  1988],  Hence,  when  using  the  fully 
coherent  measurement  configuration,  it  is  necessary  that  all  four 
components  of  the  scattering  matrix  be  measured  within  a  few 
milliseconds  in  order  to  obtain  accurate  data.  If  the  V-  and  H- 
polarized  signals  are  transmitted  sequentially  in  the  swept 
frequency  mode,  it  will  take  at  least  0.5  to  1  second  to  get  a 
complete  scattering  matrix,  including  the  data  transfer  time 
between  the  NWA  and  the  computer.  Obviously,  the  NWA-based 
MMW  radar  used  in  the  swept  frequency  mode  is  not  suited  for 
coherent  polarimetric  measurements. 

There  are  two  ways  to  overcome  the  shortcoming  of  the 
traditional  swept-frequency  NWA  based  polarimetric  radar.  The 
first  approach  is  the  incoherent  polarimetric  measurement 
technique.  With  this  technique  the  swept  frequency  mode  can  still 
be  used  for  the  NWA  operation  but  the  radar  transmitter  must 
be  modified  to  transmit  four  independent  polarizations.  The  data 
processing  and  calibration  are  substantially  more  complicated 
than  those  associated  with  the  coherent  polarimetric  technique. 
The  second  approach  is  the  coherent  polarimetric  measurement 
technique  using  Coupled/Chop  mode  and  point  by  point  external 
triggering  of  the  NWA.  We  have  developed  both  coherent  and 
incoherent  polarimetric  radars  based  on  these  techniques  at  35 
and  94  GHz.  The  radar  front  end  and  data  acquisition  system  are 
the  same  for  both  systems.  The  only  difference  is  the  operating 
mode  of  the  NWA  and  the  data  processing.  It  is,  therefore,  . .. 
possible  to  obtain  polarimetric  data  of  the  same  targets  \ 
coherently  and  incoherently.  The  block  diagram  of  the  MMW 
radar  system  and  the  35  GHz  front  end  are  shown  in  Figs,  1  and 
2.  The  block  diagram  of  the  94  GHz  is  essentially  the  same  as 
that  of  the  35  GHz  system.  In  the  following  section  the  details  of 
the  coherent  and  incoherent  systems  will  be  discussed. 

2.1  Coherent  Polarimetric  Radar 

The  coherent  polarimetric  radar  has  many  advantages  over  the 
incoherent  polarimetric  radar.  For  example,  with  the  coherent 
polarimetric  radar  the  statistical  data  including  the  phase 
difference  between  the  two  copolarized  channels,  can  be  easily 
obtained.  Another  advantage  is  the  significantly  simpler  signal 
processing  and  calibration  processes  compared  to  those  of  the 
incoherent  polarimetric  radar.  As  discussed  in  the  previous 
section,  the  NWA-based  radar  operated  in  the  swept  frequency 
mode  is  not  suited  for  coherent  polarimetric  measurements.  In 
this  section,  we  will  describe  a  new  technique  which  utilizes  a 
relatively  inexpensive  NWA  (HP8753C)  that  allows  the 
acquisition  of  coherent  polarimetric  data  at  a  much  faster  rate. 
With  this  system  it  is  possible  to  measure  the  scattering  matrix 
within  2  ms  at  35  and  94  GHz. 

The  Hewlett-Packard  network  analyzer,  HF8753C,  has  two 
independent  receiving  channels  which  can  be  used  in  the 
Coupled/Chop  mode.  It  also  has  a  point  by  point  external 
triggering  capability  in  the  swept  frequency  mode.  These 
functions  are  ideally  suited  for  the  coherent  polarimetric  radar. 
For  example,  the  simultaneous  acquisition  of  V  and  H  channels 
can  be  done  by  operating  A  and  B  inputs  in  the  Coupled/Chop 
mode.  The  point  by  point  external  triggering  can  be  used  for 
transmitting  V  and  H  sequentially  and  synchronizing  a 
polarization  control  circuit  to  create  different  polarizations.  At 
present,  HP8753C  does  not  support  the  external  point  by  point 
triggering  in  the  CW  mode  but  a  near  CW  mode  can  be  created 
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in  the  swept  frequency  mode  by  choosing  the  output  frequency 
bandwidth  to  be  1Hz.  The  minimum  time  to  get  a  complete 
scattering  matrix  is  approximately  2  ms  in  the  present  system. 
The  polarization  of  the  transmitted  MMW  signal  is  controlled  by 
a  Faraday  rotator  whose  switching  time  is  less  than  5  p a.  Using 
the  maximum  number  of  points  provided  by  the  HP8753C,  it  is 
possible  to  obtain  800  scattering  matrices  within  3.2  a  without 
transferring  data  into  a  computer. 

The  separation  of  signal  from  unwanted  noise,  such  aa  antenna 
coupling,  is  accomplished  by  the  hardware  gating  circuit  in  the 
IF  path  as  shown  in  Fig.  2.  The  transmitted  pulse  length  is  20  ns 
and  the  pulse-repetition-rate  is  5  MHz.  Although  it  is  not 
necessary  to  scan  the  RF  frequency  band  in  the  coherent 
polarimetric  mode,  additional  independent  samples  can  be 
realized  by  averaging  the  backscattering  coefficient  over  the  RF 
bandwidth  (Ulaby  et,  al.,  1988].  A  bandwidth  of  1  GHz  at  35 
GHz,  for  example,  offers  5  to  10  independent  samples  per  spatial 
observation  for  the  tree  measurements. 

The  calibration  of  the  coherent  system  is  straightforward. 
Because  the  system  has  more  than  23  dB  of  isolation  between  the 
V  and  H  channels,  a  simple  calibration  technique  that  requires  a 
sphere  and  a  depolarizing  target  is  used  [Sarabandi  et.  al.,  1990]. 

2.2  Incoherent  Polarimetric  Radar 

In  the  incoherent  polarimetric  radar  technique,  the  Mueller 
matrix  of  the  target  is  measured  directly  by  transmitting  four 
independent  polarizations  and  receiving  the  Stokes  vector  of  the 
scattered  signal.  Because  the  correlation  between  the  V —  and 
/f-  polarized  signals  is  inherently  included  in  the  received  Stokes 
vector,  the  measurement  time  between  the  different  incident 
polarizations  can  be  much  longer  than  the  decorrelation  time  of 
the  target.  The  incoherent  polarimetric  technique  also  permits 
the  use  of  MMW  sources  that  do  not  have  good  phase-stability  in 
the  transmitter  section  [Mead,  1990].  A  desired  polarization  can 
be  created  by  placing  two  quarter-wave  plates  in  front  of  the 
transmitting  antenna  and  by  adjusting  the  orientation  angle  of 
each  wave  plate  relative  to  the  incident  polarization. 

The  received  Stokes  vector  for  a  given  incident  polarization  is 
usually  obtained  by  two  different  approaches,  incoherent  and 
coherent-on-receive  techniques.  The  incoherent  receive  technique, 
which  often  is  employed  in  optics  measures  the  intensity  of  six 
different  receive  polarizations,  but  the  phase  measurement  is  not 
required.  The  Stokes  vector  is  obtained  by  taking  appropriate 
ratios  of  the  receive  intensities,  as  shown  in  Appendix  A. 

The  receiver  of  the  coherent-on-receive  technique  is  similar  to 
that  of  the  coherent  polarimetric  radar.  The  coherent-on-receive 
method  requires  the  measurement  of  the  magnitudes  of  the  V— 
and  H—  polarized  receive  signals  and  the  phase  difference 
between  them,  but  it  does  not  have  to  measure  the  phase  angle 
relative  to  the  transmitted  signal,  as  is  the  case  with  the  coherent 
polarimetric  radar.  The  Stokes  vector  can  be  computed  from  the 
magnitudes  of  the  V  and  H  components  of  the  received  signal 
and  the  phase  difference  between  them  as  shown  in  Appendix  A. 
Because  it  is  relatively  easy  to  measure  the  phase  difference 
between  the  V  and  H  channels,  our  system  is  based  on  the 
coherent-on-receive  technique. 

Calibration  of  incoherent  polarimetric  radar  systems  involves 
two  steps  [Mead,  1990].  In  the  first  step,  the  receiver  distortion 
matrix  is  obtained  by  placing  a  wire  grid  polarizer  in  front  of  the 
receiving  antenna  at  three  different  positions.  In  the  second  step, 
the  exact  polarization  properties  of  the  transmitter  are 
determined  by  measuring  the  backscatter  from  a  point  target 
with  known  scattering  matrix  using  the  calibrated  receiver. 


3  EXPERIMENTAL  DATA 


Table  1  shows  the  average  Mueller  matrices  of  the 
rhododendron  and  metallic  trees  obtained  by  the  coherent  and 
incoherent  polarimetric  radar  techniques.  The  Mueller  matrices 
are  normalized  with  respect  to  the  in  component  to  show  the 
relative  magnitude  of  the  matrix  elements.  The  average  Mueller 
matrix  of  the  coherent  polarimetric  radar  technique  was 
computed  from  the  8000  scattering  matrices  obtained  over  the 
34-35  GHz  band.  Because  of  the  slow  data-acquisition  speed  in 
the  incoherent  polarimetric  radar,  the  average  Mueller  matrix  is 
obtained  from  only  500  samples,  including  those  due  to  frequency 
averaging  over  the  1-GHz  RF  bandwidth. 

The  sum  of  the  Mueller  matrix  elements  L33  and  L44,  which  is 
a  function  of  the  correlation  between  Suv  and  Shh,  is  higher  for 
the  metallic  tree  than  for  the  rhododendron  tree.  Although  the 
polarimetric  signature  computed  from  the  average  Mueller  matrix 
is  useful  for  showing  the  characteristics  of  the  target,  it  is  not 
easy  to  directly  relate  the  target  characteristics  to  the  values  of 
the  Mueller  matrix  elements. 

Figure  4  shows  the  probability  density  function  of  the  phase 
difference  between  the  two  co-polarized  channels  (<6e  =  phase  of 
Svv  -  phase  of  Shh)  obtained  with  the  coherent  polarimetric 
radar.  As  expected,  p(^c)  of  the  metallic  tree  is  much  narrower 
than  that  of  the  rhododendron  tree,  showing  strong  correlation 
between  Svv  and  Shh- 

Unlike  the  coherent  polarimetric  radar,  the  information 
obtained  with  the  incoherent  polarimetric  radar  is  limited  to  the 
average  Mueller  matrix  and  it  is  not  possible  to  measure  the 
probability  density  function  p(<f»e)  directly.  Due  to  a  recent 
theoretical  derivation,  however,  the  phase  statistics  of  <f>e  can  be 
estimated  from  the  average  Mueller  matrix  [Sarabandi,  1991]. 

The  probability  density  function  p(<f>e)  is  given  by 
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The  function  p(^fi)  is  completely  specified  in  terms  of  the 
elements  of  the  average  Mueller  matrix  Cm. 

The  average  Mueller  matrix  given  in  Table  1  and  the 
probability  density  function  of  the  phase  difference  shown  in  Fig. 

4  are  obtained  from  the  same  target  by  two  different  polarimetric 
measurement  techniques.  If  the  probability  density  function  given 
by  Eq.  1  is  correct,  p(<f>c)  estimated  from  the  average  Mueller 
matrix  must  be  similar  to  the  one  Bhown  in  Fig.  4.  Figure  5 
shows  the  probability  density  function  computed  from  the  average 
Mueller  matrix  obtained  by  the  incoherent  polarimetric  radar. 

The  agreement  between  Figs.  4  and  5  is  excellent  for  both  trees. 


4  CONCLUSION 

The  work  described  in  this  paper  has  demonstrated  that  the 
average  Mueller  matrices  obtained  using  the  coherent  and 
incoherent  polarimetric  measurement  techniques  are  essentially 
identical.  The  advantage  of  the  coherent  polarimetric  radar  over 
the  incoherent  polarimetric  radar  is  its  ability  to  measure  the 
statistical  distributions  of  the  magnitudes  and  relative  phases  of 
the  scattering  matrix  elements.  The  incoherent  polarimetric 
radar,  however,  is  particularly  useful  if  the  target  decorrelation 
time  is  much  faster  than  the  data  acquisition  time. 


To  demonstrate  that  the  coherent  and  incoherent  polarimetric 
measurement  techniques  do  indeed  provide  identical  information 
for  distributed  targets,  experiments  were  conducted  using  a 
rhododendron  tree  and  a  metallic  structure  resembling  a  short 
tree.  Photographs  of  these  targets  are  shown  in  Fig.  3.  The 
metallic  structure  is  used  for  creating  a  target  return  with  strong 
correlation  between  the  Svv  and  Shh  components.  To  create  many 
independent  samples  and  also  to  show  that  the  incoherent 
polarimetric  technique  can  provide  accurate  results  even  if  the 
data  acquisition  time  is  much  longer  than  the  target  decorrelation 
time,  the  trees  were  rotated  at  slow  (0.67  rpm)  and  fast  (1.33 
rpm)  speeds  during  the  data* collection  process. 
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APPENDIX  A 

COHERENT-ON-RECEIVE 

TECHNIQUE 

Complete  polarimetric  characterization  of  the  scattering 
properties  of  a  distributed  target  can  be  obtained  by  measuring 
either  the  scattering  matrix  S  or  the  Mueller  matrix  Cm . 
Measurement  of  the  scattering  matrix  requires  accurate  phase 
measurements.  Also  4  elements  of  S  must  be  obtained  within  the 
decorrelation  time  of  the  target  which  is  in  the  order  of 
milliseconds  at  MMW  frequencies.  The  scattered  electric  field  Er, 
in  terms  of  the  scattering  matrix  S  and  the  incident  electric  field 
E‘ ,  is  given  by 


(XI) 

(A.2) 


To  obtain  5,  we  need  to  send  [£„,()]*  and  [0,  2?a)*,  and  measure 
Ev  and  Eh  simultaneously. 

The  polarized  wave  can  also  be  expressed  in  terms  of  the 
Stoke’s  vector  Fm  which  is  defined  as 


where  Cm  is  called  the  Mueller  matrix. 

The  totally  incoherent  method  does  not  require  phase 
measurements.  With  this  method,  the  4  elements  of  Stokes  vector 
are  obtained  by  receiving  6  polarizations  (V,  H ,  45, 135,  LHC, 
RHC).  For  example,  the  third  element  of  Stokes  vector  U  is  given 
as  a  ratio  of  intensities  at  45  linear  to  135  linear.  For  a  given 
incident  polarization,  therefore,  we  can  obtain  a  column  of  the 
Mueller  matrix.  To  get  the  complete  Mueller  matrix,  we  need  to 
repeat  this  process  for  4  independent  incident  polarizations. 
Altogether,  at  least  24  magnitude  only  measurements  are  required 
to  obtain  the  complete  Mueller  matrix.  Although  the  phase 
measurement  is  not  required  with  the  incoherent  method,  it  is 
necessary  to  receive  all  6  polarizations.  The  elements  of  the 
Stokes  vector,  in  terms  of  6  polarizations  and  a  set  of  4 
independent  incident  polarizations,  given  by: 


Wv 

~  Wh  +  Wv 

wh 

wh  +  wv 

w4i  -  Wise 
+  W13B 

_  Wlhc  -  Wrhc 
Wlhc  +  WRffc 


(A.6) 
(XT) 
(X8) 
(X  9) 


r 1 1 
0 

0 

,h  = 

r°* 

i 

0 

,  hf>  = 

ri/2] 

1/2 

1 

.  Ilhc  = 

ri/2i 

1/2 

0 

.0. 

.0, 

.  0  . 

.  1  . 

(X10) 


where  W  is  the  received  intensity  of  polarization. 

If  a  receiver  is  able  to  measure  the  phase  between  the  V  and  H 
channels,  it  is  possible  to  do  the  incoherent  method  without 
measuring  6  polarizations.  This  method  is  called  the 
coherent-on-receive  (COR)  technique.  The  elements  of  the  Stokes 
vector  can  be  expressed  as 
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where  6  is  the  phase  difference  between  V  and  H  channels. 
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then  (A. I)  in  terms  of  Stokes  vector  becomes 
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Radara 


Figure  1  Block  Diagram  o£  the  MMW  Polaximelric  Radar, 


Figure  3  Photograph*  of  the  Metallic  Tree  mnd  Rhododendron  Tree, 


35  GHz  Radar  (Fully  Polarimetrlc) 


Transmitter 
IF  bandwidlh 
Power 
Antenna 
Polarization 


up  to  2  GHz 
-1-23  dBm 

6‘  Lens  fbeamwldth  4.2  degrees) 
Any  polarization 


Receive  r 

Dual  Channel 

Mixers 

Antenna 


Vend  H 

Fundamental  mixing 

6“  Lens  (beamwldih  4.2  degrees) 


Polarimetrlc  data 

Incoherent  (coherent-on-recelve) 
Mueller  matrix 


Coherent 


Scattering  amplitude  matrix 

P&L 


Rhododendron  Tree  (Target  in  motion) 

Incoherent  Polarimetric  Measurement  Mode  (500  samples) 
Fast  Motion 


1 

0.16  -0.006 

0.007  " 

0.186 

0.828  -0.017 

0.017 

0.04 

0.059  0.735 

0.056 

-0.059 

-0.023  -0.019 

0.472  _ 

Slow  Motion 

1 

0.231  -0.052 

0.012  ' 

0.126 

1.035  -0.021 

0.049 

-0.038 

0.082  0.697 

0.003 

.  -0.013 

-0.015  -0.064 

0.619 

Coherent  Polarimetric  Measurement  Mode  (8000  samples) 
"  1  0.159  -0.002  -0.006' 

0.179  0.823  -0.003  -0.018 
0.0  -0.01  0.683  -0.023 

-0.033  -0.001  0.003  0.596 


Metallic  Tree  (Target  in  motion) 

Incoherent  Polarimetric  Measurement  Mode  (500  samples) 
Fast  Motion 


1 

0.089 

0.02 

0.03  ' 

0.094 

0,74 

0.02 

0.001 

0.026 

0.011 

0.888 

0.105 

: 

-0.028 

-0.007 

-0.126 

0.619 

Coherent  Polarimetric  Measurement  Mode  (8 

'  1 

0,06 

0.0 

0.0  ' 

0.06 

1.16 

0.009 

0.004 

0.011 

0.002 

0.973 

-0.053 

0.008 

0.002 

0.049 

0.89 

Table  1.  Average  Mueller  matrices  of  rhododendron  and 

metallic  trees  measured  with  coherent  and  incoherent 
polarimetric  measurement  modes. 
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Figure  4  Probability  Density  Function  Measured  by  the  Coherent  Figure  5  Probability  Density  Function  Computed  from  the  Average 

Polarimetric  System.  Mueller  Matrix. 


DISCUSSION 


G.  Neininger,  Germany 

What  is  the  realized  cross  polarization  isolation? 

Author's  Reply 

I  am  not  sure  but  I  believe  that  it  is  listed  in  the  preprint  article. 

E.  Schweicher,  BE 

You  mentioned  (in  your  second  talk)  a  solid  state  design.  Does  that  mean  a  solid-state 
transmitter  and  in  that  case  what  kind  of  transistors  are  used  at  36  GHz  and  94  GHz? 

Author’s  Reply 

The  presenter  did  not  know. 

G.  Brown,  US 

The  histograms  of  the  phase  differences  between  the  co-  and  cross  polarization  returns  for 
data  (from  real  and  metallic  trees)  were  relatively  narrow  in  spread.  What  range  bin  was  used  to 
generate  these  data  and  how  do  you  explain  the  rather  small  phase  difference  variance? 

Author's  Reply 

I  do  not  know  the  range  bin  used  in  the  phase  distribution  plots  since  I  did  not  work  on 
those  data.  However,  I  do  know  that  a  very  large  number  of  samples  were  used  in  the  preparation 
of  the  phase  statistics,  so  I  imagine  that  the  bin  was  rather  small. 


POLARIMETRIC  OBSERVATIONS  OF  TREES  AT  35  AND  94  GHZ 

Adib  Nashashibi*,  Yasuo  Kuga,  and  Fawwaz  T.  Ulaby 
Radiation  Laboratory 

Department  of  Electrical  Engineering  and  Computer  Science 
University  of  Michigan,  Ann  Arbor,  MI  48 109 


l.  SUMMARY 

Two  different  kinds  of  trees  were  measured  with  a  newly  developed 
polarimetric  radars  at  35  and  94  GHz.  The  system  is  based  on  the  COR  technique 
and  is  capable  of  measuring  the  Mueller  matrix  directly.  The  degree  of  polarization, 
phase  statistics,  and  polarization  signatures  are  presented  in  this  paper.  Our  results 
show  that  the  degree  of  polarization  is  sensitive  to  the  radar  frequencies  and  tree 
types. 

II.  SYSTEM  DESCRIPTION 

The  block  diagram  of  the  University  of  Michigan  MMW  polarimetric  radar  is 
shown  in  Fig.l.  The  35  and  94  GHz  radars  are  fully  polarimetric  and  are  capable  of 
measuring  the  Mueller  matrix  using  the  coherent-on-receive  (COR )  technique  [2,3]. 
With  the  COR  technique,  6  different  polarizations  (V,  H,  45,  135,  LHC,  and  RHC) 
are  transmitted  sequentially,  and  the  phase  and  magnitude  of  V  and  H  channels  are 
received  simultaneously.  Unlike  the  coherent  polarimetric  radar  which  measures  the 
scattering  amplitude  matrix  S,  the  COR  radar  measures  the  Mueller  matrix  Lm. 

The  original  MMW  system  was  based  on  the  HP8510A  and  it  did  not  have  the 
capability  of  obtaining  data  from  two  different  channels  simultaneously.  To 
overcome  this  problem,  a  delay  line  technique  was  developed.  With  this  technique, 
one  channel,  usually  H,  is  delayed  by  a  long  delay  line  (200  nsec  delay)  and 
combined  with  the  non-delayed  channel.  The  received  data  in  the  time  domain 
shows  V  and  H  channels  separated  by  the  electrical  length  of  the  delay  line.  Using 
the  time  gating  function  of  the  network  analyzer,  V  and  H  channels  can  be 
processed  separately. 

The  calibration  of  the  COR  system  is  performed  in  two  steps.  First,  the  receiver 
is  calibrated  using  an  odd  bounce  reflector  and  a  polarizing  grid  placed  in  front  of  the 
receiving  antenna.  The  receiver  distortion  matrix  is  obtained  by  positioning  the 
polarizer  at  three  different  angles.  Second,  the  actual  transmitted  polarizations  are 
obtained  using  the  calibrated  receiver  and  the  odd  bounce  reflector  [3]. 

m.  MUELLER  MATRIX  AND  DEGREE  OF  POLARIZATION 

The  polarization  of  the  EM  wave  can  be  expressed  by  the  modified  Stokes 
vector,  Fn,  which  is  defined  in  terms  of  the  vertical,  Ey,  and  horizontal,  E^,  electric 
fields  [1]. 

Fm=[Il,  I2.  U,  V]T,  I1=IEvl2,I2»IEhl2tU=2Re(EvEh*),andV=2Im(EvEh*)  (1) 
The  transmitted  and  received  modified  Stokes  vectors  are  related  to  each  other  by 

F„r  =  (1/r2)  L„  Fnl  (2) 
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where  r  is  the  distance  between  the  target  and  the  receiving  antenna,  Lm  is  the 
modified  Mueller  matrix  of  the  target,  and  Sik  (l,k=v  or  h)  are  the  four  elements  of 
the  scattering  matrix  5. 

The  COR  radar  obtains  Lm  directly  by  sending  4  to  6  independent 
polarizations.  Using  the  elements  of  Lm,  the  phase  difference  of  the  co-  and  -cross- 
polarized  components  of  S  can  be  written  as 

$vv  -hh  =  tan-1[(L43-L34)/(L33+L44)]  (4) 

<|>w  -hv  =  tan’^Lai/Lai)  (5) 

4>vh-hv=  tan-1[(L43+L34)/(L33-L44)],  (6) 


The  degree  of  polarization,  m,  which  is  directly  obtained  from  the  received 
Fm,  is  given  by 


m—  V((<Ijp>  -  <l2p>)2+<Up>2+<Vp>2)  /  (<Ijp>+<l2p>)  (7) 

whereo  is  an  ensemble  average,  and  "P”  denotes  the  transmitted  polarization. 

IV.  RESULTS 

Measurements  were  conducted  using  leafy  evergreen  and  coniferous  trees.  The 
first  was  a  rhododendron  (rhodo)  with  a  planophil  leaf  orientation  distribution  and  an 

average  leaf  area  of  40  cm2.  The  second  was  a  spruce  blue  shiner  (spruce)  with  a 
uniform  leaf  orientation  distribution  and  an  average  needle  length  of  1.5  cm.  To 
obtain  the  statistical  data,  more  than  100  independent  samples  were  measured  for 
each  incidence  angle. 

Table  1  shows  the  degree  of  polarization  for  6  different  incident  polarizations  at 
35  and  94  GHz.  The  degree  of  polarization  is  quite  sensitive  to  frequencies  and  tree 
types.  The  MMW  return  from  spruce  tree,  for  example,  shows  a  significant  change 
in  m  between  35  and  94  GHz.  The  dependence  on  the  incidence  angle  is  also  clear 
for  the  planophil  type  tree. 


Figure  2  shows  the  phase  distribution  of  the  received  signals.  In  all  cases,  <f>w 
.hv  had  a  uniform  distribution  showing  the  decorrelation  of  co-  and  cross-pol. 
components.  The  most  interesting  case  is  the  phase  difference  between  the  two  co¬ 
polarized  returns  shown  in  Fig.  2.  The  distribution  of  <t>w  _hh  seems  to  depend  on 
the  tree  type  and  on  the  incidence  angle.  $w  -hh  rhododendron  has  a  much 

narrower  spread  than  that  of  spruce.  The  phase  <JVh  -hv  °f  a  troty  monostatic  system 
should  be  a  delta  function  centered  at  zero  degrees.  Our  data,  however,  shows  a 
wide  distribution.  We  think  this  is  caused  by  the  bistatic  configuration  of  the  radar 
system . 


I 

1 


I 
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Figure  3  is  the  co-polarized  polarization  signatures  of  rhododendron  and 
spruce.  The  pedestal  which  is  related  to  depolarization  is  much  larger  for  spruce  tree 
showing  a  strong  depolarization. 
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TABLE  1.  Degree  of  polarization  m  of  rhodo  and  spruce  trees  at  35  and  94  GHz. 


mm 
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0.83 

0.87 

0.69 

0.49 

0.49 

0.82 

0.77 
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0.77 
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0.1 
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0.81 

0.73 

RHC 
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0.79 
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0.06 
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0.77 

0.73 
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Figure  1:  Block  diagram  of  ttm  COR  systi 
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igure  3:  Co-polarized  signatures  at  35  GHz  and  90°  incidence 
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SUMMARY 

.A*1  extensive  radar  clutter  database  was  generated  by  the 
University  of  Michigan’s  millimeter-wave  mobile  polarimetric 
radar  system  during  the  past  few  years  (1).  The  data  base 
includes  millimeter-wave  observations  of  snow,  trees, 
vegetation,  and  soil  and  road  surfaces  at  35, 94, 140  and  215 
GHz.  The  radar  measurements  were  often  augmented  with 
close-up  observations  of  the  targets  including  such 
measurement  as  water  contents  and  surface  roughness,  when 
appropriate.  For  each  data  set,  a  summary  of  these 
observations  and  photographs  of  the  target  scene  are  provided. 
The  millimeter-wave  system  consists  of  track-mounted  radars 
capable  of  making  observations  from  a  20-m  high  platform  at 
incidence  angles  between  0  and  70  degrees.  The  35  and  94 
GHz  radars  are  fully  polarimetric  and  capable  of  obtaining  the 
Mueller  matrix  in  situ. 

MEASUREMENT  SYSTEMS 


Radar  System 

The  block  diagram  of  the  MMW  radar  system  at  the 
University  of  Michigan  is  shown  in  Figure  1.  The  system 
consists  of  MMW  radar  front-ends  at  35,  94,  140,  and  215 
GHz,  a  receiver,  and  a  system  control  unit  The  receiver  is 
based  on  the  HP8753  vector  network  analyzer  (NWA).  The 
IF  signal  (2  to  3  GHz)  from  the  NWA  is  upconverted  to  the 
MMW  frequency  in  each  radar  front-end  and  the  received 
MMW  signal  is  down-converted  to  the  IF  signal  and  fed  into 
the  NWA. 


Ridin 


Fig.  1 .  Block  diagram  of  MMW  radar  system. 


Both  software  and  hardware'  gatiags  are  utilized  to 
eliminate  the  unwanted  noise  from  the  returned  signal. 
Because  the  system  has  one  GHz  bandwidth,  the  software 
gating  using  the  time  domain  capability  of  NWA  gives  an 
approximately  15  cm  spatial  resolution.  Unfortunately,  it 
takes  400  millisecond  to  sweep  401  frequency  points  on 
HP8753.  If  a  target  is  moving,  the  sweep  time  can  be  much 
longer  than  the  target  decoirelation  time  at  the  MMW 
frequency.  The  software  gating,  therefore,  is  mainly  used  for 
calibration  and  indoor  measurements  in  which  the  target 
motion  is  negligibly  small.  For  outdoor  tree  measurements,  a 
hardware  gating  circuit,  which  creates  20  nsec  pulses  with  a 
pulse-repetition-rate  of  5  MHz,  is  used.  Because  the 
pulse-repetition-rate  is  much  greater  than  the  NWA  IF 
bandwidth  (3KHz),  the  system  is  essentially  operating  in  the 
CWmode. 

■Hie  original  MMW  system  was  based  on  the  HP8510A 
and  it  did  not  have  the  capability  of  obtaining  data  from  two 
different  channels  simultaneously.  To  overcome  this  problem, 
a  delay  line  technique  was  developed.  With  this  technique, 
one  channel,  usually  H,  is  delayed  by  a  long  delay  line  (200 
nsec  delay)  and  combined  with  a  non-delayed  channel.  The 
received  data  in  the  time  domain  shows  V  and  H  channels 
separated  by  the  electrical  length  of  the  delay  line.  Using  the 
time  gating  function  on  the  NWA,  V  and  H  channels  can  be 
processed  separately. 

The  data  acquisition  and  the  system  control  are  performed 
by  a  HP  computer  with  a  control  program  written  in  HP 
BASIC.  An  IBM  AT-compatible  computer  controls  the 
polarization  control  unit  which  can  be  controlled  by  HP 
computer  in  the  automatic  mode  through  a  RS232  interface. 


35  and  94  Radar  Front-ends 


The  35  and  94  GHz  radars  are  fully  polarimetric  and 
capable  of  measuring  the  Mueller  matrix  using  the 
coherent-on-rcceive  (COR)  technique  [2,3].  With  the  COR 
technique,  6  different  polarizations  (V,  H,  45,  135,  LHC,  and 
RHC)  are  transmitted  sequentially,  and  the  phase  and 
magnitude  of  V  and  H  channels  are  received  simultaneously. 
Unlike  the  coherent  polarimetric  radar  which  measures  the 
scattering  amplitude  matrix  5,  the  COR  radar  measures  the 
Mueller  matrix  Lm.  The  description  of  COR  technique  is 
given  in  Appendix  I.  The  block  diagram  of  the  94  GHz  radar 
is  shown  in  Figure  2.  The  transmitted  wave  polarization  is 
controlled  by  two  independently  rotatable  quarter-wave  plates. 
A  precise  control  of  wave  plate  motion  is  provided  by  stepping 
motors  which  are  controlled  by  an  IBM  AT-compatible 
computer.  To  stabilize  the  MMW  LO,  a  stable  X-band  signal 
is  multiplied  and  injected  into  the  Gunn  oscillator.  In  our 
original  design,  the  upconverted  MMW  signal  was  transmitted 
without  amplification,  resulting  in  poor  S/N  ratio.  To  improve 
the  S/N  ratio  and  to  obtain  stable  phase  data,  MMW  amplifiers 
were  added  to  boost  the  output  power  from  +4  dBm  to  +23 
dBm  at  35  GHz  and  from  +0  dBm  to  +13  dBm  at  94  GHz. 


Fig.  2.  Block  diagram  of  94  GHz  radar.  35  GHz 
radar  has  a  similar  design. 


140  and  215  GHz  Radar  Front-ends 

The  present  140  and  215  GHz  radars  are  limited  toco-  and 
cross-pol.  measurements  only.  The  block  diagram  is  shown  in 
Figure  3.  The  polarization  of  transmitted  and  received  signals 
are  controlled  by  wave  plates  and  the  receiver  is  limited  to 
measuring  either  V  or  H  channel  at  a  time.  Because  Gunn 
oscillators  are  not  available  at  these  frequencies,  the  MMW 
signal  is  obtained  by  multiplying  a  lower  frequency  signal  by  2 
at  215  GHz  and  by  3  at  140  GHz.  The  transmitted  power  is  -4 
dBm  at  140  GHz  and  -10  dBm  at  215  GHz. 
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Fig,  3.  Block  diagram  of  140  and  215  GHz  radars. 


EXFER1 


rALDATA 


In  the  following  section  we  will  describe  two  data  sets 
obtained  with  our  MMW  system.  The  first  is  a  snow  diurnal 
data  set  obtained  in  March  1990  and  the  second  data  set  shows 
the  polarimetric  observations  of  two  kinds  of  trees  conducted 
recently. 

Snow  Diurnal  Data 


Because  of  absorption  by  water,  MMW  backscattering 
from  snow  is  quite  sensitive  to  liquid  water  content  (LWC), 
and  this  sensitivity  depends  on  the  radar  frequency.  When 
snow  is  frozen,  the  scattering  by  the  snow  ice  particles  is 
significant  and  consequently,  the  backscattering  cross-section 
is  large.  As  snow  melts  during  the  warm  daytime  hours  and 
LWC  increases,  the  radar  backscattering  cross-section 
decreases  dramatically.  Figure  4  shows  the  backscattering 
coefficient  at  35,  94,  and  140  GHz  as  a  function  of  time. 
Figure  5  is  the  air  temperature  and  LWC  measured  by  a 
freezing  calorimeter  at  the  same  time.  The  sharp  increase  of 
LWC  in  snow  is  clearly  detected  by  MMW  radars.  Among  the 
three  channels,  the  140  GHz  data  is  least  sensitive  to  LWC.  A 
detailed  theoretical  analysis  and  more  extensive  experimental 
results  are  given  in  References  4  and  5. 
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Fig.  4.  Backscattering  variation  measured  at  35, 94, 
and  140  GHz  on  1  March,  1990.  The  Incident  angle 
was  40*. 


Fig.  5.  Measured  air  temperature,  snow  surface 
temperature,  and  liquid  water  content  of  the  top  5-cm 
layer  on  1  March,  1990. 


Folarimciric  ObgCTVflrinq;i  of  Tup* 


MMW  scattering  by  tree  is  usually  dominated  by  the  tree 
canopy  during  the  spring  and  summer  time.  The  polarization 
ot  the  scattered  wave  may  contain  information  which  can  be 
used  for  extracting  tree  characteristics  such  as  leaf  size  and 
density.  Hie  35  and  94  GHz  radars  at  the  University  of 
Michigan  are  capable  of  obtaining  fully  polarimetric 

information  including  the  Mueller  matrix,  phase  statistics,  and 
degree  of  polarization. 


Measurements  were  conducted  for  leafy  evergreen  and 
coniferous  trees.  The  first  was  a  rhododendron  (rhodo)  with  a 
pianophil  leaf  orientation  distribution  and  an  average  leaf  area 
of  40 cm  •The  second  was  a  spruce  blue  shiner  (spruce)  with 
a  uniform  leaf  orientation  distribution  and  an  average  needle 
length  of  1.5  cm.  To  obtain  the  statistical  data,  more  than  100 
independent  samples  were  measured  for  each  incidence  angle. 


Table  1  shows  the  degree  of  polarization  for  6  different 
incident  polarizations  at  35  and  94  GHz.  The  degree  of 
is  quite  sensitive  to  frequency  and  tree  type.  The 
MMW  return  from  the  spruce  tree,  for  example,  shows  a 
significant  change  in  m  between  35  and  94  GHz.  The 
dependence  on  incidence  angle  is  also  clear  for  the  pianophil 
type  tree. 


TABLE  1.  Degree  of  polarization  i 

m  of  rhodo  and  spruce  trees 

at  35  and  94  GHz, 

POL 

35  GHz 

94  GHz 

RHODO 

SPRUCE 

RHODO 

SPRUCE 

30° 

50° 

90° 

30° 

90° 

50“ 

90° 

9011 

V 

0.86 

0.89 

0.80 

0.47 

0.60 

0.92 

0.82 

0.64 

45 

0.80 

0.85 

0.66 

0.48 

-  0.43 

0.91 

0.87 

0.71 

H 

0.76 

0.83 

0.69 

0.37 

0.57 

0.87 

0.87 

0.71 

135 

0.83 

0.87 

0.69 

0.49 

0.49 

0.82 

0.77 

0.58 

LHC 

0.69 

0.77 

0.61 

0.1 

0.11 

0.81 

0.73 

0.39 

RHC 

0.61 

0.79 

0.51 

0.06 

0.21 

0.77 

0.73 

0.39 

Figure  6  shows  the  phase  distribution  of  the  received 
signals.  In  alt  cases,  $yy  had  a  uniform  distribution, 
suggesting  that  the  co-  and  cross-pol.  components  arc 
decorrelated.  Tie  most  interesting  case  is  the  phase  difference 
between  the  two  copolarized  returns  shown  in  Fig.  6.  The 

distribution  of  ^  ^  seems  to  depend  on  tree  type  and  on 
the  incidence  angle.  $vv  of  rhododendron  has  a  much 

narrower  spread  than  that  of  spruce.  The  phase  4>vh  -hv  of  a 
truly  monostatic  system  should  be  a  delta  function  centered  at 
zero  degrees.  Our  data,  however,  shows  a  wide  distribution. 
We  think  this  is  caused  by  the  quasi-bistatic  configuration  of 
tne  radar  system  (small  separation  between  transmit  and 
receive  antennas). 


Fig.  6.  Distributions  of  ^  _hh  at  35  and  94  GHz. 


Figure  7  is  the  co-polarized  polarization  signatures  of 
rhododendron  and  spruce.  The  pedestal,  which  is  related  to 
depolarization,  is  much  larger  for  the  spruce  tree  than  for  the 
rhododendron  tree. 


Fig.  7.  Co-polarized  signatures  at  35  GHz  and  90’ 
incidence  angle. 


CONCLUSION 


MMW  technology  has  progressed  significantly  in  recent 
years  and  the  cost  of  making  MMW  radars  has  been 


decreasing.  To  utilize  MMW  radars  for  target  detection  and 
identification  it  is  important  to  understand  how  millimeter- 
waves  interact  with  geophysical  media.  We  believe  the  M MW 
Clutter  Handbook  published  at  the  University  of  Michigan  will 
be  a  useful  source  of  data  for  engineers  and  scientists  working 
with  MMW  radars  [1]. 
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APPENDIX.  I 

CQHERENT-QN-RECEIVE  TECHNIQUE 

The  complete  polarimetric  data  can  be  obtained  by 
measuring  either  the  scattering  amplitude  matrix  S  or  the  Mueller 
matrix  Lm-  The  measurement  of  the  scattering  amplitude  matrix 
requires  accurate  phase  measurement  Also  4  elements  of  S  must 
be  obtained  within  the  decorrelation  time  of  the  target  which  is  in 
the  order  of  millisecond  at  the  MMW  frequency.  The  scattered 
electric  field  Er  in  terms  of  the  scattering  amplitude  matrix  S  and 
the  incident  electric  field  is  given  by 


To  obtain  5,  wc  need  to  send  [Ey,0]  and  [0,Eh],  anc* 
measure  Ev  and  Ejt  simultaneously*  Because  a  long  time  is 
required  to  sweep  the  freouency  points,  the  scatterometer  based 
on  the  NWA  is  not  suited  for  the  fast  acquisition  of  scattering 
matrix.  For  example,  sweeping  401  points  takes  400  msec  and  if 
we  measure  VV,  HV,  VH,  and  HH  sequentially,  it  takes  at  least 
2  seconds  to  acquire  4  elements  Of  S.  Obviously  this  technique 
does  not  work  at  the  MMW  frequency. 

The  polarized  wave  can  also  be  expressed  in  terms  of  the 
Stoke's  vector  Fm  which  is  defined  as  [6] 


Fto=Ui.  I*  U.  V]T.  Ii=IEv  I2.  l2='Eh  1 2-  U=2Rc(EvEh*), 
and  V=2Im(EvEh*)  A*3 

Then  (A-l)  in  terms  of  Stoke’s  vector  becomes 


Fmr  =  (1/t2)  Fm‘ 


A-4 


Lm 


\s,y  is,j* 

|5J'  IS., I1  R«<p..> 

2Re(5„0  2R,ts-0  “«<*.. C  * 

2 1  m(5„5  j  2 1  I  m(S ,  ,S*»  +  5,.0 


•In(S,X.  •  S,/j 


A-5 


where  Lm  is  cnlled  a  Mueller  matrix. 


The  totally  incoherent  method  does  not  require  the  phase 
measurement  With  this  method  the  ^elements  of  Stoke's  vector 
are  obtained  by  receiving  6  polarizations  (V,  H,  45,1J3,  LnU 
RHC)  simultaneously.  For  example,  the  third  element  of  Stoke  s 
vector  U  is  given  as  a  ratio  of  intensities  at  45  linear  to  135  linear. 
For  a  given  incident  polarization,  therefore,  we  can  obtain  a 
column  of  the  Mueller  matrix.  To  get  the  complete  Mueller 
matrix,  we  need  to  repeat  this  process  for  4  independent  incident 
polarizations.  Altogether  at  least  24  magnitude  only 
measurements  are  required  to  obtain  the  complete  Mueller 
Although  the  phase  measurement  is  not  required  with  the 
incoherent  method,  it  is  necessary  to  receive  all  6  polarizations 
simultaneously.  The  Stoke's  vector  in  terms  of  6  polanzauons 
and  a  set  of  4  independent  incident  polarizations  are  shown 
below. 


h 


W, 

"  W„  +  Wy 


A-6 


_  W„ 

”  Wh  +  Wv 


A-7 


W4S  -  W|35 

W45  +  W,35 


A-8 


V  Wmc  •  Wrhc 
Wine  +  Wbhc 


A-9 


V 

V 

'V 

'V 

0 

0 

.  • 

1 

0 

.  Us  - 

vz 

1 

.  1  LHC  - 

V2 

0 

.0, 

.0. 

0 

1 

A- 10 


where  W*  is  a  received  intensity  of  *  polarization. 

If  a  receiver  is  able  to  measure  phase  between  V  pud  H 
channels,  it  is  possible  to  do  the  incoherent  method  without 
measuring  6  polarizations.  This  method  is  called  the 
coherent-on-rcceive  (COR)  technique  [2,3].  The  elements  of  the 
Stoke's  vector  can  be  expressed  as 


li  -|E»| 
li  - 1  Ewl a 

U»2|E»(lEh|cos8 
V  ■  2 1  Ev||  Eh|  sin  5 


A-l  1 
A-12 
A-13 
A-14 


where  8  is  a  phase  difference  between  V  and  H. 

The  requirements  of  the  COR  method  arc  (1)  the 
transmitter  must  be  able  to  transmit  at  least  4  independent 
polarizations  and  (2)  the  receiver  must  be  able  to  get  V  and  H 
channels  simultaneously.  The  measurement  time  between  the 
four  incident  polarizations  can  be  much  longer  than  the 
decorrelation  time  of  target  because  correlation  between  V  and  H 
is  obtained  directory.  Since  it  is  relatively  easy  to  obtain  the 
phase  difference  between  V  and  H  at  the  MMW  frequencies,  we 
chose  the  COR  method  for  our  radars.  Also  COR  is  suited  for  a 
radar  system  based  on  NWA  as  a  receiver. 


Using  the  elements  of  Lm,  the  phase  difference  of  the  co-  and 
cross-polarized  components  of  5  can  be  written  as 


-hh  *  tan*l[(U3-L34)/(L33+L44)] 

A-15 

•frw  -hv  ~  *an  ^(L4l/L3l) 

A- 16 

<t>vh-hv=  tan*l[(L43+L34)/(L33-L44)I, 

A- 17 

The  degree  of  polarization,  m,  which  is  directly  obtained 
from  the  received  Fm,  is  given  by 

trt=  V((<I1P>  -  <l2P>)2+<Up>^+<Vp>2)  /  (<I]p>+<l2p>)  A- 18 

whereo  is  an  ensemble  average,  and  "P"  denotes  the  transmitted 
polarization. 
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Abstract- 

To  answer  the  question  of  what  radar  polarimetry  has  to  offer  to  the  remote  sensing 
of  random  media,  statistics  of  the  phase  difference  of  the  scattering  matrix  elements 
must  be  studied.  Recent  polarimetric  measurements  of  rough  surfaces  have  indicated 
that  the  statistical  parameters  of  the  phase  difference  (mean,  standard  deviation,  etc.) 
are  very  sensitive  to  some  of  the  physical  parameters.  In  this  paper  the  probability 
density  function  of  the  phase  differences  is  derived  from  the  Mueller  matrix  assuming 
that  the  elements  of  the  scattering  matrix  are  jointly  Gaussian.  It  is  shown  that  the 
probability  density  functions  of  the  co-polarized  and  cross- polarized  phase  differences 
are  similar  in  form  and  each  can  be  determined  by  two  parameters  (a  and  £)  completely. 
The  expressions  for  the  probability  density  functions  axe  verified  by  comparing  the  his¬ 
tograms,  the  mean,  and  the  standard  deviations  of  phase  differences  derived  directly 
from  polarimetric  measurements  of  variety  of  rough  surfaces  to  the  probability  density 
function  and  its  mean  and  standard  deviation  derived  from  the  Mueller  matrices  of  the 
same  data.  The  expressions  for  the  probability  density  functions  are  of  special  interest 
for  non-coherent  polarimetric  radars  and  non-coherent  polarimetric  models  for  random 
media  such  as  vector  radiative  transfer. 
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1  Introduction 


In  the  past  decade  substantial  effort  within  the  microwave  remote  sensing  community  has 
been  devoted  to  the  development  and  improvement  of  polarimetry  science.  Polarimetric 
radars  are  capable  of  synthesizing  the  radar  response  of  a  target  to  any  combination  of 
the  receive  and  transmit  polarizations  from  coherent  measurements  of  the  target  with  two 
orthogonal  channels.  Polarimetric  radars  have  demonstrated  their  abilities  in  improving 
point-target  detection  and  classification  (Ioannidis  and  Hammers,  1979].  That  is,  for  a 
point  target  in  a  clutter  background  the  transmit  and  receive  polarizations  can  be  chosen 
such  that  the  target  to  clutter  response  is  maximum.  Also,  different  point  targets  in  the 
radar  scene  can  be  classified  according  to  their  optimum  polarization.  Although  radar 
polarimeters  have  shown  a  great  potential  in  point-target  detection  and  classification, 
their  capabilities  in  remote  sensing  of  distributed  targets  is  not  completely  understood 
yet. 

Considering  the  complexity  involved  in  designing,  manufacturing,  and  processing  the 
data  of  an  imaging  polarimeter  as  opposed  to  a  conventional  imaging  radar,  it  is  neces¬ 
sary  to  examine  the  advantages  that  the  imaging  polarimeter  provides  about  the  targets 
of  interest.  For  sample,  in  retrieving  the  biophysical  parameters  from  the  polarimetric 
radar  data  one  should  ask  whether  there  exists  a  dependency  between  the  parameters 
and  the  measured  phase  of  the  scattering  matrix  components.  If  the  answer  is  nega¬ 
tive,  obviously  gathering  polarimetric  data  for  inversion  of  that  parameter  is  a  waste 
of  effort.  One  way  of  confirming  this  question  is  by  collecting  data  over  a  range  of  the 
desired  parameter  while  keeping  other  influential  parameters  constant.  This  procedure, 
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if  not  impossible,  is  very  difficult  to  conduct  because  of  problems  in  repeatability  of 
the  experiment  and  difficulties  in  controlling  the  environmental  conditions.  Moreover  at 
high  frequencies  (millimeter- wave  frequencies  and  higher)  coherent  measurement  of  the 
scattering  matrix  is  impossible  because  of  instabilities  of  local  oscillators  and  relative 
movements  of  the  target  and  the  radar  platform  [Meads  and  McIntosh,  1991].  At  these 
frequencies  non-coherent  radars  are  employed  which  provide  the  Mueller  matrix  of  the 
target. 

Another  approach  to  examine  the  dependency  of  the  radar  response  to  the  desired 
parameters  of  the  targets  is  the  application  of  theoretical  models.  One  of  the  most 
successful  polarimetric  models  for  random  media  is  the  vector  radiative  transfer  theory 
[Tsang  et  al.,  1985].  This  model  is  based  on  conservation  of  energy  and  the  single 
scattering  properties  of  the  constituent  particles.  The  solution  of  the  radiative  transfer 
equation  relates  the  scattered- wave  Stokes  vector  to  the  incident-wave  Stokes  vector  via 
the  Mueller  matrix.  The  Mueller  matrix,  as  computed  by  this  method,  is  an  ensemble- 
averaged  quantity  because  of  the  inherent  nature  of  the  radiative  transfer  theory.  Since 
the  Mueller  matrix  is  related  to  the  scattering  matrix  through  a  nonlinear  process  and 
the  components  of  the  scattering  matrix  are  statistically  dependent,  the  information 
about  the  phase  difference  of  the  scattering  matrix  components  cannot  be  obtained  from 
the  Mueller  matrix  directly.  To  achieve  information  about  the  phase  statistics,  one 
may  resort  to  the  Monte  Carlo-type  models  which  are  computationally  inefficient  and  in 
general  inaccurate. 

Experimental  observations  of  phase  difference  statistics  from  a  polarimetric  SAR  at 
L-band  [Ulaby  et  al.,  1987;  Zebker  et  al.,  1987]  over  agricultural  terrain  and  bare  soil 
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surfaces  indicate  that  the  statistics  of  the  co-polarized  phase  difference  depends  on  the 
target  type  and  its  conditions.  Recent  measurements  of  bare  soil  surfaces  by  polarimetric 
scatterometers  show  that  the  variance  of  the  co-polarized  phase  difference  is  a  function 
of  the  roughness  parameters  and  incidence  angle  but  is  less  sensitive  to  moisture  content 
[Sarabandi  et  al.,  1991]. 

In  view  of  difficulties  in  measuring  the  scattering  matrix  at  high  frequencies  and 
performing  controlled  experiments,  it  is  necessary  to  establish  a  relationship  between 
the  Mueller  matrix  and  the  statistics  of  the  phase  differences  of  the  scattering  matrix 
elements.  In  the  next  section  we  derive  the  probability  density  function  of  the  co-  and 
cross- polarized  phase  difference  in  terms  of  the  Mueller  matrix  elements  assuming  that 
the  scattering  matrix  elements  are  jointly  Gaussian.  Then  the  assumptions  and  final 
results  are  compared  with  the  experimental  data  acquired  by  polarimetric  scatterometers 
in  Section  3. 

2  Theoretical  Derivation  of  Phase  Difference  Statistics 


The  polarimetric  response  of  a  point  or  distributed  target  can  be  obtained  by  simultane¬ 
ously  measuring  both  the  amplitude  and  phase  of  the  scattered  field  using  two  orthogonal 
channels.  If  the  incident  and  scattered  field  vectors  are  decomposed  into  their  horizontal 
and  vertical  components,  the  polarimetric  response  can  be  represented  by  the  scattering 
matrix  S,  which,  for  plane  wave  illumination  we  can  write 


E' 


Bikr 


Svv  $vh 
Shv  Sfxh 


E* 


(1) 
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where  r  is  the  distance  from  the  radar  to  the  center  of  the  distributed  target.  It  should 
be  noted  that  in  the  backscattering  case  reciprocity  implies  Svk  =  Each  element  of 
the  scattering  matrix,  in  general  is  a  complex  quantity  characterized  by  an  amplitude 
and  a  phase.  When  the  radar  illuminates  a  volume  of  a  random  medium  or  an  area  of  a 
random  surface,  many  point  scatterers  contribute  to  the  total  scattered  energy  received 
by  the  radar  and  therefore  each  element  of  the  scattering  matrix  may  be  represented  by 

s„  =  \s„\e<*"  =  £  KW*”  P,q  =  v,h  .  (2) 

n  =  l 

Here  N  is  the  total  number  of  scatterers  each  having  scattering  amplitude  \spq\  and 
phase  <f>pq.  It  should  be  mentioned  that  the  phase  of  each  scatterer,  as  given  in  (2), 
includes  a  phase  delay  according  to  the  location  of  the  scatterer  with  respect  to  the 
center  of  the  distributed  target.  Without  loss  of  generality  all  multiple  scattering  over 
the  surface  or  in  the  medium  can  be  included  in  (2).  Since  the  location  of  the  scatterers 
within  the  illuminated  area  (volume)  is  random,  the  process  describing  the  phasor  spq 
is  a  Wiener  process  (random  walk)  [Davenport,  1970].  If  N  is  large  enough,  application 
of  the  central  limit  theorem  shows  that  the  real  and  imaginary  parts  of  the  scattering 
matrix  element  Spq  are  independent  identically  distributed  zero-mean  Gaussian  random 
variables.  Equivalently  it  can  also  be  shown  that  |5P9|  and  4>pq  are,  respectively,  Rayleigh 
and  uniform  independent  random  variables.  The  three  elements  of  the  scattering  matrix, 
in  general,  can  be  viewed  as  a  six-element  random  vector  and  it  is  again  reasonable  to 
assume  that  the  six  components  are  jointly  Gaussian. 

Observation  of  polarimetric  data  for  a  variety  of  distributed  targets  such  as  bare 
soil  surfaces  and  different  kinds  of  vegetation-covered  terrain  all  indicate  that  the  cross- 
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polarized  component  of  the  scattering  matrix  (5/,v)  is  statistically  independent  of  the  co¬ 
polarized  terms  ( Svv  and  Shh).  Therefore  the  statistical  behavior  of  5/,v  can  be  obtained 
from  a  single  parameter,  namely  the  variance  ( a *)  of  the  real  or  imaginary  part  of 
Skv  =  +  iX6,  that  is 

fxs.x<(*  5.*«)  =  ^exp[-^2^p] 

or  equivalently  the  joint  density  function  \Svh\  and  <t>vh  Is 

=  2k<j2  !6^lexp(“  2a2  ^  ^ 

which  indicates  that  <f>vh  is  uniformly  distributed  between  (  —  7r,+7r). 

Since  measurement  of  the  absolute  phase  of  the  scattering  matrix  elements  is  very 
difficult,  it  is  customary  to  factor  out  the  phase  of  one  of  the  co-polarized  terms,  for 
example  Svvy  and  therefore  the  phase  difference  statistics  are  of  concern  as  opposed  to 
the  absolute  phases.  Since  S^v  is  assumed  to  be  independent  of  Svv  (not  a  necessary 
assumption)  and  both  <p^v  and  4>vv  are  uniformly  distributed,  it  can  be  easily  shown  that 
the  cross- polarized  phase  difference  <f>c  =  <f>vh  —  <t>w  is  also  uniformly  distributed  between 

(-7T,+7T). 

The  co-polarized  elements  of  the  scattering  matrix,  however,  are  dependent  random 
variables  which  can  be  denoted  by  a  four-component  jointly  Gaussian  random  vector  X. 
Let  us  define 

Svv  =  X\  +  1X2  ,  Shh  —  %3  +  LX4 


and  since  X\,*"X4  are  Gaussian  their  joint  probability  density  function  can  be  fully 
determined  by  a  4  x  4  symmetric  positive  definite  matrix  known  as  covariance  matrix 


(A)  whose  entries  axe  given  by  [2] 


—  ^ji  “ ^  X{Xj  >  *»  J  €  { 1,  •  *  • ,  4} 

The  joint  probability  density  function  in  terms  of  the  covariance  matrix  takes  the  fol¬ 
lowing  form: 

/x(i,’  '’i')=^expHiA',x]  •  (4) 

where  X  is  transpose  of  the  column  vector  X.  To  characterize  the  covariance  matrix  the 
following  observations  are  in  order.  First,  it  was  shown  that  the  real  and  imaginary  parts 
of  the  scattering  matrix  elements  are  mutually  independent  and  identically  distributed 
zero-mean  random  variables,  therefore 


Au  =  A22  =<  Xj  >  =  <  X%  > 

A12  - <  X\Xi  >  =  0 
A33  =  A44  =<  Xl  >=<  X l  > 
A34  =<  X3X4  >=  0 


(5) 

(6) 

(7) 

(8) 


Second,  it  was  shown  that  the  absolute  phase  <j>pp  is  uniformly  distributed  and  is  inde¬ 
pendent  of  |5pP|.  Thus  the  random  variable  <j>vv  +  4>hh  is  also  uniformly  distributed  and 
is  independent  of  |Svv||Saa|  from  which  it  can  be  concluded  that 


<  ISuvIISm*!  COS^™  +  <}>hh)  >=  0 

<  I-S'i/wII-S'a/iI sin(0vv  +  <f>hh)  >=  0 


(9) 


In  fact,  the  complex  random  variable  SvvShh  is  obtained  from  a  similar  Wiener  process 


which  led  to  the  random  variables  Svv  and  Shh ■  On  the  other  hand 

X1X3  -  X2X4  =  |Sw||SfcA|  cos (<f>vv  +  <pHk)  , 
X1X4  +  X2X3  =  |Sw||$fcA|  sin(<£„„  +  4>h.h )  • 

In  view  of  (9)  and  (10)  it  can  easily  be  seen  that 

Al3  =  A24  , 

A14  =  — A23  . 


(10) 

(11) 

(12) 


The  properties  derived  for  the  entries  of  the  covariance  matrix,  as  given  by  (5)-(8)  and 
(11)-(12),  indicates  that  there  are  only  four  unknowns  left  in  the  covariance  matrix.  The 
unknowns,  namely  An,  Aj3,  An,  and  A33,  can  be  obtained  directly  from  the  Mueller 
matrix  of  the  target  as  will  be  shown  next. 

The  Mueller  matrix  relates  the  scattered-wave  Stokes  vector  to  the  incident-wave 
Stokes  vector  by  (van  Zyl  and  Ulaby,  1990] 

f*  =  -4mf'  , 

Tl 

where  E*’*  are  the  modified  incident-  and  scattered- wave  Stokes  vector  defined  by 

|£.IJ 

|£»l! 

20[£„£;] 

The  Mueller  matrix  can  be  expressed  in  terms  of  the  elements  of  the  scattering  matrix 
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as  follows  [6] 

|Sv„|2  \Svh\2  8[  S:hSvv)  -3[5v*a5vv) 

I^AvI2  |SaaI2  X[S'hhShv}  -*[SmhhShv] 

M  = 

23?[5v„5AuJ  23?[5’fcwSAA]  JJ[SUUSAA  +  ShvS*^]  ~'*[SwSkh  ~ 

20(5wu5Av]  2<3[5’At,.5'AA ]  9[5wt,5AA  +  -^av^va]  ^[^uv^aa  —  ^Av^vaI 

In  the  case  of  a  random  medium  we  are  dealing  with  a  partially  polarized  scattered  wave 

and  the  quantity  of  interest  is  the  ensemble  averaged  Mueller  matrix.  Using  the  assump¬ 
tion  that  the  co-  and  cross-polarized  terms  of  the  scattering  matrix  are  independent  and 
employing  the  properties  given  by  (5)-(8)  and  (ll)-(12),  the  Mueller  matrix  in  terms  of 
the  entries  of  the  covariance  matrix  is  given  by 

2AU  2  a2  0  0 

2<r?  2A33  0  0 

M  =<  M  >=  .  (13) 

0  0  2Aj3  ■(-  2<r 2  2Aj4 

0  0  -2A14  2A13  -  2a2 

Equation  (13)  provides  enough  equations  to  determine  the  unknown  elements  of  the 
covariance  matrix  and  variance  of  the  cross-polarized  component,  i.e. 


<72  - 

uc  — 

2  ’ 

** 

II 

jMii 

2  » 

A  33  = 

u 

II 

4  » 

An  = 

With  the  covariance  matrix,  the  joint  density  function  of  Xi ,  •  •  • ,  X4  can  be  obtained 
as  given  by  (4).  Using  a  rectangular  to  polar  transformation,  i.e. 

Xj  —  p\  cos  1  X2  —  p\  sin  (fiw  1 

X3  =  P?  COS  d>AA  >  X4  =  P2  sin  <t>hh  . 
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the  joint  probability  density  function  of  the  amplitudes  and  phases  takes  the  following 
form 


=  -£~=pip2exp  +  a2p\  -  2a3pi^2]|  ,  (14) 


where 

A=  |Aj  =  (A11A33  —  Aj3  -  Aj4)2  , 

<*1  =  A33/\/A  ,  a2  =  An/>/A  , 

=  [A13  cos (4>hh  -  4>vv)  +  Au  s\n(4>hh  ~  4>  vv  ))/V, A  . 

To  obtain  the  co-polarized  phase  difference  statistics  the  joint  density  function  of  <j>vv 

and  <f>hh  is  needed  which  can  be  obtained  from 


noo 

fpi  4>vv^4>hh)^P\^P2  •  (15) 

Noting  that  a\  is  a  positive  real  number,  the  integration  with  respect  to  pi  can  be  carried 
out  which  results  in 


4>hk)  — 


1  r  1  r  , 

vsuX  «e  ’  ^ 


4  w2y/K 
+ 


]fs 


l1±erf(v^'^) 


c—  a^7  (ai  <*2  -«1  )p: 


'l*P2 }  , 


(16) 


where  erf(-)  is  the  error  function  and  the  plus  or  minus  sign  is  used  according  to  the  sign 
of  03.  To  evaluate  the  integrals  in  (16),  we  need  to  show  that  both  a2  and  aia2  -  a3  are 
positive  numbers.  By  definition,  a2  is  positive  and  to  show  aja2  -  a3  is  positive  we  note 
that  A  is  a  symmetric  positive  definite  matrix,  therefore  its  eigenvalues  must  be  positive. 
It  can  be  shown  that  A  has  two  distinct  eigenvalues  71  and  72  each  with  multiplicity  2 
and  their  product  is  given  by 


a\ai  -  =  7172  +  (Ai3COs(<£m  —  <f>vv)  —  Ausin (<f>hh  -  4>w) ]2 


given  by 


An  = 


Ai3  = 


A^33  +  At  44 


A33  = 


A14  = 


At  34  -  A<43 


Some  limiting  cases  can  be  considered  in  order  to  check  the  validity  of  (17).  For  example, 
when  Svv  and  are  uncorrelated,  then  both  Al3  and  Ai4  are  zero  for  which  /$(<£)  = 
1/(2tt),  as  expected.  Also,  for  the  case  of  completely  polarized  scattered  wave  where  Svv 
and  Shh  are  completely  correlated,  the  determinant  of  A  is  zero  and  so  /$(<£)  is  a  delta 
function. 

It  is  interesting  to  note  that  the  p.d.f.  of  the  phase  difference  is  only  a  function  of 
two  parameters  defined  by 

_  /A13  +  Aj4  _  -1  A14 

V  *11*33  '  <_,i"  *,3 

where  a  and  (  can  vary  from  0  to  1  and  —7 r  to  tt  respectively.  In  fact  if  the  wave  were 
completely  polarized,  (  would  have  been  the  phase  difference  between  the  co-polarized 
terms.  The  parameter  (  will,  henceforth,  be  referred  to  as  the  polarized-phase-difference. 
In  terms  of  these  parameters  (17)  can  be  written  as 


M<f>)  - 


1  —  a2 

2 x  [1  —  a2  cos2(d>  —  C)] 

f  1  a  cos (<ft-Q  *  _i  acos(</>-C) 

\  \/\  —  a2  cos2(d>  —  C)  .2  y/l  —  a2  COS^(0  -  £) 


It  can  be  shown  that  the  maximum  of  the  probability  density  function  occurs  at  d>  = 
independent  of  a.  However,  the  width  of  the  p.d.f.  (e.g.  the  3  dB  angular  width)  is  only 
a  function  of  a  which  will  be  referred  to  as  the  degree  of  correlation.  The  probability 
distribution  function  given  by  (18)  is  the  analog  of  Gaussian  distribution  for  periodic 


random  variables  where  £  and  a  are,  respectively,  the  counterparts  of  the  mean  and 


12 


variance  for  Gaussian  random  variables.  Figure  1  shows  the  p.d.f.  for  different  values 
of  (“  while  keeping  a  constant  and  Fig.  2  shows  the  p.d.f.  for  a  fixed  vaJue  of  a 
while  changing  (  as  a  parameter.  The  calculated  mean  and  standard  deviation  of  the 
phase  difference  as  a  function  of  both  the  polarized-phase-difference  and  the  degree  of 
correlation  are  depicted  in  Figs.  3  and  4  respectively. 

Lastly,  it  is  necessary  to  point  out  that  the  formulation  of  the  co polarized  phase 
difference  p.d.f.,  as  given  in  (17),  is  not  restricted  to  the  backscattering  case  or  to  the 
co-  and  cross-polarized  components  being  uncorrelated.  In  fact  we  can  derive  the  cross- 
polarized  phase  difference  statistics  in  a  similar  manner  and  the  p.d.f.  in  this  case  for 
the  backscattering  case  can  be  obtained  from  (17)  upon  the  following  substitution  for 
the  elements  of  the  cross-polarized  covariance  matrix 

A.l  =  ^  ,  ^33=  ^  , 

Ai3=  ,  A,,  =  ^  . 

3  Comparison  with  Measurements 

Using  the  polarimetric  data  gathered  by  scatterometers  from  a  variety  of  natural  targets, 
the  assumptions  leading  to  the  probability  density  function  of  phase  differences  as  derived 
in  the  previous  section  are  examined.  Also  by  generating  the  histograms,  means,  and 
standard  deviations  of  the  phase  differences  from  the  data  and  comparing  them  with 
the  results  based  on  the  p.d.f.  derived  from  the  measured  Mueller  matrices  validity  of 
the  model  is  also  examined.  The  polarimetric  radar  measurements  of  bare  soil  surfaces 
were  performed  at  L-,  C-,  and  X-band  frequencies  for  a  total  of  eight  different  soil  surface 
conditions  (four  roughness  and  two  moisture  conditions).  For  this  experiment  we  tried  to 
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preserve  the  absolute  phase  of  the  measured  scattering  matrix  by  calibrating  the  surface 
data  with  a  metallic  sphere  located  at  the  same  distance  from  the  radar  as  the  center  of 
the  surface  target.  For  each  frequency,  surface  condition,  and  incidence  angle  a  minimum 
of  700  independent  samples  were  collected.  The  detailed  procedure  of  the  data  collection 
and  calibration  is  given  in  reference  [4]. 

By  generating  the  histograms  of  the  real  and  imaginary  parts  of  the  elements  of  the 
scattering  matrix  for  all  surfaces,  it  was  found  that  they  have  a  zero- mean  Gaussian 
distribution  as  we  assumed.  Figure  5  represents  a  typical  case  where  the  histogram  of 
the  real  and  imaginary  parts  of  Svv  and  S ^  of  a  dry  surface  with  rms  height  0.32  cm 
and  correlation  length  9.9  cm  at  C-band  have  a  bell-shaped  distribution.  The  properties 
of  the  covariance  matrix  as  given  by  (5)-(8)  and  (11)-(12)  are  verified  by  calculating  the 
covariance  matrices  of  the  data  for  all  cases.  Table  1  represents  a  typical  situation  where 
the  covariance  matrix  for  the  same  rough  surface  at  C-band  possesses  the  mentioned 
properties  approximately,  that  is  An  ^  A22,  A12  &  A34  «  0,  A33  ^  A44,  A13  %  A24,  and 
Am  «  — A23.  The  small  discrepancies  are  due  to  the  fact  that  the  measurement  of  the 
scattering  matrix  with  absolute  phase  has  an  uncertainty  of  ±30  degrees. 

Table  2  gives  the  Mueller  matrix  of  the  typical  surface  (Table  1)  at  C-band  from  which 
the  co-  and  cross- polarized  phase  difference  probability  density  functions  are  calculated 
using  (17)  and  are  compared  with  the  measured  phase  histograms  in  Figs.  6  and  7 
respectively.  Similar  comparisons  were  also  made  for  the  rest  of  surfaces,  frequencies, 
and  incidence  angles  and  it  was  found  that  the  expression  (17)  predicts  the  density 
functions  very  accurately.  Some  example  of  these  comparisons  are  shown  in  Figs.  8  and 
9.  Figures  8  and  9  compare  the  mean  and  standard  deviation  of  the  co-pularized  phase 
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difference  versus  incidence  angle  at  L-  and  X-band  for  a  surface  with  rms  height  0.4 
cm  and  correlation  length  8.4  cm  in  dry  condition  using  the  results  based  on  the  direct 
calculation  and  the  results  derived  from  (17). 

4  Conclusions 

Prompted  by  the  experimental  observations  which  show  strong  dependence  of  phase 
differences  of  the  scattering  matrix  elements  on  the  physical  parameters  of  random  media, 
the  statistical  behavior  of  the  phase  differences  for  distributed  targets  is  studied.  The 
probability  density  functions  of  the  phase  differences  are  derived  from  the  Mueller  matrix 
of  the  target.  In  derivation  of  the  density  functions  it  is  assumed  that  the  real  and 
imaginary  parts  of  the  co-  and  cross* polarized  terms  of  the  scattering  matrix  are  jointly 
Gaussian  and  their  covariance  matrices  are  found  in  terms  of  the  Mueller  matrix  elements. 
The  functional  form  of  the  co-  and  cross-polarized  density  functions  are  similar  and  are 
obtained  independently.  It  is  shown  that  the  density  function  of  the  phase  difference 
is  completely  determined  in  terms  of  only  two  parameters.  The  assumptions  and  final 
expressions  are  verified  by  using  a  set  of  polarimetric  data  acquired  by  scatterometers 
from  rough  surfaces. 
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1.00 

0.03 
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-0.12 

0.03 

0.90 

0.08 

0.68 

0.75 

0.08 

0.77 

0.05 

-0.12 

0.68 

0.05 

0.69 

Table  1:  Normalized  covariance  matrix  of  co-polarized  terms  of  scattering  matrix  for  a 
surface  with  rms  height  0.3  cm  and  correlation  length  9  cm  at  C-band  and  at  30  degrees 
incidence  angle. 


1.000 

0.028 

0.000 

0.000 


0.030 

0.767 

0.000 

0.000 


0.000 

0.000 

-0.11 

0.711 


0.000 

0.000 

0.770 

0.110 


Table  2:  Normalized  Mueller  matrix  for  a  surface  with  rms  height  0.3  cm  and  correlation 
length  9  cm  at  C-band  and  at  30  degrees  incidence  angle. 
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Figure  6:  The  histogram  and  p.d.f.  of  the  co-polarized  phase  difference  for  a  rough 
with  rms  height  0.32  cm  and  correlation  length  9.9  cm  at  C-band  and  30°  incident 
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Figure  7:  The  histogram  and  p.d.f.  of  the  cross-polarized  phase  difference  for  a  rough 
surface  with  rms  height  0.32  cm  and  correlation  length  9.9  cm  at  C-band  and  30°  inci¬ 
dence  angle. 


23 


Mean  of  d>hh  -  <t>  (Degrees) 


Figure  8:  Angular  dependency  of  the  mean  of  the  co-polarized  phase  difference  for  a  dry 
rough  surface  with  rms  height  0.4  cm  and  correlation  length  8.4  cm  at  L-  and  X-band. 
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Figure  9:  Angular  dependency  of  the  standard  deviation  of  the  co-polarized  phase  dif¬ 
ference  for  a  dry  rough  surface  with  rms  height  0.4  cm  and  correlation  length  8.4  cm  at 
L-  and  X-band. 
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